Genome-wide transcriptomic and proteomic analyses of bollworm-infested developing cotton bolls revealed the genes and pathways involved in the insect pest defence mechanism

Genome-Wide Transcriptomic and Metabolomic Analyses Unveiling the Defence Mechanisms of Populus tremula against Sucking and Chewing Insect Herbivores

Plants and insects coevolved as an evolutionarily successful and enduring association. The molecular arms race led to evolutionary novelties regarding unique mechanisms of defence and detoxification in plants and insects. While insects adopt mechanisms to conquer host defence, trees develop well-orchestrated and species-specific defence strategies against insect herbivory. However, current knowledge on the molecular underpinnings of fine-tuned tree defence responses against different herbivore insects is still restricted. In the current study, using a multi-omics approach, we unveiled the defence response of Populus tremula against aphids (Chaitophorus populialbae) and spongy moths (Lymantria dispar) herbivory. Comparative differential gene expression (DGE) analyses revealed that around 272 and 1203 transcripts were differentially regulated in P. tremula after moth and aphid herbivory compared to uninfested controls. Interestingly, 5716 transcripts were differentially regulated in P. tremula between aphids and moth infestation. Further investigation showed that defence-related stress hormones and their lipid precursors, transcription factors, and signalling molecules were over-expressed, whereas the growth-related counterparts were suppressed in P. tremula after aphid and moth herbivory. Metabolomics analysis documented that around 37% of all significantly abundant metabolites were associated with biochemical pathways related to tree growth and defence. However, the metabolic profiles of aphid and moth-fed trees were quite distinct, indicating species-specific response optimization. After identifying the suitable reference genes in P. tremula, the omics data were further validated using RT-qPCR. Nevertheless, our findings documented species-specific fine-tuning of the defence response of P. tremula, showing conservation on resource allocation for defence overgrowth under aphid and moth herbivory. Such findings can be exploited to enhance our current understanding of molecular orchestration of tree responses against herbivory and aid in developing insect pest resistance P. tremula varieties.

A comprehensive review on Gossypium hirsutum resistance against cotton leaf curl virus

Cotton (Gossypium hirsutum L.) is a significant fiber crop. Being a major contributor to the textile industry requires continuous care and attention. Cotton is subjected to various biotic and abiotic constraints. Among these, biotic factors including cotton leaf curl virus (CLCuV) are dominant. CLCuV is a notorious disease of cotton and is acquired, carried, and transmitted by the whitefly (Bemisia tabaci). A cotton plant affected with CLCuV may show a wide range of symptoms such as yellowing of leaves, thickening of veins, upward or downward curling, formation of enations, and stunted growth. Though there are many efforts to protect the crop from CLCuV, long-term results are not yet obtained as CLCuV strains are capable of mutating and overcoming plant resistance. However, systemic-induced resistance using a gene-based approach remained effective until new virulent strains of CLCuV (like Cotton Leaf Curl Burewala Virus and others) came into existence. Disease control by biological means and the development of CLCuV-resistant cotton varieties are in progress. In this review, we first discussed in detail the evolution of cotton and CLCuV strains, the transmission mechanism of CLCuV, the genetic architecture of CLCuV vectors, and the use of pathogen and nonpathogen-based approaches to control CLCuD. Next, we delineate the uses of cutting-edge technologies like genome editing (with a special focus on CRISPR-Cas), next-generation technologies, and their application in cotton genomics and speed breeding to develop CLCuD resistant cotton germplasm in a short time. Finally, we delve into the current obstacles related to cotton genome editing and explore forthcoming pathways for enhancing precision in genome editing through the utilization of advanced genome editing technologies. These endeavors aim to enhance cotton's resilience against CLCuD.

An insight into the gene expression evolution in Gossypium species based on the leaf transcriptomes

BACKGROUND

Gene expression pattern is associated with biological phenotype and is widely used in exploring gene functions. Its evolution is also crucial in understanding species speciation and divergence. The genus Gossypium is a bona fide model for studying plant evolution and polyploidization. However, the evolution of gene expression during cotton species divergence has yet to be extensively discussed.

RESULTS

Based on the seedling leaf transcriptomes, this work analyzed the transcriptomic content and expression patterns across eight cotton species, including six diploids and two natural tetraploids. Our findings indicate that, while the biological function of these cotton transcriptomes remains largely conserved, there has been significant variation in transcriptomic content during species divergence. Furthermore, we conducted a comprehensive analysis of expression distances across cotton species. This analysis lends further support to the use of G. arboreum as a substitute for the A-genome donor of natural cotton polyploids. Moreover, our research highlights the evolution of stress-responsive pathways, including hormone signaling, fatty acid degradation, and flavonoid biosynthesis. These processes appear to have evolved under lower selection pressures, presumably reflecting their critical role in the adaptations of the studied cotton species to diverse environments.

CONCLUSIONS

In summary, this study provided insights into the gene expression variation within the genus Gossypium and identified essential genes/pathways whose expression evolution was closely associated with the evolution of cotton species. Furthermore, the method of characterizing genes and pathways under unexpected high or slow selection pressure can also serve as a new strategy for gene function exploration.

Pathogenesis-related protein 10 in resistance to biotic stress: progress in elucidating functions, regulation and modes of action

Introduction

The Family of pathogenesis-related proteins 10 (PR-10) is widely distributed in the plant kingdom. PR-10 are multifunctional proteins, constitutively expressed in all plant tissues, playing a role in growth and development or being induced in stress situations. Several studies have investigated the preponderant role of PR-10 in plant defense against biotic stresses; however, little is known about the mechanisms of action of these proteins. This is the first systematic review conducted to gather information on the subject and to reveal the possible mechanisms of action that PR-10 perform.

Methods

Therefore, three databases were used for the article search: PubMed, Web of Science, and Scopus. To avoid bias, a protocol with inclusion and exclusion criteria was prepared. In total, 216 articles related to the proposed objective of this study were selected.

Results

The participation of PR-10 was revealed in the plant's defense against several stressor agents such as viruses, bacteria, fungi, oomycetes, nematodes and insects, and studies involving fungi and bacteria were predominant in the selected articles. Studies with combined techniques showed a compilation of relevant information about PR-10 in biotic stress that collaborate with the understanding of the mechanisms of action of these molecules. The up-regulation of PR-10 was predominant under different conditions of biotic stress, in addition to being more expressive in resistant varieties both at the transcriptional and translational level.

Discussion

Biological models that have been proposed reveal an intrinsic network of molecular interactions involving the modes of action of PR-10. These include hormonal pathways, transcription factors, physical interactions with effector proteins or pattern recognition receptors and other molecules involved with the plant's defense system.

Conclusion

The molecular networks involving PR-10 reveal how the plant's defense response is mediated, either to trigger susceptibility or, based on data systematized in this review, more frequently, to have plant resistance to the disease.

Gene Profiling of the Ascorbate Oxidase Family Genes under Osmotic and Cold Stress Reveals the Role of AnAO5 in Cold Adaptation in Ammopiptanthus nanus

The uplift of the Qinghai Tibet Plateau has led to a drastic change in the climate in Central Asia, from warm and rainy, to dry and less rainfall. Ammopiptanthus nanus, a rare evergreen broad-leaved shrub distributed in the temperate desert region of Central Asia, has survived the drastic climate change in Central Asia caused by the uplift of the Qinghai-Tibet Plateau. Ascorbate oxidase (AO) regulates the redox status of the apoplast by catalyzing the oxidation of ascorbate acid to dehydroascorbic acid, and plays a key role in the adaptation of plants to environmental changes. Analyzing the evolution, environmental response, and biological functions of the AO family of A. nanus is helpful for understanding how plant genome evolution responds to climate change in Central Asia. A total of 16 AOs were identified in A. nanus, all of which contained the ascorbate oxidase domain, most of which contained transmembrane domain, and many were predicted to be localized in the apoplast. Segmental duplication and tandem duplication are the main factors driving the gene amplification of the AO gene family in A. nanus. Gene expression analysis based on transcriptome data and fluorescence quantitative PCR, as well as enzyme activity measurements, showed that the expression levels of AO genes and total enzyme activity decreased under short-term osmotic stress and low-temperature stress, but the expression of some AO genes (AnAO5, AnAO13, and AnAO16) and total enzyme activity increased under 7 days of cold stress. AnAO5 and AnAO11 are targeted by miR4415. Further functional studies on AnAO5 showed that AnAO5 protein was localized in the apoplast. The expression of AnAO5 in yeast cells and the transient expression in tobacco enhanced the tolerance of yeast and tobacco to low-temperature stress, and the overexpression of AnAO5 enhanced the tolerance of Arabidopsis seedlings to cold stress. Our research provides important data for understanding the role of AOs in plant adaptation to environmental change.

Screening of Key Indices and the Gene Transcriptional Regulation Analysis Related to Salt Tolerance in Salix matsudana Seedlings

Pot experiments were performed to comparatively study the differences in 16 salt tolerance indices between the seedlings of six Salix matsudana clones under the stress of various concentrations of NaCl (0, 0.1%, 0.3%, 0.5%, and 0.7%), including the salt injury index, shoot fresh weight, root fresh weight, leaf water content, relative conductivity, malondialdehyde content, and antioxidant enzyme activity. The salt-tolerant clones and key indices of salt tolerance were selected. Transcriptome sequencing analysis was performed on the selected salt-tolerant and salt-sensitive clones under salt stress, and the links between the physiological indices of salt tolerance and gene expression were analyzed. Results: (1) Superoxide dismutase (SOD), peroxidase (POD), chlorophyll, and net photosynthetic rate were closely related to the salt tolerance of Salix matsudana at the seedling stage. The regression equation was constructed as follows: salt tolerance index (y) = 0.224x10 + 0.216x11 + 0.127x12 + 0.191x7 − 0.187 (x10 = chlorophyll, x11 = SOD, x12 = POD, x7 = net photosynthetic rate). (2) The number of differentially expressed genes between the seedlings of salt-tolerant and salt-sensitive clones varied with the time of exposure (0 h, 4 h, 12 h, and 24 h) to 200 mmol·L−1 NaCl stress. The most differentially expressed genes in Sm172 were detected upon 24 h vs. 4 h of salt treatment, while the most in Sm6 were in the 24 h vs. 0 h comparison. Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes analysis showed that several differentially expressed genes were involved in carotenoid biosynthesis and plant mitogen-activated protein kinase signaling pathways. The nine highly expressed transcription factor genes (Sm172-f2p30-2392, Sm172-f2p28-2386, Sm6-f8p60-2372, Sm6-f2p39-2263, Sm6-f16p60-2374, Sm6-f3p60-931, Sm6-f2p60-1067, Sm172-f3p54-1980, and Sm172-f3p54-1980) were closely correlated with the four key indices of salt tolerance. These genes could become genetic resources for salt tolerance breeding of Salix matsudana.

Global gene expression profiling under nitrogen stress identifies key genes involved in nitrogen stress adaptation in maize (Zea mays L.)

Maize is a heavy consumer of fertilizer nitrogen (N) which not only results in the high cost of cultivation but may also lead to environmental pollution. Therefore, there is a need to develop N-use efficient genotypes, a prerequisite for which is a greater understanding of N-deficiency stress adaptation. In this study, comparative transcriptome analysis was performed using leaf and root tissues from contrasting inbred lines, viz., DMI 56 (tolerant to N stress) and DMI 81 (susceptible to N stress) to delineate the differentially expressed genes (DEGs) under low-N stress. The contrasting lines were grown hydroponically in modified Hoagland solution having either sufficient- or deficient-N, followed by high-throughput RNA-sequencing. A total of 8 sequencing libraries were prepared and 88–97% of the sequenced raw reads were mapped to the reference B73 maize genome. Genes with a p value ≤ 0.05 and fold change of ≥ 2.0 or ≤ − 2 were considered as DEGs in various combinations performed between susceptible and tolerant genotypes. DEGs were further classified into different functional categories and pathways according to their putative functions. Gene Ontology based annotation of these DEGs identified three different functional categories: biological processes, molecular function, and cellular component. The KEGG and Mapman based analysis revealed that most of the DEGs fall into various metabolic pathways, biosynthesis of secondary metabolites, signal transduction, amino acid metabolism, N-assimilation and metabolism, and starch metabolism. Some of the key genes involved in N uptake (high-affinity nitrate transporter 2.2 and 2.5), N assimilation and metabolism (glutamine synthetase, asparagine synthetase), redox homeostasis (SOD, POX), and transcription factors (MYB36, AP2-EREBP) were found to be highly expressed in the tolerant genotype compared to susceptible one. The candidate genes identified in the present study might be playing a pivotal role in low-N stress adaptation in maize and hence could be useful in augmenting further research on N metabolism and development of N-deficiency tolerant maize cultivars.

Integration of transcriptome and proteome reveals molecular mechanisms underlying stress responses of the cutworm, Spodoptera litura, exposed to different levels of lead (Pb)

Heavy metals are major environmental pollutants that affect organisms across different trophic levels. Herbivorous insects play an important role in the bioaccumulation, and eventually, biomagnification of these metals. Although effects of heavy metal stress on insects have been well-studied, the molecular mechanisms underlying their effects remain poorly understood. Here, we used the RNA-Seq profiling and isobaric tags for relative and absolute quantitation (iTRAQ) approaches to unravel these mechanisms in the polyphagous pest Spodoptera litura exposed to lead (Pb) at two different concentrations (12.5 and 100 mg Pb/kg; PbL and PbH, respectively). Altogether, 1392 and 1630 differentially expressed genes (DEGs) and 58, 114 differentially expressed proteins (DEPs) were identified in larvae exposed to PbL and PbH, respectively. After exposed to PbL, the main up-regulated genes clusters and proteins in S. litura larvae were associated with their metabolic processes, including carbohydrate, protein, and lipid metabolism, but the levels of cytochrome P450 associated with the pathway of xenobiotic biodegradation and metabolism were found to be decreased. In contrast, the main up-regulated genes clusters and proteins in larvae exposed to PbH were enriched in the metabolism of xenobiotic by cytochrome P450, drug metabolism-cytochrome P450, and other drug metabolism enzymes, while the down-regulated genes and proteins were found to be closely related to the lipid (lipase) and protein (serine protease, trypsin) metabolism and growth processes (cuticular protein). These findings indicate that S. litura larvae exposed to PbL could enhance food digestion and absorption to prioritize for growth rather than detoxification, whereas S. litura larvae exposed to PbH reduced food digestion and absorption and channelized the limited energy for detoxification rather than growth. These contrasting results explain the dose-dependent effects of heavy metal stress on insect life-history traits, wherein low levels of heavy metal stress induce stimulation, while high levels of heavy metal stress cause inhibition at the transcriptome and proteome levels.

Transcriptome analysis reveals genes potentially related to high fiber strength in a Gossypium hirsutum line IL9 with Gossypium mustelinum introgression

Cotton (Gossypium L.) is the most important fiber crop worldwide. Here, transcriptome analysis was conducted on developing fibers of a G. mustelinum introgression line, IL9, and its recurrent parent, PD94042, at 17 and 21 days post-anthesis (dpa). Differentially expressed genes (DEGs) of PD94042 and IL9 were identified. Gene Ontology (GO) enrichment analysis showed that the annotated DEGs were rich in two main biological processes and two main molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis likewise showed that the annotated DEGs were mainly enriched in metabolic pathways and biosynthesis of secondary metabolites. In total, 52 DEGs were selected as candidate genes based on comparison of the DEGs and GO function annotation information. Quantitative real-time PCR (RT-qPCR) analysis results for 12 randomly selected DEGs were consistent with transcriptome analysis. SNP identification based on G. mustelinum chromatin segment introgression showed that 394 SNPs were identified in 268 DEGs, and two genes with known functions were identified within fiber strength quantitative trait loci (QTL) regions or near the confidence intervals. We identified 52 key genes potentially related to high fiber strength in a G. mustelinum introgression line and provided significant insights into the study of cotton fiber quality improvement.

Label-free quantitative proteomics of Sorghum bicolor reveals the proteins strengthening plant defense against insect pest Chilo partellus

BACKGROUND

Spotted stem borer- Chilo partellus - a Lepidopteran insect pest of Sorghum bicolor is responsible for major economic losses. It is an oligophagous pest, which bores through the plant stem, causing 'deadheart' and hampering the development of the main cob. We applied a label-free quantitative proteomics approach on three genotypes of S. bicolor with differential resistance/ susceptibility to insect pests, intending to identify the S. bicolor's systemic protein complement contributing to C. partellus tolerance.

METHODS

The proteomes of S. bicolor with variable resistance to insect pests, ICSV700, IS2205 (resistant) and Swarna (susceptible) were investigated and compared using label-free quantitative proteomics to identify putative leaf proteins contributing to resistance to C. partellus.

RESULTS

The multivariate analysis on a total of 967 proteins led to the identification of proteins correlating with insect resistance/susceptibility of S. bicolor. Upon C. partellus infestation S. bicolor responded by suppression of protein and amino acid biosynthesis, and induction of proteins involved in maintaining photosynthesis and responding to stresses. The gene ontology analysis revealed that C. partellus-responsive proteins in resistant S. bicolor genotypes were mainly involved in stress and defense, small molecule biosynthesis, amino acid metabolism, catalytic and translation regulation activities. At steady-state, the resistant S. bicolor genotypes displayed at least two-fold higher numbers of unique proteins than the susceptible genotype Swarna, mostly involved in catalytic activities. Gene expression analysis of selected candidates was performed on S. bicolor by artificial induction to mimic C. partellus infestation.

CONCLUSION

The collection of identified proteins differentially expressed in resistant S. bicolor, are interesting candidates for further elucidation of their role in defense against insect pests.

Genetic and evolution analysis of extrafloral nectary in cotton

Extrafloral nectaries are a defence trait that plays important roles in plant-animal interactions. Gossypium species are characterized by cellular grooves in leaf midribs that secret large amounts of nectar. Here, with a panel of 215 G. arboreum accessions, we compared extrafloral nectaries to nectariless accessions to identify a region of Chr12 that showed strong differentiation and overlapped with signals from GWAS of nectaries. Fine mapping of an F2 population identified GaNEC1, encoding a PB1 domain-containing protein, as a positive regulator of nectary formation. An InDel, encoding a five amino acid deletion, together with a nonsynonymous substitution, was predicted to cause 3D structural changes in GaNEC1 protein that could confer the nectariless phenotype. mRNA-Seq analysis showed that JA-related genes are up-regulated and cell wall-related genes are down-regulated in the nectary. Silencing of GaNEC1 led to a smaller size of foliar nectary phenotype. Metabolomics analysis identified more than 400 metabolites in nectar, including expected saccharides and amino acids. The identification of GaNEC1 helps establish the network regulating nectary formation and nectar secretion, and has implications for understanding the production of secondary metabolites in nectar. Our results will deepen our understanding of plant-mutualism co-evolution and interactions, and will enable utilization of a plant defence trait in cotton breeding efforts.

Constitutive expression of Asparaginase in Gossypium hirsutum triggers insecticidal activity against Bemisia tabaci

Whitefly infestation of cotton crop imparts enormous damage to cotton yield by severely affecting plant health, vigour and transmitting Cotton Leaf Curl Virus (CLCuV). Genetic modification of cotton helps to overcome both the direct whitefly infestation as well as CLCuV based cotton yield losses. We have constitutively overexpressed asparaginase (ZmASN) gene in Gossypium hirsutum to overcome the cotton yield losses imparted by whitefly infestation. We achieved 2.54% transformation efficiency in CIM-482 by Agrobacterium-mediated shoot apex transformation method. The relative qRT-PCR revealed 40-fold higher transcripts of asparaginase in transgenic cotton line vs. non-transgenic cotton lines. Metabolic analysis showed higher contents of aspartic acid and glutamic acid in seeds and phloem sap of the transgenic cotton lines. Phenotypically, the transgenic cotton lines showed vigorous growth and height, greater number of bolls, and yield. Among six representative transgenic cotton lines, line 14 had higher photosynthetic rate, stomatal conductance, smooth fiber surface, increased fiber convolutions (SEM analysis) and 95% whitefly mortality as compared to non-transgenic cotton line. The gene integration analysis by fluorescence in situ hybridization showed single copy gene integration at chromosome number 1. Collectively, asparaginase gene demonstrated potential to control whitefly infestation, post-infestation damages and improve cotton plant health and yield: a pre-requisite for farmer's community.

RNA-Seq analysis reveals transcript diversity and active genes after common cutworm (Spodoptera litura Fabricius) attack in resistant and susceptible wild soybean lines

BACKGROUND

Common cutworm (CCW) is highly responsible for destabilizing soybean productivity. Wild soybean is a resource used by breeders to discover elite defensive genes.

RESULTS

The transcriptomes of two wild accessions (W11 and W99) with different resistance to CCW were analyzed at early- and late-induction time points. After induction, the susceptible accession W11 differentially expressed 1268 and 508 genes at the early and late time points, respectively. Compared with W11, the resistant accession W99 differentially expressed 1270 genes at the early time point and many more genes (2308) at the late time point. In total, 3836 non-redundant genes were identified in both lines. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the differentially expressed genes (DEGs) in W99 at the late time point were mostly associated with specific processes and pathways. Among the non-redundant genes, 146 genes were commonly up-regulated in the treatment condition compared with the control condition at the early- and late-induction time points in both accessions used in this experiment. Approximately 40% of the common DEGs were related to secondary metabolism, disease resistance, and signal transduction based on their putative function. Excluding the common DEGs, W99 expressed more unique DEGs than W11. Further analysis of the 3836 DEGs revealed that the induction of CCW not only up-regulated defense-related genes, including 37 jasmonic acid (JA)-related genes, 171 plant-pathogen-related genes, and 17 genes encoding protease inhibitors, but also down-regulated growth-related genes, including 35 photosynthesis-related genes, 48 nutrition metabolism genes, and 28 auxin metabolism genes. Therefore, representative defense-related and growth-related genes were chosen for binding site prediction via co-expression of transcription factors (TFs) and spatial expression pattern analyses. In total, 53 binding sites of 28 TFs were identified based on 3 defense-related genes and 3 growth-related genes. Phosphate transporter PT1, which is a representative growth-related gene, was transformed into soybean, and the transgenic soybean plants were susceptible to CCW.

CONCLUSIONS

In summary, we described transcriptome reprograming after herbivore induction in wild soybean, identified the susceptibility of growth-related genes, and provided new resources for the breeding of herbivore-resistant cultivated soybeans.

The early transcriptome response of cassava (Manihot esculenta Crantz) to mealybug (Phenacoccus manihoti) feeding

The mealybug, Phenacoccus manihoti, is a leading pest of cassava (Manihot esculenta Crantz), damaging this crop globally. Although the biological control of this mealybug using natural predators has been established, resistance breeding remains an important means of control. Understanding plant responses to insect herbivory, by determining and identifying differentially expressed genes (DEGs), is a vital step towards the understanding of molecular mechanisms of defence responses in plants and the development of resistant cultivars by gene editing. Morphological and molecular analysis confirmed the mealybug identity as Phenacoccus manihoti (Matile-Ferrero). The transcriptome response of the green mite resistant cassava genotype AR23.1 was compared to P40/1 with no known resistance at 24 and 72 hours of mealybug infestation compared to non-infested mock. A total of 301 and 206 genes were differentially expressed at 24 and 72 of mealybug infestation for AR23.1 and P40/1 genotypes respectively, using a log2 fold change and P-value ≤ 0.05. Gene ontology functional classification revealed an enrichment of genes in the secondary metabolic process category in AR23.1 in comparison with P40/1, while genes in the regulation of molecular function, cellular component biogenesis and electron carrier categories were more significantly enriched in P40/1 than in AR23.1. Biological pathway analysis, based on KEGG, revealed a significant enrichment of plant-pathogen interaction and plant hormonal signal transduction pathways for a cohort of up-regulated and down-regulated DEGs in both genotypes. Defence-related genes such as 2-oxogluterate, gibberellin oxidase and terpene synthase proteins were only induced in genotype AR23.1 and not in P40/1, and subsequently validated by RT-qPCR. The study revealed a difference in response to mealybug infestation in the two genotypes studied, with AR23.1 showing a higher number of differentially expressed transcripts post mealybug infestation at 24 and 72 hours. Candidate defence-related genes that were overexpressed in the AR23.1 genotype post mealybug infestation will be useful in future functional studies towards the control of mealybugs.