Next Article in Journal
Chemical Composition of the Cinnamomum malabatrum Leaf Essential Oil and Analysis of Its Antioxidant, Enzyme Inhibitory and Antibacterial Activities
Next Article in Special Issue
Evaluating the Translational Potential of Bacteriocins as an Alternative Treatment for Staphylococcus aureus Infections in Animals and Humans
Previous Article in Journal
Outcome of Irrigation and Debridement with Topical Antibiotic Delivery Using Antibiotic-Impregnated Calcium Hydroxyapatite for the Management of Periprosthetic Hip Joint Infection
Previous Article in Special Issue
Antimicrobial and Antibiofilm Effect of Bacteriocin-Producing Pediococcus inopinatus K35 Isolated from Kimchi against Multidrug-Resistant Pseudomonas aeruginosa
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 12
Issue 5
10.3390/antibiotics12050939
antibiotics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

From Gene to Transcript and Peptide: A Deep Overview on Non-Specific Lipid Transfer Proteins (nsLTPs)

by
Carlos André dos Santos-Silva
1,†,
José Ribamar Costa Ferreira-Neto
2,†,
Vinícius Costa Amador
2,
João Pacífico Bezerra-Neto
3,
Lívia Maria Batista Vilela
2,
Eliseu Binneck
4,
Mireli de Santana Rêgo
2,
Manassés Daniel da Silva
2,
Ana Luiza Trajano Mangueira de Melo
2,
Rahisa Helena da Silva
2 and
Ana Maria Benko-Iseppon
2,*
1
Institute for Maternal and Child Health–IRCCS Burlo Garofolo, 34137 Trieste, Italy
2
Departamento de Genética, Centro de Biociências, Universidade Federal de Pernambuco, Recife 50670-901, Brazil
3
Instituto de Ciências Biológicas, Universidade de Pernambuco, Recife 50100-010, Brazil
4
Empresa Brasileira de Pesquisa Agropecuária, Embrapa Soja, Londrina 86085-981, Brazil
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Antibiotics 2023, 12(5), 939; https://doi.org/10.3390/antibiotics12050939
Submission received: 28 April 2023 / Revised: 13 May 2023 / Accepted: 18 May 2023 / Published: 21 May 2023
(This article belongs to the Special Issue Alternatives to Antibiotics: Bacteriocins and Antimicrobial Peptides, 2nd Edition)
Browse Figures
Review Reports Versions Notes

Abstract

:
Non-specific lipid transfer proteins (nsLTPs) stand out among plant-specific peptide superfamilies due to their multifaceted roles in plant molecular physiology and development, including their protective functions against pathogens. These antimicrobial agents have demonstrated remarkable efficacy against bacterial and fungal pathogens. The discovery of plant-originated, cysteine-rich antimicrobial peptides such as nsLTPs has paved the way for exploring the mentioned organisms as potential biofactories for synthesizing antimicrobial compounds. Recently, nsLTPs have been the focus of a plethora of research and reviews, providing a functional overview of their potential activity. The present work compiles relevant information on nsLTP omics and evolution, and it adds meta-analysis of nsLTPs, including: (1) genome-wide mining in 12 plant genomes not studied before; (2) latest common ancestor analysis (LCA) and expansion mechanisms; (3) structural proteomics, scrutinizing nsLTPs’ three-dimensional structure/physicochemical characteristics in the context of nsLTP classification; and (4) broad nsLTP spatiotemporal transcriptional analysis using soybean as a study case. Combining a critical review with original results, we aim to integrate high-quality information in a single source to clarify unexplored aspects of this important gene/peptide family.

1. Introduction

The non-specific lipid transfer proteins (nsLTPs) are a plant-specific superfamily of cysteine-rich AMPs (antimicrobial peptides). They received this name due to their ability to bind to several hydrophobic molecules, such as phospholipids and fatty acids, among others. nsLTPs are characterized by their reduced size (6.5–10.5 kDa) and the presence of eight cysteine residues (8CM domain), which form four disulfide bonds [1]. They are associated with various plant biological processes, such as growth and development, abiotic stress responses, besides plant defense [2,3,4,5,6,7,8]. The mentioned plant-originated AMPs have demonstrated remarkable efficacy against bacterial and fungal pathogens [9]. The discovery of cysteine-rich AMPs such as nsLTPs has paved the way for exploring plants as potential biofactories for synthesizing antimicrobial compounds, which holds significant promise for their application as biotherapeutic agents in the field of antimicrobial drug development [9]. It has been proposed that such antimicrobial activity is due to the nsLTPs’ ability to disrupt the permeability and integrity of the pathogens’ outer membranes, similar to other plant AMPs [7,10]. However, further studies are necessary to understand all their biological roles.
Previous works indicate that nsLTPs are encoded by a large gene family, presenting more than 50 loci in many angiosperm genomes and up to 50 loci in bryophytes, ferns, and gymnosperms [1]. Some classification systems are available for nsLTPs. These, however, are heterogeneous in terms of subgroups’ numbers and nomenclature. The classification initially proposed for nsLTPs—division into ‘nsLTP1’ and ‘nsLTP2’ subfamilies—was based on their molecular mass, sequence identity (<30% of similarity between nsLTP1s and nsLTP2s), and lipid transfer efficiency [11,12]. Disulfide bond patterns in nsLTPs may also differ. ‘nsLTP1s’ displays ‘Cys1-Cys6, Cys2-Cys3, Cys4-Cys7, and Cys5-Cys8’ pattern [13], whereas ‘nsLTP2s’ exhibit that of ‘Cys1-Cys5, Cys2-Cys3, Cys4-Cys7, and Cys6-Cys8’ [14]. Such differences in disulfide bond configuration play a crucial role in the stability of these proteins, constraining their conformational dynamics [1].
Another notable difference regards the hydrophobic cavity. The subgroup nsLTP1s may have a long tunnel-like cavity [15,16], while nsLTP2s may have two adjacent hydrophobic cavities [14]. However, it is worth mentioning that there is no specific rule in this regard since the hydrophobic cavities can vary according to the sequence of amino acid residues, the number of disulfide bonds present in the structure, and the lipid binding specificity.
Boutrot et al. [17] proposed a different classification approach. These authors introduced phylogenetic grouping as a key classification criterion. The proposition includes early diverging nsLTP homologs found in mosses and liver plants. The mentioned classification established nine groups (types ‘1’ to ‘9’) stratified with their respective consensus cysteine motifs (8CM domain)—Cys-Xn-Cys-Xn-CysCys-Xn-CysXCys-Xn-Cys-Xn-Cys (X: different amino acids; n: variable number of amino acids)—and the inter-cysteine amino acid residues diversity.
Later, another proposition, by Edstam et al. [1], sought to break the limitations inherent in the sequence conservation, considering post-translational modifications in the glycosylphosphatidylinositol (GPI) anchoring sites, intron positions, and spacing in the regions between cysteine residues, besides sequence similarity. This classification system suggests five major groups (‘LTP1’, ‘LTP2’, ‘LTPc’, ‘LTPd’, and ‘LTPg’) and four minor groups (‘LTPe’, ‘LTPf’, ‘LTPh’, ‘LTPj’, and ‘LTPk’). Despite the efforts made on the new nsLTP classification systems, the conventional classification of ‘LTP1’ and ‘LTP2’ is still widely used due to a lack of consensus among different studies. Since the nsLTP gene family is so complex and diversified, no established classification guidelines are final [18].
Most publications on nsLTP studies focus on their structure, abundance, and diversity in well-studied plant clades or model organisms. In the present work, we focused on the nsLTPs’ genome-wide annotation in some plant genomes with no nsLTP data or new genome versions available. Additionally, this work covers relevant information about the referred peptide group. It is divided into four sections, with the first three covering omics studies (genomics and transcriptomics), evolution (gene expansion mechanisms and lowest common ancestor analysis), and structural proteomics (focusing on three-dimensional structure/physicochemical characteristics) in the context of their classification.
We also present a case study for a comprehensive transcriptional analysis in soybean, considering both baseline (covering different soybean tissues and developmental stages) and differential gene expressions under diverse biotic and abiotic stresses. The present work, thus, provides a critical review, together with a meta-analysis and original results, aiming to synthesize, in a single source, high-quality information about this important peptide group.

2. Omics Studies for nsLTPs

2.1. Understanding nsLTPs from Previous Studies

Previous studies (Table 1) involving nsLTP gene search and plant transcriptional expression analyses were scrutinized (mining methodology in Supplementary Material S1). Based on this search, different views on nsLTP abundance, types (Table 1), and functions were found. The seminal work by Edstam et al. [1] proposes a greater nsLTP abundance in terrestrial plants and their absence in green algae (chlorophytes and charophytes) (Table 1), suggesting that nsLTP genes evolved soon after the terrestrial environment conquest. In favor of the mentioned proposition, a limited number of representatives and types of nsLTPs are observed when comparing lower plants (as bryophytes and lichens) to spermatophytes, which indicates the emergence of new types of nsLTPs in higher plants [1,19,20].
In the barley (Hordeum vulgare) genome, 70 HvnsLTPs (Hordeum vulgare nsLTPs) were identified (Table 1), which were classified into five groups (‘1’, ‘2’, ‘C’, ‘D’, and ‘G’) [27]. Each of these genes shared common structures. Considering their expansion mechanisms, the 70 HvnsLTP genes presented 15 tandem duplication repeats (encompassing 36 genes). The HvnsLTPs’ baseline expression profiles in different tissues across developmental stages indicated that this group of genes might perform a variety of functions [27]. In addition, the differential expression profile indicated that HvnsLTP genes might have diverged in terms of the cis-regulatory elements of their promoters [27].
For the Solanaceae family, there are data on nsLTPs for potato (StnsLTPs, Solanum tuberosum nsLTPs). Li et al. [29] found 83 StnsLTP genes in potato genomes, categorized into eight types, namely, ‘1’, ‘2’, ‘4’, ‘5’, ‘7’, ‘8’, ‘12*’, and ‘13*’ (Table 1). Chromosome distribution and collinearity analyses suggested that the expansion of the StnsLTP gene family was enhanced by tandem duplications. In turn, Ka/Ks analysis showed that 47 pairs of duplicated genes have gone through purifying selection during evolution. StnsLTP genes were expressed mainly in younger tissues. Furthermore, StnsLTPs contained a large number of stress-responsive, cis-acting elements in their promoter regions. These results indicated that StnsLTPs might play significant and functionally varied roles in potato plants.
In Arachis duranensis, Song et al. [31] discovered 64 AdnsLTPs (Arachis duranensis nsLTPs) genes, which were divided into six groups (‘1’, ‘2’, ‘C’, ‘D’, ‘E’, and ‘G’; Table 1), anchored over nine chromosomes. Considering the AdnsLTPs’ expansion mechanisms, the study revealed some gene clustering by tandem duplication, while other family members showed segmental duplication in several chromosomes. Following treatments with high salt (NaCl, 250 mM), PEG, low temperature (4 °C), and abscisic acid, the AdnsLTPs’ expression levels were altered. Three AdnsLTPs were linked to nematode infection resistance. The DOF and WRI1 transcription factors were suggested as potential controllers of the AdnsLTP response to nematode infection.
Fang et al. [30] found 330 TansLTPs (Triticum aestivum nsLTPs) genes in wheat (T. aestivum) (Table 1). Such a quantitative result can be considered an update of the 461 nsLTP loci found by Kouidri et al. [26] for the same species. To date, T. aestivum is the plant with the highest number of nsLTPs. The TansLTPs clustered into five groups (‘1’, ‘2’, ‘C’, ‘D’, and ‘G’) by phenetic analysis (Table 1). Gene structure and MEME pattern analyses showed that different groups of nsLTPs had similar structural compositions. Chromosome anchoring revealed that all five groups were distributed on 21 chromosomes. Furthermore, 31 gene clusters were identified as tandem duplications, and 208 gene pairs were identified as segmental duplications. Data mining of RNA-seq libraries, covering multiple stress conditions, showed that the transcript levels of some of the nsLTP genes could be strongly up-regulated by drought and high salt (NaCl, 250 mM) stresses.
In another context, Liang et al. [35] scrutinized the Brassica napus pangenome for BnnsLTPs (B. napus nsLTPs). These authors identified 246 BnnsLTP genes, divided into five groups (‘1’, ‘2’, ‘C’, ‘D’, and ‘G’; Table 1). Different BnnLTP genes were identified among the eight studied B. napus varieties (ZS11, Gangan, Zheyou7, Shengli, Tapidor, Quinta, Westar, and No2127). BnnsLTPs showed different duplication patterns in different varieties. Cis-regulatory elements that respond to biotic and abiotic stresses were anchored at all BnnsLTP genes. Finally, RNA-Seq analysis showed that the BnnsLTP genes were involved in responses to the fungus Sclerotinia sclerotiorum infection.
Vangelisti et al. [36], studying sunflower (Helianthus annuus) HansLTPs (Helianthus annuus nsLTPs), observed the existence of four (‘1’, ‘2’, ‘3’, and ‘4’) groups (Table 1). The authors did not explicitly classify the observed groups according to the available classification systems. The HansLTPs (101 in total) were further examined by looking into potential gene duplication sources, which revealed a high prevalence of tandem- in addition to whole-genome duplication (WGD) events. This finding is consistent with polyploidization events that occurred during the evolution of the sunflower genome. Three (‘1’, ‘3’, and ‘4’) of the four HansLTP groups responded uniquely to environmental cues, including auxin, abscisic acid, and the saline environment. Interestingly, sunflower seeds were the only source of expression for HansLTP group ‘2’ genes.
In line with the reports mentioned above and other works in Table 1, it is observed that the nsLTP genes act multifunctionally and show genetic variability even within accessions of the same species. nsLTPs are present in a wide range of plants, showing gene expression in different tissues, developmental stages, and stressful conditions.

2.2. Filling the Gap: Discovering and Classifying nsLTPs in New Plant Genomes

To provide genomic information for nsLTPs in plants not yet studied in the previous topic, and to update nsLTP data for some species with improved genome versions that have been made available, the following plant genomes were scrutinized (Table 2): (1) Marchantia polymorpha; (2) Ceratopteris richardii; (3) Selaginella moellendorffii; (4) Thuja plicata; (5) Gossypium hirsutum; (6) Lactuca sativa; (7) Manihot esculenta; (8) Mimulus guttatus; (9) Populus trichocarpa; (10) Sinapis alba; (11) Solanum tuberosum; and (12) Spinacea oleracea. The mentioned species were chosen to diversify the number of analyzed clades.
Three distinct and complementary strategies were used to retrieve nsLTPs in the genomes selected (information about the applied methodologies in Supplementary Material S1). The nsLTP exhaustive mining applied to 12 evaluated genomes (Table 2) returned 258 candidate sequences by BLASTp search, 344 by the cysteine pattern-based strategy (RegEx mining), and 1191 by the machine learning approach (HMMER tool).
The machine learning approach (HMMER tool) recovered a more comprehensive number of sequences considering nsLTP domains. Strategies based on machine learning are emerging as the future of DNA/RNA/protein sequence identification and bioinformatics, in general. The search with the RegEx approach was more restrictive, on the other hand, and not as accurate as the local alignment (BLASTp) method when observing the presence of the conserved nsLTP domain. However, it is worth mentioning that the BLASTp strategy did not recover sequences from the hypothetical proteomes of C. richardii, S. moellendorffii, and M. polymorpha (Table 2). When working with sequences as diverse as nsLTPs, it is advisable to combine several mining methods to increase the chances of finding the maximum number of sequences combined with subsequent data curation. In the present work, however, almost all of the sequences retrieved by the BLASTp and RegEx approaches were also retrieved by the machine-learning-based strategy .
From the analyzed species pool (Table 2), S. moellendorffii (21) and G. hirsutum (218) presented, respectively, the lowest and highest number of nsLTPs. Our search confirmed the tendency of these peptides to be encoded by large gene families: nine (75%) of the 12 analyzed species had more than 50 nsLTP loci in their respective genomes (Table 2). There was no correlation (r = 0.10; methodology in Supplementary Material S1) between genome size and the number of nsLTPs in the analyzed species pool (Table 2). However, angiosperms and the analyzed gymnosperm have a higher amount of nsLTPs than pteridophytes and the analyzed bryophyte (Table 2). This fact may be associated with the evolution of these basal groups, as shown by Edstam et al. [1]. Pteridophytes and bryophytes are phylogenetically closer to green algae than the other clades analyzed. This scenario is possibly responsible for the nsLTPs’ reduced number in the mentioned clades since no nsLTPs were identified in green algae [1].
The present work also updated the data for the nsLTP content in cotton (G. hirsutum), potato (S. tuberosum), common liverwort (M. polymorpha), and spikemoss (S moellendorffii) from the respective updated genome versions (Ghirsutum_527_v2.1, Stuberosum_686_v6.1, Mpolymorpha_320_v3.1, and Smoellendorffii_91_v1.0 | Phytozome database). While Li et al. [6] and Li et al. [29] identified, respectively, 91 and 83 nsLTP loci for the G. hirsutum and S. tuberosum (Table 1), 218 and 105 nsLTP loci were identified in our data mining approach (Table 2). For M. polymorpha and S. moellendorffii, 21 and 36 nsLTPs were found, respectively (Table 2), in the present study (compared to 13 and 23 nsLTPs reported by Fonseca-García et al. [34]; (Table 1)). Genome assemblies are never perfect since they are models for the actual genome. It is hard to completely rule out all potential technological or algorithmic flaws, and no single assembly can accurately capture all the variety within populations of a species. Thus, published genomes that have an active research community are continuously improved. Such modified and updated versions are potential sources for changing minor paradigms.
Considering the 1191 nsLTPs’ categorization (nsLTP domains are available in Supplementary Material S2), we observed nine large groups, of which five could be identified (‘nsLTL1’, ‘nsLTL2’, ‘nsLTLG’, ‘nsLTLD’, and ‘nsLTLC’; Figure 1), in addition to four distinct groups denominated ‘Unknown 1-4’ (Table 2). The complete tree, with bootstrap values and identified sequences with assigned categories, is available in Supplementary Figure S2.
Group separation, obtained from the implemented neighbor-joining (NJ) approach (methodology in Supplementary Material S1), reprised the nine groups in the Edstam et al. [1] classification. The strategy of performing NJ analysis from 8CM domain sequences (as performed by Edstam et al. [1] and Xue et al. [38]) promoted better group discrimination. The NJ tree derived from the complete nsLTP sequences (8CM + upstream and downstream regions) did not present such a level of resolution (only six groups were formed.
As will be seen in the ‘How structural nsLTP proteomics correlates with current nsLTP classification systems?’ section, nsLTP sequences show high variability in the amino acid sequence. Outside the 8CM region, the mentioned variability is accentuated, which causes greater noise in the distance analysis, resulting in an efficiency reduction in the formation of ‘true’ groups. Although reduced, compared to other nsLTP regions, the variability of the 8CM region was also a significant factor in the analysis using the NJ method. This is evidenced by the low bootstrap values of the first branches formed (Supplementary Figure S2) and in nsLTP NJ analysis for a manifold of species (see works by Fang et al. [30] and Vangelisti et al. [36], among others, in Table 1). Bootstrap values reflect the proportion of trees/replicates in which a recovered grouping is presented (in other words, a measure of support for that group). Despite reduced bootstrap values, the obtained tree topology was in accordance with the composition of the characterized seed sequences (methodology in Supplementary Material S1) used to perform the nsLTP classification.
Regarding the nsLTPs’ composition in the scrutinized species, L. sativa, despite not having the highest amount of nsLTPs in its genome (Table 2 and Table 3), was the species which presented the greatest variety of these peptides (Table 3), with at least one member of each of the nine groups found. M. polymorpha, in turn, had the lowest variety of nsLTPs (presenting members only for the ‘LTP1’, ‘LTP2’, ‘LTPD’, ‘LTPG’, ‘Unknown 1’, and ‘Unknown 4’ groups; Table 3), a fact that is possibly associated with its nsLTPs’ small genomic amount. Our results for this species, however, indicated that it presented nsLTP groups (e.g., ‘LTP1’, ‘LTP2’, ‘Unknown 1’, and ‘Unknown 4’) not yet identified in previous studies (i.e., in Edstam et al. [1] and Fonseca-García et al. [34]).
In another context, considering the different groups, ‘nsLTPG’ was present in the 12 analyzed plant species, while the ‘Unknown 4’ group was specific to four species (C. richardii, L. sativa, M. polymorpha, and S. moellendorffii; Table 3; Supplementary Figure S2). Subgroups were formed for ‘nsLTPG’ and ‘nsLTPD’ (Figure 1; Supplementary Figure S2). In the case of ‘Unknown 1-4’ (Figure 1; Supplementary Figure S2), no classification could be assigned based on the set of used seed sequences. The topological organization and sequence composition of the ‘Unknown’ groups needs further investigation.

3. ‘nsLTP Evolution’ Section

3.1. The Landscape of the nsLTP Expansion in Plant Genomes

Gene duplication is a crucial evolutionary mechanism providing new genetic material to a genome. New proteins often arise from duplicated copies generating unprecedented capabilities by diversifying protein functionality (paralogy). We analyzed and categorized such mechanisms (methodology in Supplementary Material S1) for the 12 species presented in Table 1. In addition, soybean (G. max) was also analyzed in this context. Although already studied in terms of nsLTP mining and classification [34], this species does not present information about nsLTPs’ expansion mechanisms.
Singleton nsLTPs were a minority in the analyzed genomes (Figure 2), indicating a high duplication rate of loci associated with this gene family. In the present study, two main mechanisms with different performance levels were identified as responsible for nsLTP expansion (Figure 2): dispersed and tandem duplications. Despite different duplication mechanisms, it is clear that similar evolutionary forces acted in the nsLTPs’ expansion in these different genomes.
Dispersed duplicated genes are prevalent in different plant genomes [19]. They generate paralogs that are neither near each other on chromosomes nor show conserved synteny [39]. This is in line with the random nature of the distribution of the analyzed loci in their respective anchor species.
Regarding the second most active mechanism, gene tandem duplication has been commonly found to be important for plant adaptation to rapidly changing environments. In contrast, genes expanded by non-tandem mechanisms tend to have intracellular regulatory roles [40]. Thus, the set of duplication mechanisms observed (mostly tandem and dispersed) may be associated with the heterogeneous physiological role that nsLTPs play in plants. The mentioned peptide group plays roles ranging from plant defense (for a review, see Liu et al. [8]) to liquid secretion, cuticular wax accumulation, pollen, and seed development, among other functions [20].
Considering the genesis of the distribution pattern, the dispersed one would be more expected in polyploids than in diploids. However, this assumption was not confirmed when comparing the two pteridophytes since there was a significant prevalence of dispersed distribution in S. moellendorffii (diploid with 2n = 16), while in the polyploid C. richardii (2n = 78), the two categories dispersed (44) and tandem (41) had similar prevalence (Figure 2).
In the diploid gymnosperm T. plicata (2n = 22), the species with the largest genome among those analyzed (ca. 12.5 Gb), the tandem dispersal mechanism clearly prevailed (Figure 2), indicating a possible importance of segmental duplications in the genesis of the high number of nsLTPs in this species.
Concerning the angiosperms (Figure 2), the species with the highest number of nsLTPs (218) was cotton (G. hirsutum, 2n = 52) with a prevalence of the dispersed mechanism, indicating the importance of WGD (i.e., polyploidy) in the expansion of nsLTPs. A similar result was observed in soybean (G. max), also polyploid (2n = 40). In contrast, this trend was not observed in potato (S. tuberosum), also considered polyploid (2n = 48), where a prevalence of the tandem (50 genes) mechanism was observed, followed by the dispersed mechanism (38). Considering the above, the mechanism of distribution and expansion of nsLTPs (and other gene families) needs to be carefully analyzed for each taxon and lineage without generalizations. This is especially true considering that some polyploids undergo diploidization processes with the retention of duplicate copies of some regions and the elimination of other ones (e.g., Renny-Byfield et al. [41]).
Future actions, such as the transcriptional expression study of tandem duplicated loci, will add more information about the evolutionary impact of identified duplication events. According to Liu et al. [8], the duplicated nsLTP genes might keep certain crucial activities in the course of subsequent evolution. The identical expression patterns could be explained by the duplicated genes’ remarkably similar protein architecture. In turn, duplicated genes with a considerable variation in expression may result from duplication events with considerable alteration of gene regulation.

3.2. nsLTP Lowest Common Ancestor (LCA) Analysis

Understanding the evolution, diversity, and abundance of nsLTPs within plants is not a trivial task. On this point, evaluating available genomic data using bioinformatics approaches may help the nsLTP origin and diversification studies. Here, we use LCA analysis (methodology in Supplementary Material S1) for the first time with nsLTPs, generating an nsLTP phylogenetic tree according to a taxonomic rank obtained from the NCBI Taxonomy database (Figure 3) using the Edstam et al. [1] classification. The LCA approach corresponds to an evolutionary method based on orthologous groups, part of Ukkonen’s Algorithm to construct suffix trees, aiming to determine the common ancestor between two or more distinct organisms based on a taxonomic tree [42,43]. The branching pattern in a phylogenetic tree (Figure 3) reflects how species or other taxonomic groups evolved from a series of common ancestors.
The nsLTP sequences used to generate the phylogenetic tree were grouped in different taxonomic categories, including Magnoliophyta (914 matches), Lycopodiopsida (five matches), and Bryopsida (one match) for the Viridiplantae clade (Figure 3). The LCA analysis confirmed the presence of nsLTP classes in plants, with some of them (e.g., LTP1 and LTP2) widely distributed among higher plants (Figure 3), corroborating the works by Edstam et al. [1], Salminen et al. [20], and Missaoui et al. [44]. Additionally, no green algae were observed in our LCA tree, strengthening the hypothesis of nsLTP origin in land plants about 470 million years ago [1].
The scientific literature mentions that higher plants exhibit exclusive nsLTPs (with different cysteine patterns) when compared to non-vascular plants (such as the moss Physcomitrella patens). These nsLTPs show low sequence similarity to the nsLTP genes of non-vascular plants [1]. From this presumption, one question arises: why do higher plants present exclusive classes not present in primitive ones? Most of these nsLTP classes probably result from evolution during terrestrial environment occupation, likely linked to their role in stress adaptation and defense, as proposed by Edstam et al. [1] and Liu et al. [8].
The earliest nsLTP groups, found in all land plants, seem to be members of nsLTPd and nsLTPg, presenting probably a common origin. Furthermore, the LCA tree revealed that nsLTPs from some species tend to cluster together—based on their structural features—with members of other distinct taxa (for more information, see the ‘Structural proteomics’ section), especially when considering early diverging angiosperms [26,34,44]. This could indicate that duplication events in the nsLTP family occurred independently in the different plant groups, even at the intraspecific level, leading to multi-functionalization of nsLTPs. As reviewed by Amador et al. [7], this great diversity and apparent redundancy of nsLTPs in angiosperm genomes, and the neofunctionalization process, contribute to paralogs’ emergence and their multiplicity of functions.
The molecular evolution of plant nsLTPs evidenced by LCA (Figure 3) showed that nsLTP diversification was mainly based on gene duplications, a known driving evolutionary force and an important factor in increasing genome complexity. Duplication events are probably the main mechanism of nsLTP evolution, including tandem and dispersed duplications (as discussed in the ‘The landscape of the nsLTP expansion in plants’ topic). In addition to the results presented here, the works by Liang et al. [35] and Yang et al. [45], among others, also corroborate this statement.
The presence/occurrence of horizontal gene transfer to/from bacteria should also be addressed (Figure 3). Supplementary Figure S1 and Figure 3 point to nsLTP diversification based on multiple duplication events in Viridiplantae, possibly associated with terrestrial environment occupation by land plants, and it was directly associated with leaf cuticle formation to avoid dehydration. Possible duplication events are indicated by black dots (Figure 3 and Supplementary Figure S1), pointing out shared nsLTP subgroups among related taxa. These duplication events are thought to have affected the nsLTP gene family, allowing both conservation and divergence of gene function [27,44]. Our LCA analysis, using plant nsLTPs as seeds in alignment against RefSeq-NCBI, also identified a few proteins from bacteria (Gammaproteobacteria (28 matches) and Bacilli (one match)) as well as Metazoan (Mammalia (one match); Figure 3) that contained the nsLTP motif and domain. The mammalian protein sequence identified (accession number XP_036992713, Supplementary Figure S1) regards the Jamaican fruit-eating bat (Artibeus jamaicensis), showing high similarity (87.78%) to a plant nsLTP1 amino acid sequence. Since no other nsLTP genes were identified for Metazoans, we believe this protein sequence may be a plant contamination of the genomic Metazoan sequencing data. We also obtained hits for two different species of bacteria [Acinetobacter baumannii (WP_143045904; Supplementary Figure S1) and Paenibacillus sp. (WP_083442884; Supplementary Figure S1)], presenting similarity to true non-plant nsLTP proteins and conserved domains related to nsLTPs from plants. Here, we consider this discovery as a possible case of HGT (horizontal gene transfer) since HGT is a well-known and pervasive evolutionary mechanism in prokaryotes for different gene families that could improve the adaptive ability of prokaryotes in changing environments [29,46]. Anyway, contamination cannot be discarded in all situations. Therefore, the presence of such sequences homologous to plant nsLTPs demands further investigation and confirmation.
Our LCA results for nsLTPs do not indicate that the complete nsLTP domain is present in proteins from species outside the plant kingdom. It is clear that to elucidate all questions regarding this kind of ‘evolutionary history noise’ of the nsLTP gene family depends on a broader range of sequenced species deposited on databases that will contribute to future studies to elucidate nsLTP functions and relevant mechanisms associated with their evolutionary origin and diversification, even considering HGT events.

4. ‘Structural Proteomics’ Section

How Structural nsLTP Proteomics Correlate with Current nsLTP Classification Systems?

Although some nsLTPs’ three-dimensional structures have been determined in different plant species, most structural studies are limited to nsLTPs of types ‘1’ and ‘2’ [44,47,48]. In this analytical context, the other nsLTP types are still poorly addressed. To circumvent this problem, this work generated theoretical protein models for nsLTP sequences (methodology in Supplementary Material S1) employed in two landmark studies [1,17] and widely used as current classification systems for this protein group. Forty-nine theoretical models for nsLTPs provided (Supplementary Material S3) and classified (i.e., types ‘1’ to ‘9’, in addition to type ‘Y’) by Boutrot et al. [17] were obtained. All these referential models are shown in Figure 4A,B. Additionally, another 41 models were generated for nsLTPs provided (Supplementary Material S4) and classified (i.e., ‘LTP1’, ‘LTP2’, ‘LTPc’, ‘LTPd’, ‘LTPg’, ‘LTPe’, ‘LTPf’, ‘LTPh’, ‘LTPj’, and ‘LTPk’) by Edstam et al. [1]. All these referential models are available in Supplementary Figure S3A,B.
Besides studying the nsLTP structures, we aimed to analyze protein modeling data to test whether resemblances/differences in the three-dimensional conformation and physicochemical parameters (physicochemical variables in Supplementary Material S1) result in clusters similar to those obtained with the parameters adopted by Boutrot et al. [17] or Edstam et al. [1] in their classification systems.
In general, for both analyzed works, it was observed that the three-dimensional structures of the different nsLTP types/groups were similar (Figure 4A; Supplementary Figure S3A). Exceptions, however, were found for the theoretical models of type ‘3’ nsLTPs (Boutrot et al. [17] classification; Figure 4A) and the nsLTPc group (Edstam et al. [1] classification), which presented a different structural configuration, since that they have only two disulfide bonds (‘Cys1-Cys5’, ‘Cys2-Cys3’). These structures are mainly composed of α-helices, common in nsLTPs, which can vary from three to four and are connected by small loops and an unstructured C-terminal tail (Figure 4A; Supplementary Figure S3A).
The position of α-helices creates an internal hydrophobic cavity (Figure 4B; Supplementary Figure S3B), suitable for binding lipids and other hydrophobic molecules. In general, such a cavity is a conserved feature in nsLTPs [7]. It is noteworthy that the association between hydrophobic molecules and hydrophobic pockets occurs with different affinities, depending on the molecule size and structure. However, this association does not cause any change in the protein’s three-dimensional conformation [49].
Despite the structural conservation observed in the generated models (Figure 4A; Supplementary Figure S3A), the nsLTPs analyzed by Boutrot et al. [17] and Edstam et al. [1] presented a high diversity of amino acid residues among the conserved cysteine sites (Supplementary Tables S1 and S2). Such heterogeneity is commonly observed in nsLTPs (and other cysteine-rich peptides). It may be related to the fact that they are involved in a multitude of plant molecular mechanisms, including defense against pathogens [50].
Given the diversity in the nsLTPs’ amino acids’ linear composition, there is a need to understand the peptide ‘sequence–structure–classification’ relationship. Such information will enable the understanding of how the obtained three-dimensional structures and their physicochemical characteristics relate to the classification systems [1,17] scrutinized here.
Considering the above, when analyzing the RMSD distance, i.e., the measure of the average distance between the atoms (usually the backbone atoms) of superimposed protein models (RMSD values in Figure 4A, Supplementary Figure S3A, and Supplementary Tables S3 and S4), it was possible to observe that most representative present similar values. The exceptions, however, were type ‘3’ nsLTPs (Boutrot et al. [17] classification) and the ‘nsLTPc’ group (Edstam et al. [1] classification). As mentioned earlier, these nsLTP classes presented a different three-dimensional structure than their peers. The almost generalized similarity was also found in the structural PCA performed with the models. The results presented in Figure 5A,B are inconsistent with the Boutrot et al. [17] and Edstam et al. [1] classification systems since they showed a generalized similar three-dimensional conformation. The explanation for the clustering pattern obtained in Figure 5A,B is related to the conservation of the structure to the detriment of the amino acids’ diversity among the cysteine residues anchored in the different nsLTP types/groups. The three-dimensional structural conservation makes it challenging to differentiate the nsLTPs into distinct and isolated groups. However, it is worth noting that there is some incipient group formation, such as type ‘3’, type ‘5’ (except for 5.1, which has five amino acid residues less in the C-terminal tail compared to the others of type ‘5’), and type ‘6’ (Figure 5A). The mentioned types slightly stand out from the other structures in the principal components (PCs) area (Figure 5A). In the structural PCA considering the sequences classified by Edstam et al. [1], however, such a small separation does not occur (Figure 5B).
In another context, when analyzing the physicochemical properties with PCA (Figure 6A,B), the absence of individualized groups was also observed. As mentioned before, the nsLTP sequences are very variable in terms of linear amino acid composition. Such diversity also results in variant physicochemical properties. In the sequence similarity index (Supplementary Table S1), groups ‘4’, ‘5’, ‘8’, ‘9’, and ‘Y’ (Boutrot et al. [17] classification) presented higher indices of intragroup similarity. Considering the classification by Edstam et al. [1], there was no group formation based on the similarity index (Supplementary Table S3).
Despite the almost general conservation of tertiary structure (Figure 4A; Supplementary Figure S3A), the surface charge qualitatively varied (positive or negative, but with a preponderance of positive charges in the studied sets) among the nsLTP types/groups (Figure 7 and Supplementary Figure S4). In addition, the number of hydrophobic, acidic, basic, and neutral amino acid residues also was diverse. Such observations may be associated with the different nsLTP molecular roles.
Considering both PCA analyses (based on (i) physicochemical features of aa residues or (ii) the peptide tertiary structure), no correlation could be made with the classification systems proposed by Boutrot et al. [17] and Edstam et al. [1]. Even intragroup, there was significant variation in the physicochemical characteristics of the amino acid residues of nsLTPs. In turn, considerable conservation of the tertiary structure was observed (except for type 3, classification by Boutrot et al. [17], or type C, classification by Edstam et al. [1]). Even being not yet ideal, these classifications are less artificial than the initially proposed grouping [11,12] into ‘nsLTP1’ and ‘nsLTP2’ subfamilies—based mainly on molecular mass and sequence identity.
In summary, current classifications target mainly the linear structures of nsLTPs, without major implications considering their specific functions. A future association of their roles, in terms of functional omics, with the physical–chemical and tertiary characteristics uncovered here may help elaborate a less artificial classification system reflecting the function of nsLTPs.

5. nsLTP Transcriptional Expression: Soybean as a Case Study

Expression in Soybean Plants: From Baseline to Differential Expression

Baseline gene expression analysis in 14 tissues/developmental stages of soybean (methodology in Supplementary Material S1) revealed that GmnsTLPs (Glycine max nsLTPs) were expressed on all analyzed substrates (Figure 8). Eighty-three (~70%; Supplementary Table S5) of the 120 GmnsLTP loci (Supplementary Table S6) were expressed in at least one of the analyzed substrates. About 15% (13 GmnsLTP loci) were expressed in only one or two different tissues/developmental stages, while approximately 44% (37 GmnsLTP loci) were expressed in at least 10 different tissues or developmental stages. Eighteen percent (15 GmnsLTP loci) were expressed on all analyzed treatments. GmnsLTPs were expressed mainly in flower (59 GmnsLTP loci), followed by ‘one centimeter pod’ (57) and root (57) (Figure 8). The lowest representation of expressed GmnsLTP loci occurred in ‘seed 42DAF’ (Figure 8).
Regarding the group with the highest number of expressed loci, GmnsLTPg stood out in 13 of the 14 tissues/stages studied (Figure 8), followed by the GmLTP1 group in seven, and by the GmnsLTPd group, in six tissues/stages. Interestingly, none of the analyzed tissues/stages presented expression of GmnsLTPc or GmnsLTPe representatives. The explanation for their absence remains an open question. In another context, there was no direct relationship between tissues/developmental stages and expression level of the different tested loci; it was observed, however, that there were loci with high levels of expression (RPKM counts > 100) in several tissues (Supplementary Table S5). These may become future targets for gene manipulation, with consequent functional characterization of GmnsLTPs.
Our analysis showed that nsLTP1, nsLTPg, nsLTPd, and nsLTP2 groups are expressed in all studied situations. At the loci level, however, promiscuous expression of certain GmnsLTPs (Supplementary Table S5) was found, while others were explicitly expressed only in a given tissue/stage of development. Regarding the first-mentioned expression group, such ubiquity in soybean tissues may be associated with the defense role played by these proteins. Some nsLTPs have been classified as belonging to the pathogen-related (PR) protein group, specifically, the PR14 family [51]. A broad line of evidence indicates an nsLTP role in plant defense (for review, see Liu et al. [8]). Notably, the highest amount of expressed GmnsLTP loci was found in floral, root, and young fruit (‘one centimeter pod’) tissues. Flowers have been implicated as hubs of disease transmission [52], while root tissue is a site commonly inhabited by microorganisms (whether pathogens [53] or symbionts [54]). In turn, soybean pods are also substrates affected by important diseases (as in soybean pod blight [55]).
Regarding the GmnsLTP loci with specific expression, in addition to being able to perform defense functions, such proteins may be associated with less generalist functions, such as cuticular wax accumulation, liquid secretion, pollen and seed development, and seed germination [20]. As mentioned, nsLTPs are considered critical proteins for the plant’s survival and land colonization.
Regarding the differential expression approach, four assays were analyzed (methodology in Supplementary Material S1) covering 20 different treatments (Figure 9; Supplementary Table S7). The expression of 80 distinct GmnsLTP loci were detected (Figure 9; Supplementary Table S7).
nsLTPs are often associated with defense against pathogens (e.g., Xue et al. [38]). Interestingly, there was no detection of differential expression (−1 > Log2FC > 1 and FDR < 0.05) after inoculation with Fusarium oxysporum (pathogenic and non-pathogenic isolates) in soybean root tissue (up to 96 h post-inoculation). No differential expression was observed in response to inoculation with MAMPs (flg22 and chitin) in two parental soybean lines and two progeny lines (from 1 to 8; Figure 9). Detailed analysis of these assays showed that some GmnsLTPs exhibited transcriptional modulation (Log2FC) with associated FDR < 0.05. However, the mentioned Log2FCs did not show values above the cut-off (−1 > Log2FC > 1) limits. Thus, a possible reason for the absence of up-regulation for assays involving biotic conditions in soybean may be associated with a methodological bias of the differential expression.
The mining of transcripts encoding nsLTP proteins in assays under abiotic conditions, however, revealed the up-regulation of several GmnsTLPs, namely: 3, for the treatment ‘dehydration stress’ at ‘12 h’ vs. ‘control’ (number 10; Figure 9); 15 for the treatment ‘salt stress’ at ‘12 h’ vs. ‘control’ (number 13; Figure 9); 8 for the treatment ‘salt stress’ at ‘6 h’ vs. ‘control’ (number 14; Figure 9); and 2, for the treatment ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 0’ (number 15; Figure 9). The mentioned amount was associated with 23 non-redundant loci (indicated by red spheres in Figure 9).
The studied data demonstrated an important association between GmnsLTPs and salt stress response, highlighting the GmnsTLPd and GmnsTLPg groups. These presented, respectively, eight and five up-regulated loci for the treatment of ‘salt stress’ at ‘12 h’ vs. ‘control’ (number 13; Figure 9), as well as four and three up-regulated loci for the treatment ‘salt stress’ at ‘6 h’ vs. ‘control’ (number 14; Figure 9). The nsTLPs’ participation in the response to saline conditions has been reported in the scientific literature. Akhiyarova et al. [56] assumed that these proteins could participate in salt-induced pea root suberization or transport of phloem lipid molecules. In turn, Xu et al. [57] observed that nsLTP enhances salt and drought tolerance in N. tabacum. Such information adds importance to this peptide group, reinforcing its multifunctional role in plants.

6. Conclusions, Perspectives, and Open Questions

There is growing interest in nsLTPs due to their multifunctional roles in plant molecular physiology, which has resulted in some comprehensive reviews regarding their structure and function. Due to this, the present work has proposed a different approach. In addition, to recover published information about the omics of the mentioned gene family, we focused on investigating topics such as evolution and structure using strategies that are still little applied for nsLTPs (as is the case with LCA, expansion mechanisms, and high-throughput protein modeling).
Our investigation corroborates previous works confirming the absence of the nsLTPs and their protein domain in species outside the plant kingdom. Sequences from other groups in databanks with similarity to nsLTPs were considered contamination artifacts (as in the case of the bat Artibeus jamaicensis). Concerning the sequences found in bacteria (A. baumannii and Paenibacillus sp.), they can also be contaminants, or may result from horizontal gene transfer, which needs experimental confirmation.
Tandem and dispersed duplication events figure as the main nsLTP evolution mechanisms in plants. Interestingly the number of nsLTP representatives may vary when comparing plants with similar complexity, but they do not ever correlate with the genome size or the ploidy level.
The gene expression analysis in soybean revealed a complex transcriptional regulation of these peptides, which can be transcribed in a specific or a promiscuous way (considering tissues or developmental stages). Additionally, we emphasized nsLTPs’ up-regulation in different tissues under abiotic stresses (such as high salt and drought).
An insight regarding the ‘sequence–structure–classification’ relationship, in turn, revealed a high conservation of the nsLTPs’ three-dimensional structure, in contrast with a significant variation in its physicochemical characteristics (even intragroup). It is clear that the initial classification system (dividing them into nsLTP1 and 2) is artificial. Physico-chemical variables should be associated with phenetics and functional features to generate a functional classification system for this peptide group.
Several challenges need to be overcome before we can completely comprehend the roles that nsLTPs play in plants, despite all the efforts made to uncover the diverse spectrum of functions of these molecules. We still cannot explain why nsLTPs exert diverse activities in plants despite sharing the considerable number of structural similarities uncovered here. Research into nsLTP function and evolution shall be accelerated by systems biology studies, with further downstream biochemical analysis and transgenic (or gene edition) research, which will also clarify the precise functions in planta of these peptides. The informational body available here enriches our comprehension of this heterogenic, plant-specific peptide family.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/antibiotics12050939/s1, Figure S1: Radial presentation of lowest common ancestor (LCA) dendrogram, representing relationships from 950 different RefSeq-NCBI nsLTP sequences. Black dots represent points of nsLTP gene duplication. Legend: n(number): LCA levels; * LCA levels not found in the obtained tree; black dots indicate possible duplication events; Figure S2: Neighbor-joining (NJ) tree (based in 8CM nsLTP domains), showing bootstrap values and sequences identification of 1191 non-specific lipid transfer proteins (nsLTPs) predicted in Marchantia polymorpha, Ceratopteris richardii, Selaginella moellendorffii, Thuja plicata, Glycine max, Gossypium hirsutum, Lactuca sativa, Manihot esculenta, Mimulus guttatus, Populus trichocarpa, Sinapis alba, Solanum tuberosum, and Spinacea oleracea genomes. All amino acid sequences were aligned using ClustalX2, and the NJ tree was constructed using the bootstrap neighbor-joining method (1000 replicates). The obtained result was visualized with the iTOL program; Figure S3: Three-dimensional reference structures for different types of nLTPs classified by Edstam et al. [1] and visualized in Pymol. A. Theoretical models’ alignment of different nsLTP types, solved by Alpha-Fold2, presenting their respective RMSD (root-mean-square deviation of atomic positions) values. For each presented structure, black arrows indicate the location of the N- and C-terminal sites and α-helices. B. Hydrophobic clusters, analyzed by ProteinTools, relative to the theoretical three-dimensional reference models of the different nsLTP groups classified by Edstam et al. [1]. The colored spheres represent different clusters formed from interactions between nearby amino acid residues; Figure S4: Adaptive Poisson–Boltzmann solver (APBS)-generated electrostatic surface potential plotted on the solvent accessible surface for nsLTP theoretical three-dimensional models (sequences from Edstam et al. [1]); Supplementary Material S1: Methodology used for nsLTP manuscripts data mining and carrying out the original analyses of the present work; Supplementary Material S2: 8CM domain sequences of the different used species (12 different species + Medicago truncatula, Arabidopsis thaliana, and Oryza sativa categorized nsLTP seed-sequences) to obtain the neighbor-joining tree; Supplementary Material S3: nsLTP protein sequences from Boutrot et al. [17]; Supplementary Material S4: nsLTP protein sequences from Edstam et al. [1]; Table S1: Sequence similarity of nsLTP sequences from Boutrot et al. [17]. Highlighted in red: similarity below 50% between the correlated sequences. Highlighted in green: similarity above 50% between the correlated sequences; Table S2: Sequence similarity of nsLTP sequences from Edstam et al. [1]. Highlighted in red: similarity below 50% between the correlated sequences. Highlighted in green: similarity above 50% between the correlated sequences; Table S3: Root-mean-square deviation (RMSD) between nsLTP sequences from Boutrot et al. [17]. In red: structures presenting greater distances compared to the others; Table S4: Root-mean-square deviation (RMSD) between nsLTP sequences from Edstam et al. [1]. In red: Structures presenting greater distances from the others; Table S5: nsLTP baseline expression analysis in different soybean tissues, presenting the respective expressed loci ID (for Glyma 1.1 and Glyma 2.0 assemblies), nsLTP name, nsLTP group (Edstam et al. [1] classification), and RPKM-normalized log2-transformed counts; Table S6: Information about soybean nsLTPs, presenting the loci ID and classification by Edstam et al. [1]; Table S7: nsLTP differential expression assays (see the different Excel tabs) in soybean under different biotic challenges, presenting the respective expressed loci ID (for Glyma 2.0 assembly), nsLTP name (Edstam et al. [1] classification), Log2-fold change, and adjusted p-values.

Author Contributions

C.A.d.S.-S. and J.R.C.F.-N. analyzed the data and wrote the manuscript; V.C.A., M.D.d.S., J.P.B.-N., L.M.B.V., E.B., M.d.S.R., A.L.T.M.d.M. and R.H.d.S. analyzed the data; and A.M.B.-I. conceived the project and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

FACEPE (Fundação de Amparo à Ciência e Tecnologia do Estado de Pernambuco)—grant number: BFP-0021-21-57015; CNPq (Conselho Nacional de Desenvolvimento Científico Tecnológico)—grant numbers: 406048/2022-3, 313581/2020-7, 442019/2019-0, and 433931/2018-3; CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior)—grant number: 88881.507105/2020-01.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the Supplementary Materials section.

Acknowledgments

The authors acknowledge the Fundação de Amparo à Pesquisa do Estado de Pernambuco (FACEPE), Conselho Nacional de Desenvolvimento Científico Tecnológico (CNPq), and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) for fellowships and financial support.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Edstam, M.M.; Viitanen, L.; Salminen, T.A.; Edqvist, J. Evolutionary History of the Non-Specific Lipid Transfer Proteins. Mol. Plant 2011, 4, 947–964. [Google Scholar] [CrossRef] [PubMed]
  2. Yeats, T.H.; Rose, J.K.C. The Biochemistry and Biology of Extracellular Plant Lipid-Transfer Proteins (LTPs). Protein Sci. 2008, 17, 191–198. [Google Scholar] [CrossRef] [PubMed]
  3. Hauser, M.; Roulias, A.; Ferreira, F.; Egger, M. Panallergens and Their Impact on the Allergic Patient. Allergy Asthma Clin. Immunol. 2010, 6, 1–14. [Google Scholar] [CrossRef]
  4. Edstam, M.M.; Edqvist, J. Involvement of GPI-Anchored Lipid Transfer Proteins in the Development of Seed Coats and Pollen in Arabidopsis thaliana. Physiol. Plant. 2014, 152, 32–42. [Google Scholar] [CrossRef]
  5. Wei, K.; Zhong, X. Non-Specific Lipid Transfer Proteins in Maize. BMC Plant Biol. 2014, 14, 281. [Google Scholar] [CrossRef]
  6. Li, F.; Fan, K.; Ma, F.; Yue, E.; Bibi, N.; Wang, M.; Shen, H.; Hasan, M.M.-U.; Wang, X. Genomic Identification and Comparative Expansion Analysis of the Non-Specific Lipid Transfer Protein Gene Family in Gossypium. Sci. Rep. 2016, 6, 38948. [Google Scholar] [CrossRef]
  7. Amador, V.C.; Santos-Silva, C.A.D.; Vilela, L.M.B.; Oliveira-Lima, M.; De Santana Rêgo, M.; Roldan-Filho, R.S.; Oliveira-Silva, R.L.D.; Lemos, A.B.; De Oliveira, W.D.; Ferreira-Neto, J.R.C.; et al. Lipid Transfer Proteins (LTPs)—Structure, Diversity and Roles beyond Antimicrobial Activity. Antibiotics 2021, 10, 1281. [Google Scholar] [CrossRef]
  8. Liu, F.; Zhang, X.; Lu, C.; Zeng, X.; Li, Y.; Fu, D.; Wu, G. Non-Specific Lipid Transfer Proteins in Plants: Presenting New Advances and an Integrated Functional Analysis. J. Exp. Bot. 2015, 66, 5663–5681. [Google Scholar] [CrossRef] [PubMed]
  9. Srivastava, S.; Dashora, K.; Ameta, K.L.; Singh, N.P.; El-Enshasy, H.A.; Pagano, M.C.; Hesham, A.E.; Sharma, G.D.; Sharma, M.; Bhargava, A. Cysteine-rich Antimicrobial Peptides from Plants: The Future of Antimicrobial Therapy. Phytother. Res. 2021, 35, 256–277. [Google Scholar] [CrossRef]
  10. Oliveira-Lima, M.; Benko-Iseppon, A.; Neto, J.; Rodriguez-Decuadro, S.; Kido, E.; Crovella, S.; Pandolfi, V. Snakin: Structure, Roles and Applications of a Plant Antimicrobial Peptide. Curr. Protein Pept. Sci. 2017, 18, 368–374. [Google Scholar] [CrossRef]
  11. Kalla, R.; Shimamoto, K.; Potter, R.; Nielsen, P.S.; Linnestad, C.; Olsen, O.-A. The Promoter of the Barley Aleurone-Specific Gene Encoding a Putative 7 KDa Lipid Transfer Protein Confers Aleurone Cell-Specific Expression in Transgenic Rice. Plant J. 1994, 6, 849–860. [Google Scholar] [CrossRef] [PubMed]
  12. Douliez, J.-P.; Michon, T.; Elmorjani, K.; Marion, D. Mini Review: Structure, Biological and Technological Functions of Lipid Transfer Proteins and Indolines, the Major Lipid Binding Proteins from Cereal Kernels. J. Cereal Sci. 2000, 32, 1–20. [Google Scholar] [CrossRef]
  13. Pasquato, N.; Berni, R.; Folli, C.; Folloni, S.; Cianci, M.; Pantano, S.; Helliwell, J.R.; Zanotti, G. Crystal Structure of Peach Pru p 3, the Prototypic Member of the Family of Plant Non-Specific Lipid Transfer Protein Pan-Allergens. J. Mol. Biol. 2006, 356, 684–694. [Google Scholar] [CrossRef] [PubMed]
  14. Hoh, F.; Pons, J.-L.; Gautier, M.-F.; De Lamotte, F.; Dumas, C. Structure of a Liganded Type 2 Non-Specific Lipid-Transfer Protein from Wheat and the Molecular Basis of Lipid Binding. Acta Crystallogr. D Biol. Crystallogr. 2005, 61, 397–406. [Google Scholar] [CrossRef] [PubMed]
  15. Gincel, E.; Simorre, J.-P.; Caille, A.; Marion, D.; Ptak, M.; Vovelle, F. Three-Dimensional Structure in Solution of a Wheat Lipid-Transfer Protein from Multidimensional 1H-NMR Data. A New Folding for Lipid Carriers. Eur. J. Biochem. 1994, 226, 413–422. [Google Scholar] [CrossRef]
  16. Lerche, M.H.; Poulsen, F.M. Solution Structure of Barley Lipid Transfer Protein Complexed with Palmitate. Two Different Binding Modes of Palmitate in the Homologous Maize and Barley Nonspecific Lipid Transfer Proteins. Protein Sci. 1998, 7, 2490–2498. [Google Scholar] [CrossRef]
  17. Boutrot, F.; Chantret, N.; Gautier, M.-F. Genome-Wide Analysis of the Rice and Arabidopsis Non-Specific Lipid Transfer Protein (NsLtp) Gene Families and Identification of Wheat NsLtp Genes by EST Data Mining. BMC Genom. 2008, 9, 86. [Google Scholar] [CrossRef]
  18. Jülke, S.; Ludwig-Müller, J. Response of Arabidopsis thaliana Roots with Altered Lipid Transfer Protein (LTP) Gene Expression to the Clubroot Disease and Salt Stress. Plants 2015, 5, 2. [Google Scholar] [CrossRef]
  19. Wang, H.W.; Hwang, S.-G.; Karuppanapandian, T.; Liu, A.; Kim, W.; Jang, C.S. Insight into the Molecular Evolution of Non-Specific Lipid Transfer Proteins via Comparative Analysis Between Rice and Sorghum. DNA Res. 2012, 19, 179–194. [Google Scholar] [CrossRef]
  20. Salminen, T.A.; Blomqvist, K.; Edqvist, J. Lipid Transfer Proteins: Classification, Nomenclature, Structure, and Function. Planta 2016, 244, 971–997. [Google Scholar] [CrossRef]
  21. Liu, W.; Huang, D.; Liu, K.; Hu, S.; Yu, J.; Gao, G.; Song, S. Discovery, Identification and Comparative Analysis of Non-Specific Lipid Transfer Protein (NsLtp) Family in Solanaceae. Genom. Proteom. Bioinform. 2010, 8, 229–237. [Google Scholar] [CrossRef]
  22. Tapia, G.; Morales-Quintana, L.; Parra, C.; Berbel, A.; Alcorta, M. Study of NsLTPs in Lotus Japonicus Genome Reveal a Specific Epidermal Cell Member (LjLTP10) Regulated by Drought Stress in Aerial Organs with a Putative Role in Cutin Formation. Plant Mol. Biol. 2013, 82, 485–501. [Google Scholar] [CrossRef]
  23. Li, J.; Gao, G.; Xu, K.; Chen, B.; Yan, G.; Li, F.; Qiao, J.; Zhang, T.; Wu, X. Genome-Wide Survey and Expression Analysis of the Putative Non-Specific Lipid Transfer Proteins in Brassica rapa L. PLoS ONE 2014, 9, e84556. [Google Scholar] [CrossRef]
  24. Ji, J.; Lv, H.; Yang, L.; Fang, Z.; Zhuang, M.; Zhang, Y.; Liu, Y.; Li, Z. Genome-Wide Identification and Characterization of Non-Specific Lipid Transfer Proteins in Cabbage. PeerJ 2018, 6, e5379. [Google Scholar] [CrossRef] [PubMed]
  25. Hairat, S.; Baranwal, V.K.; Khurana, P. Identification of Triticum Aestivum NsLTPs and Functional Validation of Two Members in Development and Stress Mitigation Roles. Plant Physiol. Biochem. 2018, 130, 418–430. [Google Scholar] [CrossRef] [PubMed]
  26. Kouidri, A.; Whitford, R.; Suchecki, R.; Kalashyan, E.; Baumann, U. Genome-Wide Identification and Analysis of Non-Specific Lipid Transfer Proteins in Hexaploid Wheat. Sci. Rep. 2018, 8, 17087. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, M.; Kim, Y.; Zong, J.; Lin, H.; Dievart, A.; Li, H.; Zhang, D.; Liang, W. Genome-Wide Analysis of the Barley Non-Specific Lipid Transfer Protein Gene Family. Crop J. 2019, 7, 65–76. [Google Scholar] [CrossRef]
  28. D’Agostino, N.; Buonanno, M.; Ayoub, J.; Barone, A.; Monti, S.M.; Rigano, M.M. Identification of Non-Specific Lipid Transfer Protein Gene Family Members in Solanum Lycopersicum and Insights into the Features of Sola l 3 Protein. Sci. Rep. 2019, 9, 1607. [Google Scholar] [CrossRef]
  29. Li, G.; Hou, M.; Liu, Y.; Pei, Y.; Ye, M.; Zhou, Y.; Huang, C.; Zhao, Y.; Ma, H. Genome-Wide Identification, Characterization and Expression Analysis of the Non-Specific Lipid Transfer Proteins in Potato. BMC Genom. 2019, 20, 375. [Google Scholar] [CrossRef]
  30. Fang, Z.; He, Y.; Liu, Y.; Jiang, W.; Song, J.; Wang, S.; Ma, D.; Yin, J. Bioinformatic Identification and Analyses of the Non-Specific Lipid Transfer Proteins in Wheat. J. Integr. Agric. 2020, 19, 1170–1185. [Google Scholar] [CrossRef]
  31. Song, X.; Li, E.; Song, H.; Du, G.; Li, S.; Zhu, H.; Chen, G.; Zhao, C.; Qiao, L.; Wang, J.; et al. Genome-Wide Identification and Characterization of Nonspecific Lipid Transfer Protein (NsLTP) Genes in Arachis Duranensis. Genomics 2020, 112, 4332–4341. [Google Scholar] [CrossRef] [PubMed]
  32. Song, S.; You, J.; Shi, L.; Sheng, C.; Zhou, W.; Dossou, S.S.K.; Dossa, K.; Wang, L.; Zhang, X. Genome-Wide Analysis of NsLTP Gene Family and Identification of SiLTPs Contributing to High Oil Accumulation in Sesame (Sesamum Indicum L.). Int. J. Mol. Sci. 2021, 22, 5291. [Google Scholar] [CrossRef]
  33. Duo, J.; Xiong, H.; Wu, X.; Li, Y.; Si, J.; Zhang, C.; Duan, R. Genome-Wide Identification and Expression Profile under Abiotic Stress of the Barley Non-Specific Lipid Transfer Protein Gene Family and Its Qingke Orthologues. BMC Genomics 2021, 22, 674. [Google Scholar] [CrossRef] [PubMed]
  34. Fonseca-García, C.; Solis-Miranda, J.; Pacheco, R.; Quinto, C. Non-Specific Lipid Transfer Proteins in Legumes and Their Participation During Root-Nodule Symbiosis. Front. Agron. 2021, 3, 660100. [Google Scholar] [CrossRef]
  35. Liang, Y.; Huang, Y.; Chen, K.; Kong, X.; Li, M. Characterization of Non-Specific Lipid Transfer Protein (NsLtp) Gene Families in the Brassica Napus Pangenome Reveals Abundance Variation. BMC Plant Biol. 2022, 22, 21. [Google Scholar] [CrossRef]
  36. Vangelisti, A.; Simoni, S.; Usai, G.; Mascagni, F.; Ventimiglia, M.; Natali, L.; Cavallini, A.; Giordani, T. In Silico Genome-Wide Characterisation of the Lipid Transfer Protein Multigenic Family in Sunflower (H. Annuus L.). Plants 2022, 11, 664. [Google Scholar] [CrossRef]
  37. De Oliveira Silva, L.; Da Silva Pereira, L.; Pereira, J.L.; Gomes, V.M.; Grativol, C. Divergence and Conservation of Defensins and Lipid Transfer Proteins (LTPs) from Sugarcane Wild Species and Modern Cultivar Genomes. Funct. Integr. Genom. 2022, 22, 235–250. [Google Scholar] [CrossRef]
  38. Xue, Y.; Zhang, C.; Shan, R.; Li, X.; Tseke Inkabanga, A.; Li, L.; Jiang, H.; Chai, Y. Genome-Wide Identification and Expression Analysis of NsLTP Gene Family in Rapeseed (Brassica Napus) Reveals Their Critical Roles in Biotic and Abiotic Stress Responses. Int. J. Mol. Sci. 2022, 23, 8372. [Google Scholar] [CrossRef] [PubMed]
  39. Ganko, E.W.; Meyers, B.C.; Vision, T.J. Divergence in Expression between Duplicated Genes in Arabidopsis. Mol. Biol. Evol. 2007, 24, 2298–2309. [Google Scholar] [CrossRef]
  40. Hanada, K.; Zou, C.; Lehti-Shiu, M.D.; Shinozaki, K.; Shiu, S.-H. Importance of Lineage-Specific Expansion of Plant Tandem Duplicates in the Adaptive Response to Environmental Stimuli. Plant Physiol. 2008, 148, 993–1003. [Google Scholar] [CrossRef]
  41. Renny-Byfield, S.; Kovarik, A.; Kelly, L.J.; Macas, J.; Novak, P.; Chase, M.W.; Nichols, R.A.; Pancholi, M.R.; Grandbastien, M.-A.; Leitch, A.R. Diploidization and Genome Size Change in Allopolyploids Is Associated with Differential Dynamics of Low- and High-Copy Sequences. Plant J. 2013, 74, 829–839. [Google Scholar] [CrossRef] [PubMed]
  42. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.-H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A Toolkit for Detection and Evolutionary Analysis of Gene Synteny and Collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [PubMed]
  43. Bezerra-Neto, J.P.; De Araújo, F.C.; Ferreira-Neto, J.R.C.; Da Silva, M.D.; Pandolfi, V.; Aburjaile, F.F.; Sakamoto, T.; De Oliveira Silva, R.L.; Kido, E.A.; Barbosa Amorim, L.L.; et al. Plant Aquaporins: Diversity, Evolution and Biotechnological Applications. Curr. Protein Pept. Sci. 2019, 20, 368–395. [Google Scholar] [CrossRef] [PubMed]
  44. Missaoui, K.; Gonzalez-Klein, Z.; Jemli, S.; Garrido-Arandia, M.; Diaz-Perales, A.; Tome-Amat, J.; Brini, F. Identification and Molecular Characterization of a Novel Non-Specific Lipid Transfer Protein (TdLTP2) from Durum Wheat. PLoS ONE 2022, 17, e0266971. [Google Scholar] [CrossRef] [PubMed]
  45. Yang, Y.; Li, P.; Liu, C.; Wang, P.; Cao, P.; Ye, X.; Li, Q. Systematic Analysis of the Non-Specific Lipid Transfer Protein Gene Family in Nicotiana Tabacum Reveal Its Potential Roles in Stress Responses. Plant Physiol. Biochem. 2022, 172, 33–47. [Google Scholar] [CrossRef] [PubMed]
  46. Soucy, S.M.; Huang, J.; Gogarten, J.P. Horizontal Gene Transfer: Building the Web of Life. Nat. Rev. Genet. 2015, 16, 472–482. [Google Scholar] [CrossRef] [PubMed]
  47. Abdullah, S.U.; Alexeev, Y.; Johnson, P.E.; Rigby, N.M.; Mackie, A.R.; Dhaliwal, B.; Mills, E.N.C. Ligand Binding to an Allergenic Lipid Transfer Protein Enhances Conformational Flexibility Resulting in an Increase in Susceptibility to Gastroduodenal Proteolysis. Sci. Rep. 2016, 6, 30279. [Google Scholar] [CrossRef]
  48. Megeressa, M.; Siraj, B.; Zarina, S.; Ahmed, A. Structural Characterization and in Vitro Lipid Binding Studies of Non-Specific Lipid Transfer Protein 1 (NsLTP1) from Fennel (Foeniculum Vulgare) Seeds. Sci. Rep. 2020, 10, 21243. [Google Scholar] [CrossRef]
  49. Prinz, W.A. The Lipid Trade. Nat. Rev. Mol. Cell Biol. 2014, 15, 79. [Google Scholar] [CrossRef]
  50. Santos-Silva, C.A.D.; Zupin, L.; Oliveira-Lima, M.; Vilela, L.M.B.; Bezerra-Neto, J.P.; Ferreira-Neto, J.R.; Ferreira, J.D.C.; Oliveira-Silva, R.L.D.; Pires, C.D.J.; Aburjaile, F.F.; et al. Plant Antimicrobial Peptides: State of the Art, In Silico Prediction and Perspectives in the Omics Era. Bioinforma. Biol. Insights 2020, 14, 117793222095273. [Google Scholar] [CrossRef]
  51. Sels, J.; Mathys, J.; De Coninck, B.M.A.; Cammue, B.P.A.; De Bolle, M.F.C. Plant Pathogenesis-Related (PR) Proteins: A Focus on PR Peptides. Plant Physiol. Biochem. 2008, 46, 941–950. [Google Scholar] [CrossRef] [PubMed]
  52. Figueroa, L.L.; Blinder, M.; Grincavitch, C.; Jelinek, A.; Mann, E.K.; Merva, L.A.; Metz, L.E.; Zhao, A.Y.; Irwin, R.E.; McArt, S.H.; et al. Bee Pathogen Transmission Dynamics: Deposition, Persistence and Acquisition on Flowers. Proc. R. Soc. B Biol. Sci. 2019, 286, 20190603. [Google Scholar] [CrossRef]
  53. Termorshuizen, A.J. Root Pathogens. In Interactions in Soil: Promoting Plant Growth; Springer: Dordrecht, The Netherlands, 2014; pp. 119–137. [Google Scholar]
  54. O’Brien, A.M.; Harrison, T.L. Host Match Improves Root Microbiome Growth. Nat. Microbiol. 2021, 6, 1103–1104. [Google Scholar] [CrossRef]
  55. Chen, C.H.; Wang, T.C.; Seo, M.J. First Report of Soybean Pod and Stem Blight Caused by Diaporthe Phaseolorum Var. Sojae in Taiwan. Plant Dis. 2009, 93, 202. [Google Scholar] [CrossRef] [PubMed]
  56. Akhiyarova, G.R.; Finkina, E.I.; Ovchinnikova, T.V.; Veselov, D.S.; Kudoyarova, G.R. Role of Pea LTPs and Abscisic Acid in Salt-Stressed Roots. Biomolecules 2019, 10, 15. [Google Scholar] [CrossRef] [PubMed]
  57. Xu, Y.; Zheng, X.; Song, Y.; Zhu, L.; Yu, Z.; Gan, L.; Zhou, S.; Liu, H.; Wen, F.; Zhu, C. NtLTP4, a Lipid Transfer Protein That Enhances Salt and Drought Stresses Tolerance in Nicotiana Tabacum. Sci. Rep. 2018, 8, 8873. [Google Scholar] [CrossRef] [PubMed]
Antibiotics 12 00939 g001 550
Figure 1. Neighbor-joining tree (based on 8CM nsLTP domains) of 1191 non-specific lipid transfer proteins (nsLTPs) predicted in Marchantia polymorpha, Ceratopteris richardii, Selaginella moellendorffii, Thuja plicata, Glycine max, Gossypium hirsutum, Lactuca sativa, Manihot esculenta, Mimulus guttatus, Populus trichocarpa, Sinapis alba, Solanum tuberosum, and Spinacea oleracea genomes. All amino acid sequences were aligned using ClustalX2. The obtained result was visualized with the iTOL program. Legend: the numbers (1–4) inside the circles on the edge of the tree indicate different nsLTP subgroups of a given nsLTP group; at the rectangles, the nsLTP group classification is available.
Figure 1. Neighbor-joining tree (based on 8CM nsLTP domains) of 1191 non-specific lipid transfer proteins (nsLTPs) predicted in Marchantia polymorpha, Ceratopteris richardii, Selaginella moellendorffii, Thuja plicata, Glycine max, Gossypium hirsutum, Lactuca sativa, Manihot esculenta, Mimulus guttatus, Populus trichocarpa, Sinapis alba, Solanum tuberosum, and Spinacea oleracea genomes. All amino acid sequences were aligned using ClustalX2. The obtained result was visualized with the iTOL program. Legend: the numbers (1–4) inside the circles on the edge of the tree indicate different nsLTP subgroups of a given nsLTP group; at the rectangles, the nsLTP group classification is available.
Antibiotics 12 00939 g001
Antibiotics 12 00939 g002 550
Figure 2. Categorization and quantification of the nsLTPs’ main genomic expansion mechanisms in Marchantia polymorpha (Mp); Ceratopteris richardii (Cr); Selaginella moellendorffii (Sm); Thuja plicata (Tp); Glycine max (Gm); Gossypium hirsutum (Gh); Lactuca sativa (Ls); Manihot esculenta (Me); Mimulus guttatus (Mg); Populus trichocarpa (Pt); Sinapis alba (Sa); Solanum tuberosum (St); and Spinacea oleracea (So) genomes. Legend: Data presented in percentage values.
Figure 2. Categorization and quantification of the nsLTPs’ main genomic expansion mechanisms in Marchantia polymorpha (Mp); Ceratopteris richardii (Cr); Selaginella moellendorffii (Sm); Thuja plicata (Tp); Glycine max (Gm); Gossypium hirsutum (Gh); Lactuca sativa (Ls); Manihot esculenta (Me); Mimulus guttatus (Mg); Populus trichocarpa (Pt); Sinapis alba (Sa); Solanum tuberosum (St); and Spinacea oleracea (So) genomes. Legend: Data presented in percentage values.
Antibiotics 12 00939 g002
Antibiotics 12 00939 g003 550
Figure 3. Lowest common ancestor (LCA) dendrogram, representing relationships from 950 different RefSeq-NCBI nsLTP sequences. Black dots represent points of nsLTP gene duplication. Legend: n (number): LCA levels; * LCA levels not found in the obtained tree; black dots indicate possible duplication events.
Figure 3. Lowest common ancestor (LCA) dendrogram, representing relationships from 950 different RefSeq-NCBI nsLTP sequences. Black dots represent points of nsLTP gene duplication. Legend: n (number): LCA levels; * LCA levels not found in the obtained tree; black dots indicate possible duplication events.
Antibiotics 12 00939 g003
Antibiotics 12 00939 g004 550
Figure 4. Three-dimensional reference structures for different types of nLTPs classified by Boutrot et al. [17] and visualized in Pymol. (A) Theoretical model alignment of different nsLTP types, solved by Alpha-Fold2, presenting their respective RMSD (root-mean-square deviation of atomic positions) values. For each given structure, black arrows indicate the location of the N- and C-terminal sites and α-helices. (B) Hydrophobic clusters, analyzed by ProteinTools, relative to the theoretical three-dimensional reference models of the different nsLTP types classified by Boutrot et al. [17]. The green and red spheres represent different clusters formed from interactions between nearby amino acid residues.
Figure 4. Three-dimensional reference structures for different types of nLTPs classified by Boutrot et al. [17] and visualized in Pymol. (A) Theoretical model alignment of different nsLTP types, solved by Alpha-Fold2, presenting their respective RMSD (root-mean-square deviation of atomic positions) values. For each given structure, black arrows indicate the location of the N- and C-terminal sites and α-helices. (B) Hydrophobic clusters, analyzed by ProteinTools, relative to the theoretical three-dimensional reference models of the different nsLTP types classified by Boutrot et al. [17]. The green and red spheres represent different clusters formed from interactions between nearby amino acid residues.
Antibiotics 12 00939 g004
Antibiotics 12 00939 g005 550
Figure 5. Principal components analysis (PCA) of the different nsLTP types/groups’ three-dimensional structures. (A) Philomorpho-spatial graph of the 49 structural theoretical models of nsLTPs classified by Boutrot et al. [17]. (B) Philomorpho-spatial graph of the 41 structural theoretical models of nsLTPs classified by Edstam et al. [1]. Structures are colored based on the group corresponding to the principal component (PC).
Figure 5. Principal components analysis (PCA) of the different nsLTP types/groups’ three-dimensional structures. (A) Philomorpho-spatial graph of the 49 structural theoretical models of nsLTPs classified by Boutrot et al. [17]. (B) Philomorpho-spatial graph of the 41 structural theoretical models of nsLTPs classified by Edstam et al. [1]. Structures are colored based on the group corresponding to the principal component (PC).
Antibiotics 12 00939 g005
Antibiotics 12 00939 g006 550
Figure 6. Principal components analysis (PCA) of different nsLTP types/groups. (A) Philomorpho-spatial chart of the amino acid physicochemical characteristics of 216 nsLTPs classified by Boutrot et al. [17]. (B) Philomorpho-spatial chart of the amino acid physicochemical characteristics of 263 nsLTPs classified by Edstam et al. [1]. On the left, a two-dimensional PCA representation (each dot represents an nsLTP sequence); on the right, a three-dimensional PCA representation (each dot represents an nsLTP sequence; axes are the principal components ‘1’, ‘2’, and ‘3’ of the dimensional space). Structures are colored based on the group corresponding to the principal component (PC).
Figure 6. Principal components analysis (PCA) of different nsLTP types/groups. (A) Philomorpho-spatial chart of the amino acid physicochemical characteristics of 216 nsLTPs classified by Boutrot et al. [17]. (B) Philomorpho-spatial chart of the amino acid physicochemical characteristics of 263 nsLTPs classified by Edstam et al. [1]. On the left, a two-dimensional PCA representation (each dot represents an nsLTP sequence); on the right, a three-dimensional PCA representation (each dot represents an nsLTP sequence; axes are the principal components ‘1’, ‘2’, and ‘3’ of the dimensional space). Structures are colored based on the group corresponding to the principal component (PC).
Antibiotics 12 00939 g006
Antibiotics 12 00939 g007 550
Figure 7. Adaptive Poisson–Boltzmann solver (APBS)-generated electrostatic surface potential plotted on the solvent accessible surface for nsLTP theoretical three-dimensional models (classification of Boutrot et al. [17]), showing the predominance of surfaces with cationic charges for most of the peptides.
Figure 7. Adaptive Poisson–Boltzmann solver (APBS)-generated electrostatic surface potential plotted on the solvent accessible surface for nsLTP theoretical three-dimensional models (classification of Boutrot et al. [17]), showing the predominance of surfaces with cationic charges for most of the peptides.
Antibiotics 12 00939 g007
Antibiotics 12 00939 g008 550
Figure 8. Schematic representation of soybean plants indicating the 14 treatments (tissues/stages of development) in which the GmnsLTPs (Glycine max nsLTPs; Edstam et al. [1] classification) expression was monitored. The box centered at the bottom of the figure indicates, through color code, the different searched GmnsLTP groups. Quantities within colored circles represent the number of expressed coding loci for the corresponding group/treatment. Legend: * DAF (days after flowering).
Figure 8. Schematic representation of soybean plants indicating the 14 treatments (tissues/stages of development) in which the GmnsLTPs (Glycine max nsLTPs; Edstam et al. [1] classification) expression was monitored. The box centered at the bottom of the figure indicates, through color code, the different searched GmnsLTP groups. Quantities within colored circles represent the number of expressed coding loci for the corresponding group/treatment. Legend: * DAF (days after flowering).
Antibiotics 12 00939 g008
Antibiotics 12 00939 g009 550
Figure 9. Heat map showing the GmLTP‡ (for Glyma 2.0 assembly) expression under different stressful conditions: 1. ‘Fusarium oxysporum FO36 non-pathogenic’ vs. ‘control’ at ‘72 h’; 2. ‘Fusarium oxysporum FO36 non-pathogenic’ vs. ‘control’ at ‘96 h’; 3. ‘Fusarium oxysporum FO40 pathogenic’ vs. ‘control’ at ‘72 h’; 4. ‘Fusarium oxysporum FO40 pathogenic’ vs. ‘control’ at ‘96 h’; 5. ‘MAMP solution; LD’ vs. ‘mock; LD’; 6. ‘MAMP solution; LDX’ vs. ‘mock; LDX’; 7. ‘MAMP solution; RIL-11268’ vs. ‘mock; RIL-11268’; 8. ‘MAMP solution; RIL-11272’ vs. ‘mock; RIL-11272’; 9. ‘dehydration stress’ at ‘1 h’ vs. ‘control’ at ‘0 h’; 10. ‘dehydration stress’ at ‘12 h’ vs. ‘control’ at ‘0 h’; 11. ‘dehydration stress’ at ‘6 h’ vs. ‘control’ at ‘0 h’; 12. ‘salt stress’ at ‘1 h’ vs. ‘control’ at ‘0 h’; 13. ‘salt stress’ at ‘12 h’ vs. ‘control’ at ‘0 h’; 14. ‘salt stress’ at ‘6 h’ vs. ‘control’ at ‘0 h’; 15. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 0’; 16. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 12’; 17. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 16’; 18. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 20’; 19. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 4’; 20. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 8’. Legend: red spheres indicate up-regulated GmLTP loci (*Log2FC > 1 and FDR < 0.05) in at least one of the analyzed treatments; ‡ in chromosomal order, with abbreviated nomenclature (without the Glyma epithet), and followed by their respective nsLTP categorization.
Figure 9. Heat map showing the GmLTP‡ (for Glyma 2.0 assembly) expression under different stressful conditions: 1. ‘Fusarium oxysporum FO36 non-pathogenic’ vs. ‘control’ at ‘72 h’; 2. ‘Fusarium oxysporum FO36 non-pathogenic’ vs. ‘control’ at ‘96 h’; 3. ‘Fusarium oxysporum FO40 pathogenic’ vs. ‘control’ at ‘72 h’; 4. ‘Fusarium oxysporum FO40 pathogenic’ vs. ‘control’ at ‘96 h’; 5. ‘MAMP solution; LD’ vs. ‘mock; LD’; 6. ‘MAMP solution; LDX’ vs. ‘mock; LDX’; 7. ‘MAMP solution; RIL-11268’ vs. ‘mock; RIL-11268’; 8. ‘MAMP solution; RIL-11272’ vs. ‘mock; RIL-11272’; 9. ‘dehydration stress’ at ‘1 h’ vs. ‘control’ at ‘0 h’; 10. ‘dehydration stress’ at ‘12 h’ vs. ‘control’ at ‘0 h’; 11. ‘dehydration stress’ at ‘6 h’ vs. ‘control’ at ‘0 h’; 12. ‘salt stress’ at ‘1 h’ vs. ‘control’ at ‘0 h’; 13. ‘salt stress’ at ‘12 h’ vs. ‘control’ at ‘0 h’; 14. ‘salt stress’ at ‘6 h’ vs. ‘control’ at ‘0 h’; 15. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 0’; 16. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 12’; 17. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 16’; 18. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 20’; 19. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 4’; 20. ‘water deficit, 30% soil field capacity’ vs. ‘control, 70% soil field capacity’ at ‘Zeitgeber time 8’. Legend: red spheres indicate up-regulated GmLTP loci (*Log2FC > 1 and FDR < 0.05) in at least one of the analyzed treatments; ‡ in chromosomal order, with abbreviated nomenclature (without the Glyma epithet), and followed by their respective nsLTP categorization.
Antibiotics 12 00939 g009
Table
Table 1. Some previous studies involving nsLTP mining in plant omics data, including data sources, identification strategy, amount per species, and classes retrieved.
Table 1. Some previous studies involving nsLTP mining in plant omics data, including data sources, identification strategy, amount per species, and classes retrieved.
Analyzed SpeciesClassMining MethodologynsLTPs AmountnsLTPs ClassificationReferences
Oryza sativa (Os), Arabidopsis thaliana (At) Monocot/dicottBLASTn/BLASTnOs (52), At (49)1, 2, 3, 4, 5, 6, 7, 8, and YBoutrot et al. [17]
Solanum tuberosum (St), Solanum lycopersicum (Sh), Nicotiana tabacum (Nt), Nicotiana benthamiana (Nb), Capsicum annuum (Ca), and Petunia hybrida (Ph)DicotsBLASTnSt (28), Sl (28), Nt (33), Nb (17), Ca (19), Ph (10)1, 2, 3, 4, 5, 8, and 9Liu et al. [21]
Adiantum capillus-veneris (Ac-v), Marchantia polymorpha (Mp), Physcomitrella patens (Pp), Pinus taeda (Pt), Selaginella moellendorffii (Sm), and green algae (ga)Chlorophyta, bryophyta, dicots, and monocotstBLASTn/BLASTn/HMMMp (14), Pp (40), Sm (43), Ac-v (6), Pt (40), ga (0)1, 2, C, D, E, F, G, H, J, and K * Edstam et al. [1]
Lotus japonicusDicotBLAST241, 2, 3, 4, 5, 8, and 9Tapia et al. [22]
Brassica rapaDicotBLASTp631, 2, 3, 4, 5, 6, 8, 9, and 11Li et al. [23]
Zea maysMonocotBLASTp/HMM631, 2, C, D, and GWei e Zong et al. [5]
Gossypium arboreum (Ga), Gossypium raimondii (Gr), and Gossypium hirsutum (Gh)DicotBLASTpGa (51), Gr (47), Gh (91)1, 2, 3, 4, 5, 6, 8, and 9Li et al. [6]
Brassica oleraceaDicotHMM/BLASp891, 2, C, D, E, and GJi et al. [24]
Triticum aestivumMonocotBLAST Search1051 and 2Hairat et al. [25]
Triticum aestivumMonocottBLASTn4611, 2, C, D, and GKouidri et al. [26]
Hordeum vulgareDicotBLAST Search701, 2, C, D, and GZhang et al. [27]
Solanum lycopersicumDicotHMM641, 2, 3, 4, 10, and 11D’Agostino et al. [28]
Solanum tuberosumDicotBLASTp/tBLASTn/HMM831, 2, 4, 5, 7, 8, 12, and 13 **Li et al. [29]
Triticum aestivumMonocotBLASTp/HMM3301, 2, C, D, and GFang et al. [30]
Arachis duranensisDicotHMM641, 2, C, D, E, and GSong et al. [31]
Sesamum indicumDicot BLASTp/HMM521, 2, 3, 4, 5, 6, 8, 9, and 11Song et al. [32]
Hordeum vulgare (Hv) and highland barley (hb)DicotBLASTp/HMMHv (40), hb (35)1, 2, C, D, and GDuo et al. [33]
Chlamydomonas reinhardtii (Cr), Marchantia polymorpha (Mp), Physcomitrella patens (Pp), Selaginella moellendorffii (Sm), Zea mays (Zm), Sorghum bicolor (Sb), Oryza sativa (Os), Arabidopsis thaliana (At), Phaseolus vulgaris (Pv), Glycine max (Gm), Medicago truncatula (Mt), Trifolium pratense (Tp), Lotus japonicus (Lj), Lupinus albus (La), and Pisum sativum (Ps)Chlorophyta, bryophyta, dicots, and monocotsBLASTpCr (1), Mp (13), Pp (28), Sm (23), Zm (68), Sb (63), Os (73), At (82), Pv (77), Gm (120), Mt (95), Tp (85), Lj (72), La (87), Os (73)1, 2, C, D, E, and G *Fonseca-García et al. [34]
Brassica napusDicotBLASTx2461, 2, C, D, and GLiang et al. [35]
Helianthus annuusDicotBLASTp/HMM1011, 2, 3, and 4Vangelisti et al. [36]
Sorghum spontaneumDicotBLAST, TBLASTN71 and 2de Oliveira Silva et al. [37]
Brassica napusDicotHMM2381, 2, 3, 4, 5, 6, and 7 Xue et al. [38]
* All nsLTP groups found in the studied species pool; the nsLTP groups varied among the analyzed species. ** According to the similarity of the eight cysteine motif domains in amino acid sequences, Solanum tuberosum sequences were classified into two new types (12 and 13) by Li et al. [29].
Table
Table 2. Studied species and number of recovered nsLTPs from the three applied mining approaches.
Table 2. Studied species and number of recovered nsLTPs from the three applied mining approaches.
Plant CategorizationCommon NameGenome Version *2n Genome Size (Gb)Mining Method **Number of NR Loci
Higher ClassificationFamilySpeciesBLASTpRegExHMM
BryophytaMarchantiaceaeMarchantia polymorphaCommon liverwortMpolymorpha_320_v3.118~0.29 022121
PteridophytaPteridaceaeCeratopteris richardiiTriangle waterfernCrichardii_676_v2.178~11.25 0646565
SelaginellaceaeSelaginella moellendorffiiSpikemossSmoellendorffii_91_v1.016~0.1003636
GymnospermCupressaceaeThuja plicataWestern redcedarTplicata_572_v3.122~12.5412112112
AngiospermMalvaceaeGossypium hirsutumCottonGhirsutum_527_v2.152~2.436754218218
AsteraceaeLactuca sativaLettuceLsativa_467_v518~2.55544105105
EuphorbiaceaeManihot esculentaCassavaMesculenta_671_v8.136~0.729289898
PhrymaceaeMimulus guttatusMonkeyflowerMguttatusTOL_551_v5.028~0.41422114114
SalicaceaePopulus trichocarpaBlack cottonwoodPtrichocarpa_533_v4.119~0.525228585
BrassicaceaeSinapsis albaWhite mustardSalba_584_v3.124~0.52364189189
SolanaceaeSolanum tuberosumPotatoStuberosum_686_v6.148~0.83524105105
AmaranthaceaeSpinacia oleraceaSpinachSoleracea_575_Spov312~0.9684343
* Phytozome genome version used to infer the conceptual proteome; ** number of nsLTP sequences with validated eight-cysteine domain (8CM). Legend: RegEX (regular expression); HMM (Hidden Markov Model); Gb (gigabase); NR (non-redundant).
Table
Table 3. nsLTP quantification by categories in the 12 analyzed plant genomes.
Table 3. nsLTP quantification by categories in the 12 analyzed plant genomes.
Plant CladePlant SpeciesnsLTP CategoryTotal 2
LTP1LTP2LTPdLTPgLTPcUnknown 1Unknown 2Unknown 3Unknown 4
BryophytaMp111040300221
PteridophytaSm0112611000636
Cr32002514021065
GymnospermTp320941242220112
AngiospermSo1407171013043
Pt173153321113085
Me245112911225098
Ls26107311210171105
St414728024190105
Mg40410282510150114
As251531702212320189
Gh30173476535480218
Total 128260153365224747196191191
Legend: Marchantia polymorpha (Mp), Ceratopteris richardii (Cr), Selaginella moellendorffii (Sm), Thuja plicata (Tp), Glycine max (Gm), Gossypium hirsutum (Gh), Lactuca sativa (Ls), Manihot esculenta (Me), Mimulus guttatus (Mg), Populus trichocarpa (Pt), Sinapis alba (Sa), Solanum tuberosum (St), Spinacea oleracea (So), and Unk. (unknown); Total 1 (total of loci by category); Total 2 (total of loci encoding the nsLTPs in each plant genome).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Santos-Silva, C.A.d.; Ferreira-Neto, J.R.C.; Amador, V.C.; Bezerra-Neto, J.P.; Vilela, L.M.B.; Binneck, E.; Rêgo, M.d.S.; da Silva, M.D.; Mangueira de Melo, A.L.T.; da Silva, R.H.; et al. From Gene to Transcript and Peptide: A Deep Overview on Non-Specific Lipid Transfer Proteins (nsLTPs). Antibiotics 2023, 12, 939. https://doi.org/10.3390/antibiotics12050939

AMA Style

Santos-Silva CAd, Ferreira-Neto JRC, Amador VC, Bezerra-Neto JP, Vilela LMB, Binneck E, Rêgo MdS, da Silva MD, Mangueira de Melo ALT, da Silva RH, et al. From Gene to Transcript and Peptide: A Deep Overview on Non-Specific Lipid Transfer Proteins (nsLTPs). Antibiotics. 2023; 12(5):939. https://doi.org/10.3390/antibiotics12050939

Chicago/Turabian Style

Santos-Silva, Carlos André dos, José Ribamar Costa Ferreira-Neto, Vinícius Costa Amador, João Pacífico Bezerra-Neto, Lívia Maria Batista Vilela, Eliseu Binneck, Mireli de Santana Rêgo, Manassés Daniel da Silva, Ana Luiza Trajano Mangueira de Melo, Rahisa Helena da Silva, and et al. 2023. "From Gene to Transcript and Peptide: A Deep Overview on Non-Specific Lipid Transfer Proteins (nsLTPs)" Antibiotics 12, no. 5: 939. https://doi.org/10.3390/antibiotics12050939

APA Style

Santos-Silva, C. A. d., Ferreira-Neto, J. R. C., Amador, V. C., Bezerra-Neto, J. P., Vilela, L. M. B., Binneck, E., Rêgo, M. d. S., da Silva, M. D., Mangueira de Melo, A. L. T., da Silva, R. H., & Benko-Iseppon, A. M. (2023). From Gene to Transcript and Peptide: A Deep Overview on Non-Specific Lipid Transfer Proteins (nsLTPs). Antibiotics, 12(5), 939. https://doi.org/10.3390/antibiotics12050939

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop