- Research article
- Open Access
Non-specific lipid transfer proteins in maize
© Wei and Zhong; licensee BioMed Central Ltd. 2014
- Received: 3 July 2014
- Accepted: 9 October 2014
- Published: 28 October 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.
- Gene Ontology
- Salt Stress
- Drought Stress
- Cold Stress
- Seminal Root
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.
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
Some characteristics for the different types of non-specific lipid transfer proteins found in maize
Number of members
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
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
Ka/Ks analysis and estimate of the absolute dates for the duplication events between the duplicated ZmLTP genes
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.
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
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
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
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 v188.8.131.52 (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.
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.
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Williamson MP: The structure and function of proline-rich regions in proteins. Biochem J. 1994, 297 (Pt 2): 249-260.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Wanner LA, Junttila O: Cold-induced freezing tolerance in Arabidopsis. Plant Physiol. 1999, 120: 391-400. 10.1104/pp.120.2.391.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011, 39: W29-W37. 10.1093/nar/gkr367.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Pierleoni A, Martelli PL, Casadio R: PredGPI: a GPI-anchor predictor. BMC Bioinformatics. 2008, 9: 392-10.1186/1471-2105-9-392.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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-773PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.View ArticlePubMedGoogle Scholar
- Lynch M, Conery JS: The evolutionary demography of duplicate genes. J Struct Funct Genomics. 2003, 3: 35-44. 10.1023/A:1022696612931.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.