- Research article
- Open Access
Non-specific lipid transfer proteins in maize
BMC Plant Biology volume 14, Article number: 281 (2014)
In plant, non-specific lipid transfer proteins (nsLTPs) are small, basic proteins that have been reported to be involved in numerous biological processes such as transfer of phospholipids, reproductive development, pathogen defence and abiotic stress response. To date, only a tiny fraction of plant nsLTPs have been functionally identified, and even fewer have been identified in maize [Zea mays (Zm)].
In this study, we carried out a genome-wide analysis of nsLTP gene family in maize and identified 63 nsLTP genes, which can be divided into five types (1, 2, C, D and G). Similar intron/exon structural patterns were observed in the same type, strongly supporting their close evolutionary relationship. Gene duplication analysis indicated that both tandem and segmental duplication contribute to the diversification of this gene family. Additionally, the three-dimensional structures of representative nsLTPs were studied with homology modeling to understand their molecular functions. Gene ontology analysis was performed to obtain clues about biological function of the maize nsLTPs (ZmLTPs). The analyses of putative upstream regulatory elements showed both shared and distinct transcriptional regulation motifs of ZmLTPs, further indicating that ZmLTPs may play roles in diverse biological processes. The dynamic expression patterns of ZmLTPs family across the different developmental stages showed that several of them exhibit tissue-specific expression, indicative of their important roles in maize life cycle. Furthermore, we focused on the roles of maize nsLTPs in biotic and abiotic stress responses. Our analyses demonstrated that some ZmLTPs exhibited a delayed expression pattern after the infection of Ustilago maydis and differentially expressed under drought, salt and cold stresses, and these may be a great help for further studies to improve the stress resistance and tolerance in maize breeding.
Our results provide new insights into the phylogenetic relationships and characteristic functions of maize nsLTPs and will be useful in studies aimed at revealing the global regulatory network in maize development and stress responses, thereby contributing to the maize molecular breeding with enhanced quality traits.
As its name implies, plant lipid transfer protein (LTPs) were termed because of their functions that transfer phospholipids and fatty acids between membranes in vitro. They were also named non-specific LTPs (nsLTPs) due to the character of non-specific binding to different lipids. Plant nsLTPs are small, basic proteins, usually about 6.5 ~ 10.5 kDa in size, characterized by an eight cysteine motif (8CM) backbone with the general form C-Xn-C-Xn-CC-CXC-Xn-C-Xn-C . Almost all nsLTPs carry an N-terminal signal peptide in their nascent polypeptides. Thus, these proteins are likely secreted to the cell exterior for functioning. Many nsLTPs also possess a sequence for the post-translational addition of a glycosylphosphatidylinositol (GPI)-anchor, which attaches the protein to the exterior side of the plasma membrane . The conserved features of plant nsLTPs include four defined disulfide bonds formed by eight Cys residues. Furthermore, the crystal structures of plant nLTPs were comprised of four or five alpha helices (α-helices), with a central hydrophobic cavity where the lipid binding takes place .
Based on the molecular weight (Mw) of the mature proteins, plant nsLTPs can be classified into two main types, nsLTP1 (9 kDa) and nsLTP2 (7 kDa) . According to the sequence similarity, Boutrot et al. classified 49 out of the 52 rice nsLTPs and 45 out of the 49 Arabidopsis nsLTPs into nine types (I, II, III, IV, V, VI, VII, VIII, and Y) . Recently, nsLTPs have been categorized into four major and several minor types by sequence similarity, intron position and spacing between the cysteine residues, as well as a potential glycophosphatidylinositol modification site within their encoded proteins . In vitro studies showed that plant nsLTPs have the ability to facilitate the inter-membrane exchange and transfer of various amphiphilic molecules including phospholipids, glycolipids, steroids, acyl-CoAs and fatty acids. Several structures of plant nsLTPs have been resolved by X-ray and NMR spectroscopic techniques ,,. The typical fold of nsLTP1 is characterized by four to five α-helices connected by four disulfide bridges, partly wrapped by a long C-terminal segment. The overall structure delimits a large central hydrophobic cavity, where the alkyl moiety of lipids is inserted into. The backbone folds of nsLTP1 and nsLTP2 show structural similarities, however, drastic changes in their central hydrophobic cavity. Compared with these two earliest nsLTP types, the three-dimensional model of Arabidopsis DIR1 (AtDIR1) follows the general nsLTPs fold with five α-helices stabilized by four disulfide bonds around a central tunnel-shaped cavity. It was reported that most of the putative functions of nsLTPs are related to their ability to bind lipids in their hydrophobic cavity -.
In the ensuing years, some members of nsLTPs have been functionally identified in plant species, including the involvement of cuticular waxes and cutin syntheses . A recent study demonstrated that LTPg1, a Type G nsLTP from Arabidopsis, contribute either directly or indirectly to cuticular lipid accumulation . In addition, nsLTPs are expressed in diverse organs and cells, including callus, germinating and maturing seeds, leaves, roots, stems, ovaries, anthers, and pollens ,. The localization of nsLTP transcripts in anthers has been well reported in Arabidopsis and rice, and abundant Type III nsLTPs (also termed Type C nsLTPs) were expressed specifically in the anther tapetum, with levels peaking at the developmental stage of maximal pollen-wall exine synthesis . A lipid transfer protein, OsC6, is widely distributed in anther tissues such as the tapetal cytoplasm, required for postmeiotic anther development in rice . With respect to plant defense, nsLTPs are also recognized to be pathogenesis-related proteins and constitute the PR-14 family . Some of them were demonstrated to share structural similarities with oomycetous elicitins and act as competitors with elicitins for specific receptors on membrane ,. In addition, nsLTPs participate in long-distance signaling during pathogen defense . Studies on purified nsLTPs confirmed their roles in abiotic stress tolerance . Transcript levels of nsLTPs increased in response to drought, salt and cold stresses in many cases . Stabilization of membranes, cuticle deposition and/or changes in cell wall organization have been claimed as their putative roles in the responses to these stress factors. In tobacco, nsLTP genes up-regulated during drought-induced cuticular wax deposition . Overexpression of the pepper CALTP1 gene in Arabidopsis increased its tolerance to NaCl and drought stresses at various vegetative growth stages . Arabidopsis LTP3 act as a target of MYB96 to be involved in plant tolerance to freezing and drought stresses .
In previous studies, only a small portion of nsLTPs from maize have been well characterized. These include, MZm3-3, which is expressed specifically in the tapetum during male gametogenesis ; Zm-LTP (ZmLTP1.2; GRMZM2G010868), which binds to calmodulin (CaM) in a Ca2+-independent manner ; ZmLTP3 (ZmLTP1.1; GRMZM2G126397), which is induced by mannitol, salt and SA treatments ; BETL-9 (ZmLTPd6; GRMZM2G087413), which is transcribed in the outer surface of the developing endosperm . As the genome of inbred line B73 has been sequenced completely, it is the time to initiate functional annotation and to perform more comprehensive analysises of maize nsLTP gene family. To date, genome-wide overview of the maize nsLTP gene family was yet to be reported. Therefore, as the first step to elucidate the functions of ZmLTPs, a genome-wide study for this gene family is necessary. In this study, we identified 63 genes encoding 77 putative nsLTPs in the maize genome that can be classified into five types. Detailed analyses including primary sequence, phylogenetic relationships, gene structure, chromosome location, gene duplication and divergence, three-dimensional structure, gene ontology, promoter and expression profilings were performed. Based on microarray, RNA-sequencing and real-time quantitative PCR (qRT-PCR), we analyzed the expression patterns of ZmLTP genes in different tissues and developmental stages, and took a further step towards understanding plant responses to biotic and abiotic stresses. Our analysis provides novel insights on ZmLTP gene family to support further functional research on nsLTP gene family in plants, particularly the proteins that may have important functions in response to biotic and abiotic stresses.
Results and discussion
Identification of nsLTP members in maize and sorghum
In order to identify the complete and non-redundant nsLTPs in maize, an accurate scan of the maize proteome was performed. Initially, 130 potential ZmLTPs were identified in the maize genome following the removal of those sequences with an incomplete 8CM domain. Then, each of the deduced protein sequences was manually assessed through the analysis of the cysteine residue patterns. Subsequently, 40 proline-rich or hybrid proline-rich proteins, which were characterized by a high proportion of proline, histidine and glycine residues in the sequence comprised between the signal peptide and the 8CM, were removed (Additional file 1: Table S1). Next, two α-amylase/trypsin inhibitors, three prolamin storage proteins and three 2S albumin storage proteins were also discarded. Additionally, five transcript forms lacking the GPI-anchor signal peptide (GRMZM2G174680_P02, GRMZM2G083725_P02, GRMZM2G006047_P01, GRMZM2G005991_P02, GRMZM2G116167_P02) were not taken into consideration. As a result, 77 proteins were confirmed as maize nsLTPs which were encoded by 63 nsLTP genes (Additional file 2: Table S2). Additionally, the same approach was used to identify nsLTPs in sorghum, which resulted in the identification of 58 nsLTP genes (Additional file 3: Table S3 and Additional file 4: Table S4).
Sequence analysis and classification of putative maize nsLTPs
Previously, Edstam et al. have proposed that the plant nsLTPs can be divided into four major and several minor types according to sequence similarity, intron position and spacing between the cysteine residues . In one of the major types, Type G, the transcripts encode a C-terminal signal sequence in addition to the N-terminal one, leading to a post-translational modification where a GPI-anchor is added to the protein. The classification of the identified type G ZmLTPs was based on the presence of a GPI modification site, as well.. In the second round of classification, the remaining sequences were grouped according to the identity matrix calculated from the multiple sequence alignments. When compared with the classification proposed by Edstam et al., we found that 73 out of the 77 ZmLTPs could be classified into five types (1, 2, C, D and G). Among the nsLTPs that were present in four flowering plants (maize, sorghum, rice and Arabidopsis), the variation in the number of nsLTPs across types and species was detected (Additional file 5: Table S5). In type G, 26, 24, 27, and 29 nsLTPs were found in maize, sorghum, rice and Arabidopsis, respectively. Whereas, only 2 nsLTP genes were distributed to Type E in Arabidopsis. Noticeably, a similar number of nsLTP genes were found in Type 1 and Type 2. It is thus conceivable that this difference may be due to gene duplication and loss. The main characteristic of plant nsLTPs is the presence of eight cysteine residues in a highly conserved region (C-Xn-C-Xn-CC-Xn-CXC-Xn-C-Xn-C). In order to establish a specific 8CM consensus for each obtained nsLTP type, we conducted a multiple sequence alignment using the 8CMs from the 77 ZmLTPs (Additional file 6: Figure S1). The amino acid sequence alignment of the 8CMs of ZmLTPs reveals a variable number of inter-cysteine amino acid residues (Table 1).
In this study, all the characteristics of the 77 ZmLTPs were summarized in Additional file 2: Table S2. 96% of the nsLTP precursors were initially synthesized with a signal peptide of 16-46 amino acids. The putative subcellular targeting of the 77 ZmLTP pre-protein sequences was analyzed. As expected, most of the proteins are predicted to be secreted except for ZmLTPg21, which have been predicted to be cytoplasm protein (Additional file 2: Table S2). At the pre-protein level, the ZmLTPd1 and ZmLTPd4 deduced proteins are identical, as are the ZmLTPd2 and ZmLTPd10 deduced proteins. After cleavage of their signal peptide, the ZmLTPg20.1 and ZmLTPg20.2 mature proteins are identical. Therefore, before potential post-translational modifications, the 63 ZmLTP genes encode 74 different mature proteins. 29 ZmLTPs were found to have one or more phosphorylation sites and most of the phosphorylated sites were located on serine residues at the C-terminal (Additional file 7: Table S6). To clearly understand the sequence characteristics of ZmLTPs, we further analyzed the pI (isoelectric point) values, Mw values, and CXC motifs of all available ZmLTPs. As shown in Additional file 2: Table S2, maize nsLTPs are small and their molecular masses usually range from 6,854 Da to 11,107 Da except for Type G. Type 1 and Type D nsLTPs were mostly 9 kDa proteins and Type 2 and Type C nsLTPs were 7 kDa proteins. The Mw value of Type G nsLTPs was much higher than that of other types due to the presence of supernumerary amino acid residues located at the C-terminal of the deduced mature proteins. Judging from the pI value, Type 1, 2, C and D nsLTPs are mostly alkaline proteins. The majority of Type G nsLTPs are weakly alkaline or acidic (Additional file 2: Table S2). The average molecular mass and the theoretical pI are 8,963 Da and 9.27 pI, respectively. As for the CXC motif, most residues at the X position in Type 1 nsLTPs are hydrophilic, while in Type 2, C, D and G, the X position is usually occupied by a hydrophobic residue (Additional file 6: Figure S1). These conserved hydrophobic or hydrophilic residues may play significant roles in the biological functions of ZmLTPs .
Phylogenetic analysis of the maize, sorghum, rice and Arabidopsis nsLTPs
To analyze the phylogenetic relationship of the nsLTPs among maize, sorghum, rice and Arabidopsis, 274 nsLTPs from these four species were analyzed (Additional file 8: Figure S2). We performed a multiple sequence alignment of the 8CM domain sequences from maize, sorghum, rice and Arabidopsis and then generated a phylogenetic tree by the neighbor-joining method. Previously, Edstam et al. has divided the plant nsLTPs into ten types . On the basis of the comparison between the previous dataset and ours, the six groups classified here were in agreement with the Type 1, 2, C, D, E and G of nsLTPs. As shown in Additional file 8: Figure S2, members in Type 1 and 2 formed specific clades, indicating that the genes in these major nsLTP types share a common ancestor. The sequences from the minor Type C and E also formed separated clades. Since relatively high bootstrap values were observed in the internal branches of Type 1 and Type C, it clearly showed the derivation of statistically reliable pairs of possible homologous proteins sharing similar origin from a common ancestor. It was also worth noting that Type E nsLTPs seem to be monophyletic, and this may mean that the monocotyledon plants discarded these genes during the evolutionary divergence between monocots and dicots.
Intron-exon structure of the maize nsLTP gene family
As a kind of evolutionary relic, the intron-exon arrangement carries the imprint of the evolution of a gene family. Investigation of ZmLTP gene structures revealed low diverse distribution of intronic regions amid the exonic sequences, and 31 ZmLTP genes (36 Alternative splicing forms) were predicted to be interrupted by 1-2 introns positioned 4 to 67 bp downstream of the codon encoding the eighth cysteine in the 8CM (Additional file 9: Figure S3). Alternative splicing forms of the same gene usually had similar intron–exon structures, indicating that the distinct proteins from a single transcript may share similar functions. Additionally, it is interesting to find a similar exon/intron pattern in each group. For instance, the nsLTP genes in Type 1, 2 and C lacked an intron, while all members in Type G contained 1-3 introns. Except for ZmLTPd6, ZmLTPd13, ZmLTPd14 and ZmLTPd15, no intron was detected in the coding regions of Type D genes. Our comparative analysis with the gene structures of AtLTPs and OsLTPs indicated that the exon/intron structures of ZmLTPs are similar to those of Arabidopsis and rice (Additional file 10: Figure S4).
Chromosomal localization and gene duplication of ZmLTP genes
In silico mapping of the gene loci showed that 63 ZmLTPs were unevenly assigned to maize ten chromosomes (Figure 1A). Chromosome 1 contained the maximum number of ZmLTPs (11), while the minimum number (2) were presented on Chromosome 9. However, several chromosomes lacked ZmLTPs in some regions, such as the long arm of chromosome 5 as well as the short arms of Chromosome 3, 6 and 9. The exact position (in bp) of each nsLTP on maize chromosome is given in Additional file 2: Table S2.
Gene duplication is generally believed to be a major driving force in evolutionary innovation, giving rise to genomic complexity. Maize originated from an ancient allotetraploid and has undergone several rounds of whole-genome duplication events during its gene evolution ,. Segmental and tandem duplications are known to be the main causes leading to gene family expansions. In this study, a total of seven tandem repeats (ZmLTPg6/ZmLTPg7/ZmLTPg8; ZmLTP2.7/ZmLTP2.8; ZmLTPg20/ZmLTPg21; ZmLTPd10/ZmLTPd11; ZmLTPd12/ZmLTPd13; ZmLTP1.6/ZmLTP1.7 and ZmLTPd15/ZmLTPd16) were identified in the maize genome. One significant cluster of three Type G genes was found on Chromosome 1. What's more, six direct repeat tandems were identified on chromosome 6, 7, 8 and 10. Genes in the same cluster were closely related to one another. For instance, ZmLTPd15 and ZmLTPd16 protein sequences shared 95% similarity and the two genes were physically located next to each other in the same chromosomal region. In addition, five sister pairs (ZmLTPg6/ZmLTPg16, ZmLTPg11/ZmLTPx2, ZmLTPg13/ZmLTPg21, ZmLTPd7/ZmLTPd13 and ZmLTP2.7/ZmLTP2.9) appeared to be generated from segmental duplication events due to their positions on the same duplicated gene blocks (Figure 1A and Additional file 11: Table S7). Furthermore, we analyzed the evolution of nsLTP genes among maize, sorghum and rice (Figure 1B and Additional file 11: Table S7). 22 out of 63 maize nsLTP genes had collinear genes in rice, while 35 had syntenic members in sorghum. According to this analysis, two paired nsLTP genes (ZmLTPg21/OsLTPg17 and ZmLTPg13/OsLTPg17) were located in genomic regions with synteny between the maize and rice genomes. Seven paired nsLTP genes (ZmLTPg16/SbLTPx1, ZmLTPg11/SbLTPg18, ZmLTPx2/SbLTPg18, ZmLTP2.9/SbLTP2.5, ZmLTPg21/SbLTPg10, ZmLTPg13/SbLTPg10 and ZmLTPg13/SbLTPg4) were located in genomic regions with synteny between the maize and sorghum genomes, indicating that these genes may be derived from a common ancestor. Interestingly, one of the sister pairs (ZmLTPg6/ZmLTPg16) had a unique, syntenic nsLTPs in sorghum, which indicated that they might be generated from segmental duplication after the divergence of maize and sorghum.
Divergence rate of the maize nsLTP genes
The ratio (Ka/Ks) of nonsynonymous substitution rate (Ka) versus synonymous substitution rate (Ks) is widely used as an indicator of selective pressure . As shown in Table 2, the Ka/Ks ratios of 5 duplicated gene pairs were less than 1, indicating that they seemed to be under purifying selection; however, the Ka/Ks ratios of 9 duplicated gene pairs were more than 1, suggesting that they subjected to positive selection. Positive selection was thought to be one of the major forces for the emergency of new motifs/functions in protein after gene duplication, and the divergence of the duplicated genes was driven by positive selection . Based on the substitution rates analysis, we found that 64.29% of the duplicated gene pairs had strong positive selection pressure, implying that positive selection contributed to further gene diversification in the ZmLTP family. In addition, the divergent times between the duplicated gene pairs were analyzed (Table 2). The duplication events for the 5 segmental duplications were estimated to have occurred approximately between 9.43 and 34.06 million years ago (Mya), while the duplication events for the 9 tandem duplications occurred approximately between 9.13 and 50.46 Mya. Fossil data along with the phylogenetic studies estimated that the different grass families diverged from a common ancestor 50 to 70 Mya . The divergence times on the nodes of the tree were estimated with the non-parametric rate smoothing method of Sanderson, assuming that maize and rice diverged 50 Mya . Duplication events for one tandem duplicated gene pair (ZmLTPg20/ZmLTPg21) occurred around 50.46 Mya, after origin of grasses and before divergence of rice and maize (within last 70 to 50 million years). The ancestor of maize and sorghum diverged about 12 Mya and subsequently a whole genome triplication (WGT) event occurred in the maize approximately 5 Mya . Duplication events for 10 gene pairs (ZmLTPg11/ZmLTPx2, ZmLTPg13/ZmLTPg21, ZmLTPd10/ZmLTPd11, ZmLTP2.7/ZmLTP2.9, ZmLTPd7/ZmLTPd13, ZmLTP1.6/ZmLTP1.7, ZmLTP2.7/ZmLTP2.8, ZmLTPd12/ZmLTPd13, ZmLTPg6/ZmLTPg7 and ZmLTPg6/ZmLTPg8) occurred within last 39.91 to 18.12 million years, after divergence of rice and maize, but before maize and sorghum were separated from each other. Duplication events for the other two gene pairs (ZmLTPg6/ZmLTPg16 and ZmLTPd15/ZmLTPd16) occurred around last 9 million years, after maize and sorghum were separated, and before a whole genome duplication (WGD) event occurred in the maize. Therefore, evolutionary origin of the maize nsLTP genes might undertake three evolutionary stages.
The features of the three-dimensional structures of major nsLTPs
In order to better understand the non-specific binding between the hydrophobic ligands and ZmLTPs, ZmLTP1.2, ZmLTP2.1 and ZmLTPd5 were selected as representative sequences of Type 1, 2, and D for structural modeling (Additional file 2: Table S2). The crystal structures of maize ZmLTP1.6 (PDB ID: 1FK7) , wheat nsLTP2 (PDB ID: 1TUK) , AtDIR1 (PDB ID: 2RKN)  were selected as templates for structural modeling based on searches against the PDB using the Basic Local Alignment Search Tool (BLAST) at NCBI (http://blast.ncbi.nlm.nih.gov/) with the target Type 1, 2 and D protein sequences as baits.
Based on the ZmLTP1.6, nsLTP2 and AtDIR1 structures, the amino acids 29-121 of ZmLTP1.2, 32-98 of ZmLTP2.1, and 37-115 of ZmLTPd5 could be modeled with the sequence identities of 88.17%, 62.69%, and 51.95%, respectively. Our structural analysis showed that ZmLTP1.2 and ZmLTP1.6 have typical features of plant nsLTPs, including two conserved pentapeptides, T-T/A-A-D-R (positions 46-50) and P-Y-T-I-S (positions 87-91). It has been reported that these two consensus pentapeptides (T/S-X-X-D-R/K and P-Y-X-I-S) were important in catalysis or binding . ZmLTPd5 and AtDIR1 contain high proline content (positions 25-32) that is not highlighted in nsLTP2 and ZmLTP2.1. It is noteworthy that proline-rich regions are involved in protein-protein interactions . As aforementioned, X is a hydrophilic residue in the CXC motif of Type 1 nsLTPs. However, a hydrophobic residue was found at the X position in Type 2 and Type D nsLTPs. As illustrated in Figure 2A, in the CXC motif, asparagines between the two cysteines in ZmLTP1.2 and ZmLTP1.6 are replaced by a hydrophobic amino acid, phenylalanine, in ZmLTP2.1 and nsLTP2, whereas in ZmLTPd5 and AtDIR1, the X position are occupied by a hydrophobic residue, leucine. The hydrophobic residues in the CXC motif of ZmLTP2.1, nsLTP2, ZmLTPd5 and AtDIR1 are buried inside the molecule, whereas the hydrophilic residue of ZmLTP1.2 and ZmLTP1.6 are at the surface (Figure 2B, C and D). In the template protein ZmLTP1.6, the Cys residues 1-6, 2-3, 4-7, and 5-8 are paired, whereas the Cys residues 1-5, 2-3, 4-7, and 6-8 are paired in the template protein AtDIR1, similarly to the Type 2 nsLTPs. These observations suggested that the central residue of the CXC motif may govern the cysteine pairing and influence the overall fold of the protein. The crystal structure homology models were validated in terms of stereochemical quality by Ramachandran plot and the ERRAT prediction (Figure 2E).
From the superposition structure, we find that ZmLTP1.2 and ZmLTP1.6 share a common structural fold stabilized by four disulfide bonds, and the four prominent α-helices are packed against a flexible C-terminal arm formed by a series of turns, which forms a cap over the hydrophobic cavity (Figure 2B). The three-dimensional structure of ZmLTP2.1 showed that two adjacent hydrophobic cavities seem able to expand the hydrophobic cavity to accommodate the alkyl ligand (Figure 2C). The structure of ZmLTPd5 is in accordance with AtDIR in the aspects of the order and orientation. Their α-helix topology is well conserved and the four cysteine bridges superimpose well, showing high similarity in the three-dimensional structure (Figure 2D). Interestingly, the Arabidopsis mutant defective (dir1-1) was compromised in the production or transmission of an essential mobile signal to promote long-distance signaling . Whether other Type D nsLTPs have the functions related to systemic resistance signalling remain to be further investigated.
Gene ontology analysis of the maize nsLTPs
In order to comprehend the unique aspects of maize nsLTPs, Gene Ontology (GO) enrichment analysis of the functional significance was performed. Out of 77 ZmLTP proteins, annotation could not be performed for 24 proteins. Besides, 53 ZmLTPs were defined in 22 significant GO terms (Additional file 12: Table S8). The analysis showed that 53 ZmLTPs were separated into two main categories (biological process and molecular function), which included 18 and 4 significant GO terms, respectively (Figure 3A). For the enriched biological processes, the common categories are "response to stimulus" followed by "response to abiotic stimulus" and "response to stress", "multicellular organismal process", "biological regulation" and "establishment of localization". Noteworthy, about 47 ZmLTPs were shown to participate in "lipid transport" (GO: 0006869), which is concordance with the molecular role of nsLTP in transporting hydrophobic molecules in vitro, suggesting that ZmLTPs play an important role in carrying membrane components to the growing site of elongating cells. The enriched GO terms also include the parent term "response to abiotic stimulus" with enriched children term "response to cold", and 22.64% (12 of 53) ZmLTPs were exhibited to participate in "response to stress stimulus" (particularly in cold). This highlights the putative association of nsLTP proteins in stress tolerance behavior of maize. In case of molecular functions, about 88.68% (47 of 53) ZmLTPs were shown to participate in "lipid binding" (GO: 0008289). One of the interesting observations was that some significantly enriched terms are involved in "water binding" and "ice binding". To summarize, the GO analysis indicated that ZmLTPs may be involved in diverse biological processes.
Key cis-elements within the promoter regions of ZmLTP genes
In-silico analysis of 1 kb upstream region (from translation start site) of 63 maize nsLTP genes revealed the presence of various regulatory elements, which are associated with development, abiotic or biotic stress signaling, and hormone signaling (Additional file 13: Table S9). To avoid the biases in analysis, only cis-elements identified in at least ten different genes were considered. As shown in Figure 3B, except for TATA-box and CAAT-box, G-box is the most frequently found type of cis-elements, and the regulatory element involving light responsiveness seems to be enriched in ZmLTPs. Nearly 80% of ZmLTPs contained Skn-1_motifs which might account for endosperm expression. Stress-responsive cis-regulatory elements identified in this study included fungal elicitor responsive element (Box-W1), defense and stress responsive element (TC-rich repeats), heat stress responsive element (HSE), low temperature responsive element (LTR), anaerobic-response element (ARE and GC-motif), MYB binding site (MBS) involving in drought-inducibility and phytohormone-responsive elements, like auxin responsive element (TGA-element and AuxRR-core), methyl jasmonate responsive element (CGTCA-motif and TGACG-motif), salicylic acid responsive element (TCA-element), abscisic acid responsive element (ABRE) and gibberellic acid-responsive element (GARE). With a few exceptions, ZmLTPs contained circadian elements, which may be responsible for its distinct diurnal expression pattern. Other regulatory elements, such as O2-site (involved in zein metabolism regulation), DRE (involved in responses to dehydration, low temperatures and salt), CAT-box (related to meristem expression) and CCGTCC-box (related to meristem specific activation) were also presented. Especially, the promoters of 39 ZmLTPs contained MBS elements ranging from 1 to 5 copies (Additional file 13: Table S9), indicating the important role of MYB transcription factors in regulating ZmLTPs. Some recent reports have shown that nsLTP can act as targets of MYB to regulate plant tolerance to freezing and drought stresses ,, suggested that the ZmLTP members participate in some abiotic stress signaling.
The diverse roles of nsLTPs in maize development
As for multigene family, the analysis of gene expression patterns often provides useful clues to decide genes function. To investigate the temporal and spatial expression patterns of the nsLTP genes in the maize life cycle, hierarchical clustering was performed to visualize a global transcription profile of the ZmLTP genes across the 11 organs during diverse developmental stages. As depicted in Figure 4, the heatmap can be apparently divided into three clusters. Cluster I contains 21 members (excluding the E2 enzyme as the reference). Cluster II has 9 members, and 45 members belong to Cluster III. Genes in Cluster I obviously have relatively high expression levels, with the mean log-signal values of each gene ranging from 9.0 to 14.9 (Additional file 14: Table S10). Conversely, Cluster III contains those genes with relatively low expression levels and the average log-signal value range is between 5.2 and 9.2. In Cluster II, genes are expressed at a moderate level with the mean log-signal values ranging from 6.1 to 8.5. To elucidate putative differentially expressed genes in specific organs or stages, the coefficient of variation (CV value; CV = S/Xmean, where S represents the standard deviation and Xmean describes the mean expression level of a gene across all the tissues) of each gene in the three clusters was calculated (Additional file 14: Table S10). A house keeping gene in maize, which encodes a ubiquitin-conjugating enzyme (E2), was used as internal reference in our expression analysis . The results showed a huge variation among all the genes with the CV values ranging from 1.01% to 52.97%. Cluster I has the least expression variability from 1.01% to 26.80%, indicating a most stable expression pattern relative to other ZmLTP genes. ZmLTP1.1, the most changeable genes in cluster I, whose homolog in rice, Photoperiod-sensitive dwarf 1 (Psd1, Os01g60740), was previous demonstrated to encode a lipid transfer protein that may participate in regulation of plant cell division and elongation . One gene (ZmLTP1.2), which can bind to calmodulin in a Ca2+-independent manner , displayed high expression levels at nearly all of the maize organs and/or stages of development analyzed except for root, immature cob and endosperm. A total of 12 nsLTP protein-encoding genes were found to be highly expressed in anthers. 7 ZmLTPs were found to be highly expressed in leaf at different developmental stages, whereas, 5 ZmLTPs specially accumulated at V5_Base of stage-2 Leaf, V7_Base of stage-2 Leaf, V9_Immature Leaves, V9_Eleventh Leaf and V9_Thirteenth Leaf. Cluster II showed the highest expression variability, ranging from 11.95% to 52.97%. There are 8 genes with CV values more than 15% in Cluster II, among which 3 genes (ZmLTPd14, ZmLTP2.3 and ZmLTPd6) have CV values more than 45% and display transcript accumulation at late seed developmental stages (10 to 24 DAP) and endosperm. It has been reported that BETL9 (ZmLTPd6) has its expression restricted to the endosperm transfer cells, which share their position as entrance gates to the developing seed for nutrients coming from the maternal tissues . There is an increasing evidence that this protein might be involved in the defense of the developing seed against mother plant-borne pathogens ,. Genes in Cluster III showed an apparent fluctuation, with the CV values ranging from 4.82% to 31.83%. In this context, Cluster III can be further divided into two subclusters, with genes in the first subcluster having higher signal intensity values than those in the second subcluster. In the first subcluster, four ZmLTPs (ZmLTP2.6, ZmLTP2.7, ZmLTPd1 and ZmLTPd4) displayed high expression levels during root development and this finding suggested that ZmLTPs may participate in maize root development. The organ-specific expression dynamics revealed the distinct expression patterns of ZmLTP genes throughout the entire life cycle in maize. For example, two ZmLTP genes, ZmLTPg20 and ZmLTP2.5, were specifically expressed in the anthers and coleoptile, respectively. There are 6 genes (ZmLTP1.4, ZmLTP1.5, ZmLTP2.9, ZmLTPc1, ZmLTPc2 and ZmLTPd9), with CV values among 15.54% to 23.92% in the second subcluster, distinctly increasing their expression levels in meiotic tassel. Among them, ZmLTPc2 is highly homologous to MZm3-3, which is expressed specifically in the tapetum during male gametogenesis of maize . Furthermore, four nsLTP genes (Os09g35700, Os01g49650, Os01g12020 and Os08g43290) were reported to be expressed in rice anther, which are putatively related to exine synthesis in sporophytic anthers of progressive developmental stages . Interestingly, these four rice nsLTPs were the homologs of ZmLTPc2, ZmLTP2.9, ZmLTP1.4 and ZmLTPc1, respectively. Synthesis of lipidic components in anthers, including the pollen exine, is essential for plant male reproductive development. Therefore, this finding further suggested that nsLTPs play significance roles in male reproductive development.
Pathogen responses of ZmLTPsto the fungal infection
Plant nsLTPs are capable of inhibiting bacterial and fungal pathogens and are, therefore, thought to play an important role in plant defense . To discover the ZmLTP genes involved in maize pathogen response, microarray data collected after the treatments of 12, 24 hours, 2, 4, 4.5, 8 Days of Ustilago maydis (U. maydis) infection were analyzed. U. maydis is a ubiquitous pathogen of maize which depends on living tissue for proliferation and development. A total of 41 probe sets found on the maize 18 k GeneChip could be assigned to 36 different ZmLTP genes (Additional file 15: Table S11). The log2 ratio values (Treated/Control) were illustrated by a heatmap, showing the fold change in expression of each ZmLTP. As shown in Additional file 16: Figure S5, most of the ZmLTPs exhibited a delayed expression pattern after the infection of U. Maydis. 23 ZmLTPs distinctly increased their expression levels and accumulated mostly on the 4 and 8 Days post infection, especially of ZmLTPg1, ZmLTP1.1 and ZmLTP1.7, suggesting that these genes might participate in the pathogen response. Interestingly, some ZmLTP genes show up-regulated expression levels over time, especially of ZmLTP2.8, demonstrating a similar response pattern to this kind of fungus infection. Besides, the expression levels of ZmLTPd2 and ZmLTP2.1 show a rapid increase after infection, suggesting their high sensitivity to U. maydis infection. One nsLTP from Arabidopsis (AtDIR1) has been suggested to be involved in long-distance signaling during pathogen defense . Therefore, it is reasonable to find a transcript accumulation of the ZmLTPd5, a AtDIR1 homologue from maize, after the infection of U. maydis. In addition to the genes up-regulated after U. maydis infection, there are also some ZmLTP genes show a decline of expression levels as time goes by, such as ZmLTPg23 and ZmLTPd6. These genes might function in other biological processes or respond to other kinds of pathogen attacks. Taken together, these results presented here indicated that many genes in this family might participate in the pathogen response.
Transcriptional responses of maize nsLTP genes against abiotic stresses
Apart from their inducibilities upon pathogen infection, nsLTP genes are also responsive to abiotic stresses like drought, cold and salt ,. To investigate the possible role of ZmLTPs under drought stress, the microarrary data of two inbred lines, drought-tolerant line Han21 and drought-sensitive line Ye478, under different drought treatments were used . Following whole-chip data processing, 57.14% (36 of 63) ZmLTP genes, about 41 probes (some genes have two or more probe sets) in the microarray, were extracted for analysis (Additional file 15: Table S11 and Figure 5A). Using fold change >2 and P value <0.05, ZmLTP genes with changeable expression profiles in at least one conditions were regarded as differentially expressed genes (Figure 5B). The result showed that 14 and 13 ZmLTPs were differentially regulated by drought and/or re-watering treatment(s) in Han21 and Ye478 maize inbred lines, respectively. Interestingly, some ZmLTPs appeared to be up-regulated in Han21 under drought stress, suggesting that they may have a role in drought stress response in maize. For example, the transcript levels of ZmLTP2.4, ZmLTP2.7 and ZmLTP2.8, were significantly up-regulated under moderate and severe drought stresses (Figure 5A). Furthermore, ZmLTP2.6 and ZmLTPd8 showed high-abundance transcript levels under moderate drought condition compared to untreated samples. To verify the expression patterns of these predicted genes under drought stress, 8 differentially expressed genes were selected for qRT-PCR analysis and the primers used in this assay are listed in Additional file 17: Table S12. As shown in Figure 5C, ZmLTP2.4, ZmLTP2.6, ZmLTP 2.8 and ZmLTPd8 were obviously up-regulated in water deficient, whereas ZmLTPd2, ZmLTPd9, ZmLTPg16 and ZmLTPg18 were down-regulated under drought stress. Taken together, qRT-PCR-based expression profiling for these selected genes confirmed the outcome of microarray analysis.
In addition to DNA microarray analysis, RNA-seq analysis was performed to investigate the functions of ZmLTPs involved in the differential drought sensitivity of maize reproductive (ovary) and vegetative tissues (leaf) under both drought and well-watered conditions . Nearly 66.67% (42) of ZmLTPs were reliably detected by RNA-seq data (Additional file 18: Table S13). Out of 42 expressed genes, 3 genes (ZmLTP2.8, ZmLTPg14 and ZmLTP2.2) were distinctly upregulated and 17 genes were down-regulated in the drought-stressed ovary tissue, while 10 (ZmLTPd3, ZmLTP2.8, ZmLTPd4, ZmLTPd1, ZmLTP2.6, ZmLTPg2, ZmLTP1.1, ZmLTPg3, ZmLTPg1 and ZmLTPd5) were significantly up-regulated and 2 (ZmLTP2.9 and ZmLTPg24) were down-regulated in the drought-stressed leaf meristem, suggesting that these ZmLTP genes may exhibit tissue-specific expression under drought stress.
Many nsLTPs in plants have been reported to be closely related to drought stress. Plants have a cuticle on the surface of their leaves that provides a protective barrier against the environmental adversities, such as drought stress. In Arabidopsis, AtLTPg4, which have close phylogenetic relationship with ZmLTPg14 in maize, play a role in wax and/or cutin monomer transport . In our study, the expression level of ZmLTPg14 was increased in drought-stressed ovary tissue, implying a possible role in drought resistance. ZmLTP1.2, which was reported as the target of calmodulin, was significantly down-regulated in drought-stressed ovary tissue but up-regulated in drought-stressed leaf meristem. As expected, a stresses-response gene, ZmLTP1.1, also increases its expression level in the drought-stressed leaf meristem, whereas it was down-regulated in the drought-stressed ovary tissue. In any case, the analysis of the expression profiles based on both microarray and RNA-seq methods proved that ZmLTPs may play important roles in maize tolerance to drought stress. While further work is required for understanding how ZmLTPs are mechanistically linked to drought stress responses in maize.
Apart from understanding the effect of drought stress on the transcript levels of various maize nsLTP genes, it is also important to conceive about the expression level of these genes under the salt stress. Root is the first organ directly exposed to salt stress. To better understand the roles of ZmLTP genes in salt stress responses, RNA-seq analysis was performed to explore the expression profiles of ZmLTPs in three root types of maize [primary root (PR), seminal root (SR) and crown root (CR)] under salt stress. Nearly 60.32% (38 of 63) of ZmLTPs were detected by RNA-seq. As shown in Additional file 19: Figure S6, 3 ZmLTP genes (ZmLTP2.4, ZmLTPd3 and ZmLTPg17) and 4 genes (ZmLTP2.8, ZmLTPd3, ZmLTPd11 and ZmLTPg17) were found to be at least two-fold up-regulated in the PR and SR under salt stress, respectively. Under salt stress, 11 ZmLTP genes (ZmLTP1.1, ZmLTP1.5, ZmLTP2.4, ZmLTP2.8, ZmLTPd3, ZmLTPd9, ZmLTPd11, ZmLTPg2, ZmLTPg4, ZmLTPg5 and ZmLTPg10) significantly increased their expression levels in CR. It can be concluded that the gene expression profiles of ZmLTPs are more variable in the CR than those in PR and SR under salt stress. One ZmLTP gene, ZmLTPd3, was highly induced by salt stress in PR, SR and CR, suggesting that it may play an important role in salt stress response in maize. In addition, ZmLTP1.1 was specifically expressed in CR, which is accordance with previous report that ZmLTP1.1 may participate in salt resistance .
It is well-known that plants are constantly exposed to a variety of environmental stresses. Among them, the low temperature constitutes a key factor that can influence plant growth, development and crop productivity. The accumulation of soluble sugar and proline are important to protect plants against cold stress . In previous research, the expression of rice cold-inducible nsLTP gene in transgenic Phalaenopsis amabilis improved its adaptive responses to cold stress. The examination of transgenic plants showed the increased accumulation of several compatible solutes such as total soluble sugars and proline . To further dissect the biological function of ZmLTPs in cold tolerance, we investigate the possible roles of ZmLTPs under cold stress condition in two inbred maize lines, chilling-sensitive ETH-DL3 and chilling-tolerant ETH-DH7 . Following whole-chip data processing, 54% (34 of 63) ZmLTPs have corresponding probe sets in the microarray, and a total of 36 probe sets on the maize 58 k NSF Array (some genes have two probe sets) were extracted for analysis (Additional file 20: Table S14). Three ZmLTPs (ZmLTPd12, ZmLTPg9 and ZmLTP1.5) were defined as differentially expressed genes in response to cold stress (Additional file 21: Figure S7). Surprisingly, all of the three differentially expressed genes were up-regulated under severe cold stress. Among them, ZmLTPd12 and ZmLTPg9 were differentially expressed under the severe cold stress conditions in ETH-DH7, and ZmLTP1.5 and ZmLTPd12 was up-regulated under severe cold stress conditions in ETH-DL3.
Our results revealed effects of drought, salt and cold stress on transcript abundance of several members of maize nsLTP gene family, indicative of their important roles in abiotic stresses. Moreover, the promoters for most of ZmLTPs were found to contain one or more important cis-regulatory elements such as ABRE, MBS, TC-rich repeat and C-repeat/DRE (Additional file 13: Table S9), which might be involved in regulation of expression of ZmLTPs under these stress conditions. Thus, expression profiling against abiotic stresses provided novel insights into specific and/or overlapping expression patterns of nsLTP genes. These analyses will contribute to the studies for unraveling the molecular mechanisms of signal transduction, stress response, and developmental processes.
We have identified 63 nsLTP genes in maize genome, and the members in this family can be divided into five types (1, 2, C, D and G). Comprehensive study of ZmLTPs provided some important features of this gene family such as gene structure, duplicated event, and three-dimensional structure. Additionally, gene ontology analysis was performed to obtain clues about biological function of the ZmLTPs. The analyses of putative upstream regulatory elements further indicated that ZmLTPs may be involved in diverse biological processes and provided new insights into the shared and different transcriptional regulation machineries of ZmLTPs. The expression profiles of ZmLTPs across the different developmental stages showed that some members of this family exhibit tissue-specific expression, indicative of their important roles in performing diverse developmental and physiological functions during the maize life cycle. Further, we focused on the roles of maize nsLTPs in stress responses. The result indicated that some ZmLTPs exhibited a delayed expression pattern after the infection of Ustilago maydis and differentially expressed under drought, salt and cold stress, which may be used in further studies for improving stress tolerance in maize. In a word, our results are significant contributions to maize molecular breeding with enhanced quality traits.
Retrieval and identification of nsLTPs in maize and sorghum genomes
Our search strategy used successive iterations of hidden Markov model (HMM) and BLAST searches to identify small, secreted nsLTPs in maize. The complete genome sequence of maize was downloaded from the maize sequence database (release version 5b; http://www.maizesequence.org/index.html). Starting with protein sequences annotated with the Pfam domain PF00234 (plant lipid transfer/seed storage/trypsin-alpha amylase inhibitor), we scanned maize proteome sequences using the hmmsearch from the HMMER 3.0 package  with default parameters. To increase the probability of detecting putative nsLTPs in maize, BLASTP searches were further performed against maize proteome sequences using all known nsLTP amino acid sequences from Arabidopsis and rice as baits. The results were merged, and redundant sequences were deleted. Then, the deduced protein sequences of candidate nsLTPs were manually examined to harbor the 8CM and proteins lacking the essential cysteine residues were removed. Additionally, the proteins encoding for the putative proline-rich proteins were excluded from further analysis. The remaining candidate proteins were submitted to the InterProScan (http://www.ebi.ac.uk/InterProScan) to verify the presence of nsLTP domains. The protein sequences of the 2S-albumins (At2S1 to At2S4) , and alpha amylase inhibitor (RATI)  were BLAST-searched against the rest of candidate nsLTP proteins to exclude possible cereal storage proteins. In addition, sorghum genome sequences were downloaded from the Joint Genome Institute plant genomics database (release version 1.4; http://www.Phytozome.net) and the same approach was used to identify nsLTPs in sorghum.
Primary sequence analysis
All identified nsLTPs were analyzed for the presence of potential signal peptide cleavage sites using SignalP v4.0 program . The subcellular localization of nsLTP was predicted by LocTREE 3 (https://rostlab.org/services/loctree3/). Following signal peptide removal, molecular weight and isoelectric point were computed for mature proteins with ExPasy's pI/Mw tool (http://web.expasy.org/compute_pi/). The eight Cys residues of the 8CM domain were excluded during pI calculation, since Cys residues in a disulfide bond have no charge. Three prediction tools, PredGPI , big-PI Plant Predictor  and GPI-SOM  were used to determine the presence of GPI-anchor sites. Phosphorylation sites of the nsLTPs were predicted by P3DB (http://www.p3db.org/). A sequence identity matrix of the mature nsLTP sequences was computed using BioEdit v22.214.171.124 (http://www.mbio.ncsu.edu/bioedit/bioedit.html), which enable us to determine the gene subfamily assignment. The nomenclature was based on the guidelines proposed by Edstam et al. .
Sequence alignment and phylogenetic reconstruction
The nsLTP amino acid sequences of Arabidopsis (A. thaliana) and rice (O. sativa) were downloaded from the TAIR (http://www.Arabidopsis.org/) and TIGR (http://rice.plantbiology.msu.edu/) databases, respectively. Multiple alignments of the 8CM part of the mature proteins were carried out using the ClustalW implemented in MEGA 5 software. Manually refined alignments were used as input in ProtTest v2.4 , which was set to all candidate models, a BIONJ tree, and the slow optimization strategy. The best model (WAG + I + G) was selected by the Akaike information criterion implemented in ProtTest. Then, the amino acid sequences of 8CM of maize, rice and Arabidopsis, totally 274 sequences were loaded into PhyML v3.0 software . The confidence level of each node was estimated by the bootstrap procedure using 100 resampling repetitions of the data .
Gene structure and conserved intron splicing site analysis
To explore the diverse exon-intron organization of maize nsLTP genes, we compared the predicted coding sequence (CDS) of maize nsLTP genes with their corresponding genomic sequences using GSDS software (http://gsds.cbi.pku.edu.cn). For a better visualization and comparison, the 5' UTR sequences were removed beforehand. The positions of the 8CM were retrieved with PAL2NAL (http://www.bork.embl.de/) by converting a multiple sequence alignment of 8CM sequences and the corresponding gene sequences into a codon alignment, and then marked in red .
Genome distribution and gene duplication analysis
The genes were plotted separately onto the ten maize chromosomes according to their ascending order of physical position (bp), from the short arm telomere to the long arm telomere. As a gene family may be expanded through tandem and segmental duplication events, we intended to identify the mechanisms involved in the expansion of ZmLTP gene family. ZmLTP genes separated by no more than 10 intervening genes were considered as tandem duplication . The putative othologs from rice, sorghum and Arabidopsis were assigned to corresponding maize nsLTP proteins according to searches from the TIGR database (http://rice.plantbiology.msu.edu/annotation_pseudo_apk.shtml). Syntenic information from all examined genes was collected from the Plant Genome Duplication Database (PGDD, http://chibba.agtec.uga.edu/duplication). The ideograms were created using Circos software .
Estimation of synonymous and non-synonymous substitution rates
Synonymous and nonsynonymous substitution rates of duplicated genes were calculated by DnaSP v5.0 . The the ratio of nonsynonymous to synonymous nucleotide substitution rates between paralogs was analyzed to detect the mode of selection. Ka/Ks ratio >1, <1 and =1 represent positive selection, purifying selection and neutral selection, respectively. The dates of the duplication events were calculated by the equation T = Ks/2λ × 10-6 Mya, (λ =6.5 × 10-9) ,.
Homology modeling of maize nsLTPs
All the maize nsLTP proteins were searched against the Protein Data Bank (PDB) by BLASTP with the default parameters to identify the best template. The data was fed in SWISS-MODEL  for predicting the protein structure by homology modeling. To visualize and edit the PDB models, interactive molecular graphics program Chimera v1.8 was used . Ramachandran plots for the evaluation of the model were calculated using Chimera, and the amino acid environment was evaluated using ERRAT at the UCLA-DOE Structure Analysis and Verification Server (http://nihserver.mbi.ucla.edu/ERRAT/). A sequence alignment of various nsLTPs was made using the ClustalW implemented in MEGA 5 software, and the alignment was manually inspected and adjusted via Jalview (http://www.jalview.org/).
Gene ontology analysis
The GO enrichment analysis was performed using the web-based tool AgriGO (http://bioinfo.cau.edu.cn/agriGO/index.php). The singular enrichment analysis (SEA) method was used to identify enriched GO terms, and the significantly enriched functional clusters were selected with an adjusted q-value less than 0.05. The program provides the output defining three categories of GO classification, namely, biological processes, cellular components and molecular functions.
Promoter analysis of maize nsLTP genes
To investigate cis-elements in promoter sequences of maize nsLTP genes, the upstream sequences (~1000 bp) of each identified ZmLTPs were retrieved from B73 maize genome sequence using a Perl script. The PlantCARE website (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) was used to identify cis-elements in the promoters . To avoid biases in analysis, only cis-elements identified in at least ten different genes were considered further.
Microarray-based gene expression analysis
To gain an insight into the family-wide expression profile, microarray-based data analyses of both development and stress responses were carried out. In order to achieve the expression profilings of ZmLTP genes among different organs and development stages, transcriptome data of the genome-wide gene expression atlas of maize inbred line B73 were obtained from publicly available database PLEXdb with the accession number ZM37. Gene expression data of drought stress and pathogen infection were downloaded from GEO database with accession numbers GSE16567 and GSE10023, respectively. The microarray data of cold stress were downloaded from ArrayExpress database with accession number E-MTAB-1315. To identify probe sets for ZmLTP genes, probe match against the maize GeneChip assemblies were performed on the Affymetrix Netaffx website (http://www.affymetrix.com/estore/) and PLEXdb. The data were normalized using a robust multi-chip average (RMA) algorithm . Limma package was applied to data processing, and heatmaps representing log2-transformed probe intensities were generated with the gplots package . A hierarchical clustering algorithm was applied to determine similar patterns in expression profiles. Using fold change >2 and P value <0.05, ZmLTP genes having changed expression profiles in at least one of the three conditions examined were regarded as differentially expressed genes under drought stress, whereas the response must be at least 2.8-fold in cold stress.
Differential gene expression profiles based on RNA-seq
Previous maize drought transcriptome performed by RNA-seq analysis was used to investigate the expression patterns of drought-responsive genes in drought-treated and well-watered fertilized reproductive (ovaries) and vegetative tissue (basal leaf meristem) . The ovaries and the basal leaf meristem each have four libraries corresponding to well-watered and drought-stressed tissue, with one biological replicate respectively. The transcript abundance of each gene was estimated by fragments per kilobase of exon per million fragments mapped (FPKM). To achieve the expression patterns of ZmLTP genes in three root types (PR, SR and CR) under salt stress, the RNA-seq data were downloaded from the ArrayExpress database with accession number E-GEOD-53995. The transcript abundance of each gene was estimated by reads per kilobase per million mapped reads (RPKM).
Real-time PCR analysis
Maize seeds were incubated at 30°C for germination and sown in plastic pots containing sandy loam soil of 0.1 m3 each. Pots of each plant were divided into 3 groups with 2 replicates. Plant growth was carried out in a controlled condition of 14 h photoperiod, 28 × 2°C Day/night temperature and 55% air humidity. At 34 Days from the date of sowing (at this time, all plants had 3-4 fully developed leaves), one set of pots was kept as control (well-watered) and two other sets were used for drought by withholding the water for 6 Days. At the end of this period, the volumetric soil water content (measured at the 5°Cm depth from the top of soil level) was approximately 12.5% for the drought-stressed plants, compared with approximately 30% for the control plants. The leaves of the seedlings were harvested, immediately frozen in liquid nitrogen and stored at -80°C. Then, total RNA was extracted from leaves after subjecting to drought stress using Trizol reagent, followed by DNase I treatment to remove any genomic DNA contamination. About 2 μg of total RNA was reverse-transcribed using Takara RNA PCR Kit to generate cDNA. QRT-PCR was performed on iCycler iQ5 Multicolor real-time PCR detection system (Bio-Rad) by using Power SYBR Green PCR Master Mix (Applied Biosystems). Two biological replicates of each sample and three technical replicates of each biological replicate were used. The thermal cycling conditions were as follows: 50°C for 2 min, 95°C for 10 min, 40°Cycles of 95°C for 15's and 60°C for 1 min. Quantification of gene expression was performed by 2-ΔCt method. Maizeubiquilin-1 (UBQ1) gene was used as internal control for normalization.
Availability of supporting data
All the supporting data are included as additional files.
KW conceptualized, designed the project. XZ and KW carried out computational analysis, microarray and RNA-seq expression profiling. KW and XZ performed the experiments. XZ and KW wrote the manuscript. Both authors read and approved the final manuscript.
Kader JC: Lipid-Transfer Proteins in Plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 627-654. 10.1146/annurev.arplant.47.1.627.
Boutrot F, Chantret N, Gautier MF: 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 Genomics. 2008, 9: 86-10.1186/1471-2164-9-86.
Debono A, Yeats TH, Rose JK, Bird D, Jetter R, Kunst L, Samuels L: Arabidopsis LTPG is a glycosylphosphatidylinositol-anchored lipid transfer protein required for export of lipids to the plant surface. Plant Cell. 2009, 21: 1230-1238. 10.1105/tpc.108.064451.
Lascombe MB, Bakan B, Buhot N, Marion D, Blein JP, Larue V, Lamb C, Prange T: The structure of "defective in induced resistance" protein of Arabidopsis thaliana, DIR1, reveals a new type of lipid transfer protein. Protein Sci. 2008, 17: 1522-1530. 10.1110/ps.035972.108.
Douliez JP, Michon T, Elmorjani K, Marion D: 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. 10.1006/jcrs.2000.0315.
Edstam MM, Viitanen L, Salminen TA, Edqvist J: Evolutionary history of the non-specific lipid transfer proteins. Mol Plant. 2011, 4: 947-964. 10.1093/mp/ssr019.
Han GW, Lee JY, Song HK, Chang C, Min K, Moon J, Shin DH, Kopka ML, Sawaya MR, Yuan HS, Kim TD, Choe J, Lim D, Moon HJ, Suh SW: Structural basis of non-specific lipid binding in maize lipid-transfer protein complexes revealed by high-resolution X-ray crystallography. J Mol Biol. 2001, 308: 263-278. 10.1006/jmbi.2001.4559.
Hoh F, Pons JL, Gautier MF, 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. 10.1107/S0907444905000417.
Maldonado AM, Doerner P, Dixon RA, Lamb CJ, Cameron RK: A putative lipid transfer protein involved in systemic resistance signalling in Arabidopsis. Nature. 2002, 419: 399-403. 10.1038/nature00962.
Buhot N, Gomes E, Milat ML, Ponchet M, Marion D, Lequeu J, Delrot S, Coutos-Thevenot P, Blein JP: Modulation of the biological activity of a tobacco LTP1 by lipid complexation. Mol Biol Cell. 2004, 15: 5047-5052. 10.1091/mbc.E04-07-0575.
Nieuwland J, Feron R, Huisman BA, Fasolino A, Hilbers CW, Derksen J, Mariani C: Lipid transfer proteins enhance cell wall extension in tobacco. Plant Cell. 2005, 17: 2009-2019. 10.1105/tpc.105.032094.
Sterk P, Booij H, Schellekens GA, Van Kammen A, De Vries SC: Cell-specific expression of the carrot EP2 lipid transfer protein gene. Plant Cell. 1991, 3: 907-921. 10.1105/tpc.3.9.907.
Lee SB, Go YS, Bae HJ, Park JH, Cho SH, Cho HJ, Lee DS, Park OK, Hwang I, Suh MC: Disruption of glycosylphosphatidylinositol-anchored lipid transfer protein gene altered cuticular lipid composition, increased plastoglobules, and enhanced susceptibility to infection by the fungal pathogen Alternaria brassicicola. Plant Physiol. 2009, 150: 42-54. 10.1104/pp.109.137745.
Edstam MM, Blomqvist K, Eklof A, Wennergren U, Edqvist J: Coexpression patterns indicate that GPI-anchored non-specific lipid transfer proteins are involved in accumulation of cuticular wax, suberin and sporopollenin. Plant Mol Biol. 2013, 83: 625-649. 10.1007/s11103-013-0113-5.
Huang MD, Chen TL, Huang AH: Abundant type III lipid transfer proteins in Arabidopsis tapetum are secreted to the locule and become a constituent of the pollen exine. Plant Physiol. 2013, 163: 1218-1229. 10.1104/pp.113.225706.
Zhang D, Liang W, Yin C, Zong J, Gu F, Zhang D: OsC6, encoding a lipid transfer protein, is required for postmeiotic anther development in rice. Plant Physiol. 2010, 154: 149-162. 10.1104/pp.110.158865.
Van Loon LC, Van Strien EA: The families of pathogenesis-related proteins, their activities, and comparative analysis of PR-1 type proteins. Physiol Mol Plant Pathol. 1999, 55: 85-97. 10.1006/pmpp.1999.0213.
Buhot N, Douliez JP, Jacquemard A, Marion D, Tran V, Maume BF, Milat ML, Ponchet M, Mikes V, Kader JC, Blein JP: A lipid transfer protein binds to a receptor involved in the control of plant defence responses. FEBS Lett. 2001, 509: 27-30. 10.1016/S0014-5793(01)03116-7.
Blein JP, Coutos-Thevenot P, Marion D, Ponchet M: From elicitins to lipid-transfer proteins: a new insight in cell signalling involved in plant defence mechanisms. Trends Plant Sci. 2002, 7: 293-296. 10.1016/S1360-1385(02)02284-7.
Hincha DK, Neukamm B, Sror HA, Sieg F, Weckwarth W, Ruckels M, Lullien-Pellerin V, Schroder W, Schmitt JM: Cabbage cryoprotectin is a member of the nonspecific plant lipid transfer protein gene family. Plant Physiol. 2001, 125: 835-846. 10.1104/pp.125.2.835.
Jung HW, Kim W, Hwang BK: Three pathogen-inducible genes encoding lipid transfer protein from pepper are differentially activated by pathogens, abiotic, and environmental stresses. Plant Cell Environ. 2003, 26: 915-928. 10.1046/j.1365-3040.2003.01024.x.
Cameron KD, Teece MA, Smart LB: Increased accumulation of cuticular wax and expression of lipid transfer protein in response to periodic drying events in leaves of tree tobacco. Plant Physiol. 2006, 140: 176-183. 10.1104/pp.105.069724.
Sarowar S, Kim YJ, Kim KD, Hwang BK, Ok SH, Shin JS: Overexpression of lipid transfer protein (LTP) genes enhances resistance to plant pathogens and LTP functions in long-distance systemic signaling in tobacco. Plant Cell Rep. 2009, 28: 419-427. 10.1007/s00299-008-0653-3.
Guo L, Yang H, Zhang X, Yang S: Lipid transfer protein 3 as a target of MYB96 mediates freezing and drought stress in Arabidopsis. J Exp Bot. 2013, 64: 1755-1767. 10.1093/jxb/ert040.
Lauga B, Charbonnel-Campaa L, Combes D: Characterization of MZm3-3, a Zea mays tapetum-specific transcript. Plant Sci. 2000, 157: 65-75. 10.1016/S0168-9452(00)00267-3.
Li C, Xie W, Bai W, Li Z, Zhao Y, Liu H: Calmodulin binds to maize lipid transfer protein and modulates its lipids binding ability. FEBS J. 2008, 275: 5298-5308. 10.1111/j.1742-4658.2008.06660.x.
Zou HW, Tian XH, Ma GH, Li ZX: Isolation and Functional Analysis of ZmLTP3, a Homologue to Arabidopsis LTP3. Int J Mol Sci. 2013, 14: 5025-5035. 10.3390/ijms14035025.
Royo J, Gomez E, Sellam O, Gerentes D, Paul W, Hueros G: Two maize END-1 orthologs, BETL9 and BETL9like, are transcribed in a non-overlapping spatial pattern on the outer surface of the developing endosperm. Front Plant Sci. 2014, 5: 180-10.3389/fpls.2014.00180.
Lehmann MS, Pebay-Peyroula E, Cohen-Addad C, Odani S: Crystallographic data for soybean hydrophobic protein. J Mol Biol. 1989, 210: 235-236. 10.1016/0022-2836(89)90303-3.
Wei F, Coe E, Nelson W, Bharti AK, Engler F, Butler E, Kim H, Goicoechea JL, Chen M, Lee S, Fuks G, Sanchez-Villeda H, Schroeder S, Fang Z, McMullen M, Davis G, Bowers JE, Paterson AH, Schaeffer M, Gardiner J, Cone K, Messing J, Soderlund C, Wing RA: Physical and genetic structure of the maize genome reflects its complex evolutionary history. PLoS Genet. 2007, 3: e123-10.1371/journal.pgen.0030123.
Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, Liang C, Zhang J, Fulton L, Graves TA, Minx P, Reily AD, Courtney L, Kruchowski SS, Tomlinson C, Strong C, Delehaunty K, Fronick C, Courtney B, Rock SM, Belter E, Du F, Kim K, Abbott RM, Cotton M, Levy A, Marchetto P, Ochoa K, Jackson SM, Gillam B: The B73 maize genome: complexity, diversity, and dynamics. Science. 2009, 326: 1112-1115. 10.1126/science.1178534.
Wei K, Chen Y, Xie D: Genome-scale evolution and phylodynamics of H5N1 influenza virus in China during 1996-2012. Vet Microbiol. 2013, 167: 383-393. 10.1016/j.vetmic.2013.09.030.
Yang X, Tuskan GA, Cheng MZ: Divergence of the Dof gene families in poplar, Arabidopsis, and rice suggests multiple modes of gene evolution after duplication. Plant Physiol. 2006, 142: 820-830. 10.1104/pp.106.083642.
Salse J, Bolot S, Throude M, Jouffe V, Piegu B, Quraishi UM, Calcagno T, Cooke R, Delseny M, Feuillet C: Identification and characterization of shared duplications between rice and wheat provide new insight into grass genome evolution. Plant Cell. 2008, 20: 11-24. 10.1105/tpc.107.056309.
Swigonova Z, Lai J, Ma J, Ramakrishna W, Llaca V, Bennetzen JL, Messing J: Close split of sorghum and maize genome progenitors. Genome Res. 2004, 14: 1916-1923. 10.1101/gr.2332504.
Williamson MP: The structure and function of proline-rich regions in proteins. Biochem J. 1994, 297 (Pt 2): 249-260.
Seo PJ, Lee SB, Suh MC, Park MJ, Go YS, Park CM: The MYB96 transcription factor regulates cuticular wax biosynthesis under drought conditions in Arabidopsis. Plant Cell. 2011, 23: 1138-1152. 10.1105/tpc.111.083485.
Wei KF, Wang YM, Xie DX: Identification and expression profile analysis of the protein kinase gene superfamily in maize development. Mol Breeding. 2014, 33: 155-172. 10.1007/s11032-013-9941-x.
Li R, Xia J, Xu Y, Zhao X, Liu YG, Chen Y: Characterization and genetic mapping of a Photoperiod-sensitive dwarf 1 locus in rice (Oryza sativa L.). Theor Appl Genet. 2014, 127: 241-250. 10.1007/s00122-013-2213-7.
Serna A, Maitz M, O°Connell T, Santandrea G, Thevissen K, Tienens K, Hueros G, Faleri C, Cai G, Lottspeich F, Thompson RD: Maize endosperm secretes a novel antifungal protein into adjacent maternal tissue. Plant J. 2001, 25: 687-698. 10.1046/j.1365-313x.2001.01004.x.
Gomez E, Royo J, Muniz LM, Sellam O, Paul W, Gerentes D, Barrero C, Lopez M, Perez P, Hueros G: The maize transcription factor myb-related protein-1 is a key regulator of the differentiation of transfer cells. Plant Cell. 2009, 21: 2022-2035. 10.1105/tpc.108.065409.
Huang MD, Wei FJ, Wu CC, Hsing YI, Huang AH: Analyses of advanced rice anther transcriptomes reveal global tapetum secretory functions and potential proteins for lipid exine formation. Plant Physiol. 2009, 149: 694-707. 10.1104/pp.108.131128.
Jang CS, Lee HJ, Chang SJ, Seo YW: Expression and promoter analysis of the TaLTP1 gene induced by drought and salt stress in wheat (Triticum aestivum L.). Plant Science. 2004, 167: 995-1001. 10.1016/j.plantsci.2004.05.019.
Zheng J, Fu J, Gou M, Huai J, Liu Y, Jian M, Huang Q, Guo X, Dong Z, Wang H, Wang G: Genome-wide transcriptome analysis of two maize inbred lines under drought stress. Plant Mol Biol. 2010, 72: 407-421. 10.1007/s11103-009-9579-6.
Kakumanu A, Ambavaram MM, Klumas C, Krishnan A, Batlang U, Myers E, Grene R, Pereira A: Effects of drought on gene expression in maize reproductive and leaf meristem tissue revealed by RNA-Seq. Plant Physiol. 2012, 160: 846-867. 10.1104/pp.112.200444.
Kim H, Lee SB, Kim HJ, Min MK, Hwang I, Suh MC: Characterization of Glycosylphosphatidylinositol-Anchored Lipid Transfer Protein 2 (LTPG2) and Overlapping Function between LTPG/LTPG1 and LTPG2 in Cuticular Wax Export or Accumulation in Arabidopsis thaliana. Plant Cell Physiol. 2012, 53: 1391-1403. 10.1093/pcp/pcs083.
Wanner LA, Junttila O: Cold-induced freezing tolerance in Arabidopsis. Plant Physiol. 1999, 120: 391-400. 10.1104/pp.120.2.391.
Qin X, Liu Y, Mao S, Li T, Wu H, Chu C, Wang Y: Genetic transformation of lipid transfer protein encoding gene in Phalaenopsis amabilis to enhance cold resistance. Euphytica. 2011, 177: 33-43. 10.1007/s10681-010-0246-4.
Sobkowiak A, Jonczyk M, Jarochowska E, Biecek P, Trzcinska-Danielewicz J, Leipner J, Fronk J, Sowinski P: Genome-wide transcriptomic analysis of response to low temperature reveals candidate genes determining divergent cold-sensitivity of maize inbred lines. Plant Mol Biol. 2014, 85: 317-331. 10.1007/s11103-014-0187-8.
Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011, 39: W29-W37. 10.1093/nar/gkr367.
Guerche P, Tire C, De Sa FG, De Clercq A, Van Montagu M, Krebbers E: Differential Expression of the Arabidopsis 2S Albumin Genes and the Effect of Increasing Gene Family Size. Plant Cell. 1990, 2: 469-478. 10.1105/tpc.2.5.469.
Strobl S, Maskos K, Wiegand G, Huber R, Gomis-Ruth FX, Glockshuber R: A novel strategy for inhibition of alpha-amylases: yellow meal worm alpha-amylase in complex with the Ragi bifunctional inhibitor at 2.5 A resolution. Structure. 1998, 6: 911-921. 10.1016/S0969-2126(98)00092-6.
Petersen TN, Brunak S, Von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8: 785-786. 10.1038/nmeth.1701.
Pierleoni A, Martelli PL, Casadio R: PredGPI: a GPI-anchor predictor. BMC Bioinformatics. 2008, 9: 392-10.1186/1471-2105-9-392.
Eisenhaber B, Wildpaner M, Schultz CJ, Borner GH, Dupree P, Eisenhaber F: Glycosylphosphatidylinositol lipid anchoring of plant proteins. Sensitive prediction from sequence- and genome-wide studies for Arabidopsis and rice. Plant Physiol. 2003, 133: 1691-1701. 10.1104/pp.103.023580.
Fankhauser N, Maser P: Identification of GPI anchor attachment signals by a Kohonen self-organizing map. Bioinformatics. 2005, 21: 1846-1852. 10.1093/bioinformatics/bti299.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
Wei K, Chen J, Wang Y, Chen Y, Chen S, Lin Y, Pan S, Zhong X, Xie D: Genome-wide analysis of bZIP-encoding genes in maize. DNA Res. 2012, 19: 463-476. 10.1093/dnares/dss026.
Wei K, Chen J, Chen Y, Wu L, Xie D: Multiple-strategy analyses of ZmWRKY subgroups and functional exploration of ZmWRKY genes in pathogen responses. Mol Biosyst. 2012, 8: 1940-1949. 10.1039/c2mb05483c.
Wei K, Pan S: Maize protein phosphatase gene family: identification and molecular characterization. BMC Genomics. 2014, 15: 773-10.1186/1471-2164-15-773. doi:10.1186/1471-2164-15-773
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res. 2009, 19: 1639-1645. 10.1101/gr.092759.109.
Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.
Lynch M, Conery JS: The evolutionary demography of duplicate genes. J Struct Funct Genomics. 2003, 3: 35-44. 10.1023/A:1022696612931.
Yang Z, Gu S, Wang X, Li W, Tang Z, Xu C: Molecular evolution of the CPP-like gene family in plants: insights from comparative genomics of Arabidopsis and rice. J Mol Evol. 2008, 67: 266-277. 10.1007/s00239-008-9143-z.
Arnold K, Bordoli L, Kopp J, Schwede T: The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006, 22: 195-201. 10.1093/bioinformatics/bti770.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE: UCSF Chimera–a visualization system for exploratory research and analysis. J Comput Chem. 2004, 25: 1605-1612. 10.1002/jcc.20084.
Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouze P, Rombauts S: PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30: 325-327. 10.1093/nar/30.1.325.
Wei KF, Chen J, Chen YF, Wu LJ, Xie DX: Molecular phylogenetic and expression analysis of the complete WRKY transcription factor family in maize. DNA Res. 2012, 19: 153-164. 10.1093/dnares/dsr048.
Wei K, Chen Y, Lin Y, Pan Y: Genetic dynamic analysis of the influenza A H5N1 NS1 gene in China. PLoS One. 2014, 9: e101384-10.1371/journal.pone.0101384.
We are grateful to the providers who submitted the microarray and RNA-seq data to the public expression databases which can be applied freely. We sincerely thank Dr. Yanmei Wang in Tsinghua University for revising our manuscript. We also thank Yanhui Chen, Yina Lin, Si Pan, and Xiaoyao Liu in our laboratory for useful discussions.
The authors declare that they have no competing interests.
Electronic supplementary material
Additional file 9: Figure S3.: The map of intron-exon arrangement of maize nsLTP genes. Introns and exons are drawn to scale with the full encoding regions of their respective. Exons are depicted as green boxes, introns as connecting thin lines and 8CM as red boxes. Non-translated regions, when supported by full-length cDNA sequences, are shown in blue boxes. (PDF 41 KB)
Additional file 10: Figure S4.: The map of intron-exon arrangement of Arabidopsis and rice nsLTP genes. Introns and exons are drawn to scale with the full encoding regions of their respective. Exons are depicted as green boxes, introns as connecting thin lines and 8CM as red boxes. Non-translated regions, when supported by full-length cDNA sequences, are shown in blue boxes. (PDF 69 KB)
Additional file 16: Figure S5.: Normalized signal intensities of 41 probe sets representing 36 ZmLTP genes displayed for the microarray experiment, which was treated with Ustilago maydis. Log2 ratio values (shown by a green-red gradient) represents the fold change of the gene expression in the fungal infections. (PDF 24 KB)
About this article
Cite this article
Wei, K., Zhong, X. Non-specific lipid transfer proteins in maize. BMC Plant Biol 14, 281 (2014). https://doi.org/10.1186/s12870-014-0281-8
- Gene Ontology
- Salt Stress
- Drought Stress
- Cold Stress
- Seminal Root