- Open Access
START lipid/sterol-binding domains are amplified in plants and are predominantly associated with homeodomain transcription factors
Genome Biology volume 5, Article number: R41 (2004)
In animals, steroid hormones regulate gene expression by binding to nuclear receptors. Plants lack genes for nuclear receptors, yet genetic evidence from Arabidopsis suggests developmental roles for lipids/sterols analogous to those in animals. In contrast to nuclear receptors, the lipid/sterol-binding StAR-related lipid transfer (START) protein domains are conserved, making them candidates for involvement in both animal and plant lipid/sterol signal transduction.
We surveyed putative START domains from the genomes of Arabidopsis, rice, animals, protists and bacteria. START domains are more common in plants than in animals and in plants are primarily found within homeodomain (HD) transcription factors. The largest subfamily of HD-START proteins is characterized by an HD amino-terminal to a plant-specific leucine zipper with an internal loop, whereas in a smaller subfamily the HD precedes a classic leucine zipper. The START domains in plant HD-START proteins are not closely related to those of animals, implying collateral evolution to accommodate organism-specific lipids/sterols. Using crystal structures of mammalian START proteins, we show structural conservation of the mammalian phosphatidylcholine transfer protein (PCTP) START domain in plants, consistent with a common role in lipid transport and metabolism. We also describe putative START-domain proteins from bacteria and unicellular protists.
The majority of START domains in plants belong to a novel class of putative lipid/sterol-binding transcription factors, the HD-START family, which is conserved across the plant kingdom. HD-START proteins are confined to plants, suggesting a mechanism by which lipid/sterol ligands can directly modulate transcription in plants.
The StAR-related lipid transfer (START) domain, named after the mammalian 30 kDa steroidogenic acute regulatory (StAR) protein that binds and transfers cholesterol to the inner mitochondrial membrane , is defined as a motif of around 200 amino acids implicated in lipid/sterol binding . Ligands have been demonstrated for a small number of START-domain proteins from animals. The mammalian StAR and metastatic lymph node 64 (MLN64) proteins both bind cholesterol , the phosphatidylcholine transfer protein (PCTP) binds phosphatidylcholine , and the carotenoid-binding protein (CBP1) from silkworm binds the carotenoid lutein . In addition, a splicing variant of the human Goodpasture antigen-binding protein (GPBP) called CERT was recently shown to transport ceramide via its START domain .
The structure of the START domain has been solved by X-ray crystallography for three mammalian proteins: PCTP , MLN64  and StarD4 . On the basis of the structural data, START is classified as a member of the helix-grip fold superfamily, also termed Birch Pollen Allergen v1 (Bet v1)-like, which is ubiquitous among cellular organisms . Iyer et al.  used the term 'START superfamily' as synonymous with the helix-grip fold superfamily. Here we use the nomenclature established in the Protein Data Bank (PDB)  and Structural Classification of Proteins (SCOP)  databases, restricting the use of the acronym 'START' to members of the family that are distinguished by significant amino-acid sequence similarity to the mammalian cholesterol-binding StAR protein. Members of the START family are predicted to bind lipids or sterols [2, 11], whereas other members of the helix-grip fold superfamily are implicated in interactions with a wide variety of metabolites and other molecules such as polyketide antibiotics, RNA or antigens .
The presence of START domains in evolutionarily distant species such as animals and plants suggests a conserved mechanism for interaction of proteins with lipids/sterols . In mammalian proteins such as StAR or PCTP, the START domain functions in transport and metabolism of a sterol or phospholipid, respectively. START domains are also found in various multidomain proteins implicated in signal transduction , suggesting a regulatory role for START-domain proteins involving lipid/sterol binding.
To investigate the evolutionary distribution of the START domains in plants in comparison to other cellular organisms and to study their association with other functional domains, we applied a BLASTP search to identify putative START-containing protein sequences (see Materials and methods). We focused our study on proteins from the sequenced genomes of Arabidopsis thaliana (Table 1), rice (Table 2), humans, Drosophila melanogaster and Caenorhabditis elegans, as well as Dictyostelium discoideum (Table 3), in addition to sequences from bacteria and unicellular protists (Table 4). CBP1 from the silkworm Bombyx mori was also included in our analysis (Table 3). Figure 1 presents a phylogenetic tree comparing the START domains from the plant Arabidopsis to those from the animal, bacterial and protist kingdoms.
Results and discussion
Evolution of START domains in multicellular organisms
Our findings show that START domain-containing proteins are amplified in plant genomes (Arabidopsis and rice) relative to animal genomes (Figures 1,2). Arabidopsis and rice contain 35 and 29 START proteins each, whereas the human and mouse genomes contain 15 each , and C. elegans and D. melanogaster encode seven and four, respectively. In comparison, bacterial and protist genomes appear to encode a maximum of two START proteins (see below).
START-domain minimal proteins comprising the START domain only, as well as START proteins containing additional sequence of unknown or known function appear to be conserved across plants, animals, bacteria and protists (Tables 1,2,3,4, Figure 1). However, only in plants, animals and multicellular protists (D. discoideum) are START domains found in association with domains having established functions in signal transduction or transcriptional control, consistent with the idea that START evolved as a regulatory domain in multicellular eukaryotes. The cellular slime mold D. discoideum, which progresses from unicellular to multicellular developmental stages, contains an unusual START-domain protein  which has so far not been found in any other organism: FbxA/CheaterA (ChtA), an F-Box/WD40 repeat-containing protein [12, 13]. FbxA/ChtA is thought to encode a component of an SCF E3 ubiquitin ligase implicated in cyclic AMP metabolism and histidine kinase signaling during development . Mutant analysis shows that FbxA/ChtA function is required to generate the multicellular differentiated stalk fate .
Functional domains that were found associated with START in animals include pleckstrin homology (PH), sterile alpha motif (SAM), Rho-type GTPase-activating protein (RhoGAP), and 4-hydroxybenzoate thioesterase (4HBT) (Table 3), consistent with a previous report . The RhoGAP-START configuration is absent from plants, but is conserved across the animal kingdom from mammals to insects and nematodes. The RhoGAP-START combination in addition to an amino-terminal SAM domain is apparent only in proteins from humans, mouse, and rat, indicating that SAM-RhoGAP-START proteins are specific to mammals. Similarly, the 4HBT-START combination, also referred to as the acyl-CoA thioesterase subfamily , is found exclusively in proteins from humans, mouse and rat, and therefore seems to have evolved in the mammalian lineage.
In humans, about half of the START domain-containing proteins (6/15) are multidomain proteins, whereas in Arabidopsis and rice approximately three-quarters (26/35; 22/29) of START proteins contain an additional domain. The largest proportion of Arabidopsis and rice multidomain START proteins (21/26; 17/22) contain a homeodomain (HD), while a smaller group of proteins (4/26; 4/22) contain a PH domain together with a recently identified domain of unknown function 1336 (DUF1336) motif. In addition, a single START-DUF1336 protein of about the same size, but lacking strong sequence similarity to PH at its amino terminus, is present in both Arabidopsis and rice. It is striking that the sequence of the START domain correlates with the type of START protein, an indication that evolutionary speciation through duplication and subsequent sequence evolution of START domains took place after initial manifestation of novel protein architecture by domain shuffling.
The position of the START domain in proteins larger than 300 amino acids varies between plant, animal and protist kingdoms. For example, in human proteins, START is always near the carboxy terminus (1-55 amino acids from the end) (Table 3). In plant proteins, however, the START domain is not strictly confined to the carboxy terminus. In both Arabidopsis and rice the START domain can be positioned as much as approximately 470 amino acids from the carboxy terminus (HD-ZIP START proteins: Tables 1,2). Moreover, in a subset of plant proteins (PH-START DUF1336 proteins), the START domain is positioned centrally between two different domains. However, defined functional domains are typically amino terminal of the START domain in both animals and plants. By contrast, in the sole example of a START-domain protein in D. discoideum, FbxA/ChtA, the START domain is present at the amino terminus, with F-Box and WD40 domains positioned after it.
HD-START transcription factors are unique to plants
The START-domain proteins from Arabidopsis were classified into seven subfamilies according to their structures and sizes (Figure 2a). The majority of START domains are found in transcription factors of the HD family. HDs are DNA-binding motifs involved in the transcriptional regulation of key developmental processes in eukaryotes. However, only within the plant kingdom do HD transcription factors also contain START domains (Figure 1). Among around 90 HD family members in Arabidopsis , approximately one-quarter (21) contain a START domain. All HD-START proteins contain a putative leucine zipper, a dimerization motif that is not found in HD proteins from animals or yeast. Nuclear localization has been demonstrated for two HD-START proteins: GLABRA2 (GL2)  and REVOLUTA/INTERFASCICULAR FIBERLESS1 (REV/IFL1) . Furthermore, canonical DNA-binding sites are reported for GL2  and two other HD-START transcription factors, A. thaliana MERISTEM LAYER1 (ATML1) , and PROTODERMAL FACTOR2 (PDF2) .
A similar spectrum of START domain-containing proteins is found in Arabidopsis and rice, suggesting their origin in a common ancestor (Figure 2b). The size of the rice genome (430 Mb) is roughly four times that of Arabidopsis (120 Mb). Despite a twofold difference between the total number of predicted genes in rice (ssp. indica: 53,398 ) versus Arabidopsis (~28,000), the number of START domains per genome appears to be relatively constant: Arabidopsis and rice contain 35 and 29 START genes, respectively. Thus, START-domain genes belong to the subset of Arabidopsis genes (estimated at two-thirds) that are present in rice . However, one intriguing exception is the apparent absence of rice proteins orthologous to two unusual Arabidopsis START proteins (At4g26920 and At5g07260), which share sequence similarity to each other and to members of the HD-ZLZ START subfamily, but lack HD and zipper-loop-zipper (ZLZ) domains (Figure 2b; Tables 1,2). Their absence from rice makes them candidate dicot-specific START proteins.
Screening for expressed sequence tags (ESTs) by BLASTN was conducted to determine whether the types of START sequences from Arabidopsis and rice are also present in other plants (see Materials and methods). The screen detected 185 START domain-encoding sequences from a wide assortment of plants representing 25 different species. Consistent with our findings in Arabidopsis and rice (Tables 1,2), START domains were found in the plant-specific combinations (HD-START and PH-START) in both dicot and monocot members of the angiosperm division. ESTs for HD-START transcription factors were also identified from the gymnosperm Picea abies (AF328842 and AF172931), as well as from a representative of the most primitive extant seed plant, the cycad Cycas rumphii (CB093462). Furthermore, a HD-START sequence is expressed in the moss Physcomitrella patens (AB032182). Thus it appears that the HD-START plant-specific configuration evolved in the earliest plant ancestor, or alternatively has been retained in the complete plant lineage.
Two different HD-associated leucine zippers are found in HD-START proteins
Sequence alignments and phylogenetic analysis revealed two distinct classes of HD-START proteins, which differ substantially in their leucine zippers and START domains (Figures 1,2,3). Both types of leucine zipper are unrelated in sequence to the homeobox-associated leucine zipper (HalZ), which is a plant-specific leucine zipper found in other HD proteins lacking START .
Most HD-START proteins (16/21 in Arabidopsis; 12/17 in rice) contain a leucine zipper with an internal loop (defined here as zipper-loop-zipper, ZLZ; also termed 'truncated leucine zipper motif' ) immediately following a conserved HD domain (Figure 3a). The ZLZ motif appears to be less conserved than the classic basic region leucine zipper and seems to be plant specific. It was shown to be functionally equivalent to the HalZ leucine zipper domain for dimerization in an in vitro DNA binding assay .
The other HD-START proteins (5/21 in Arabidopsis; 5/17 in rice) contain a classic leucine zipper DNA-binding motif fused to the end of the HD, designated here as ZIP (Figure 3b). This leucine zipper shows strong sequence similarity to the basic region leucine zipper domains (bZIP and BRLZ) [25, 26], which have overlapping consensus sequences and are found in all eukaryotic organisms.
Despite these differences, it is likely that both types of HD-START transcription factors originated from a common ancestral gene. They share a common structural organization in their amino-terminal HD, leucine zipper (ZLZ or ZIP) and START domains (Figure 2a). Moreover, the carboxy terminus of HD-ZLZ START proteins (approximately 250 amino acids) shares sequence similarity with the first 250 amino acids of the approximately 470 amino acids at the carboxy terminus of HD-ZIP START proteins. This is exemplified by a comparison between the carboxy-terminal sequences of ATML1 (HD-ZLZ START) and REV (HD-ZIP START), which are 20% identical and 39% similar.
HD-START proteins are implicated in cell differentiation during plant development
Several HD-ZLZ START genes correspond to striking mutant phenotypes in Arabidopsis, and for numerous HD-ZLZ START genes, functions in the development of the epidermis have been implicated. Proteins of the HD-ZLZ START subfamily share strong sequence similarity to each other along their entire lengths, including the carboxy-terminal sequence (approximately 250 amino acids) of unknown function that follows the START domain. The HD-ZLZ transcription factors ATML1 and PDF2 appear to be functionally redundant: double-mutant analysis shows that the corresponding genes are required for epidermal differentiation during embryogenesis . The rice HD-ZLZ START protein RICE OUTERMOST CELL-SPECIFIC GENE1 (ROC1) seems to have an analogous function to ATML1 in that its expression is restricted to the outermost epidermal layer from the earliest stages in embryogenesis . Another HD-ZLZ gene from rice, Oryza sativa TRANSCRIPTION FACTOR 1 (OSTF1), appears to be developmentally regulated during early embryogenesis and is also expressed preferentially in the epidermis . Mutations in Arabidopsis ANTHOCYANINLESS2 (ANL2) affect anthocyanin accumulation and the cellular organization in the root, indicating a role in subepidermal cell identity . The GL2 gene is expressed in specialized epidermal cells and mutant analysis reveals its function in trichome and non-root hair cell fate determination [24, 29]. GL2 functions as a negative regulator of the phosopholipid signaling in the root , raising the possibility that the activity of GL2 itself is regulated through a feedback mechanism of phospholipid signaling through its START domain.
The HD-ZIP START genes characterized thus far are implicated in differentiation of the vasculature. Members of this subfamily are typically large proteins (more than 830 amino acids) that display strong sequence similarity to each other along their entire lengths, including the carboxy-terminal 470 or so residues of unknown function that follow the START domain. Mutations affecting PHABULOSA (PHB) and PHAVOLUTA (PHV), which have redundant functions, abolish radial patterning from the vasculature in the developing shoot, and perturb adaxial/abaxial (upper/lower) axis formation in the leaf . Mutant analysis reveals that REV [31, 32], isolated independently as IFL1 [17, 33], is also involved in vascular differentiation. Although a mutant phenotype for A. thaliana HOMEOBOX-8 (ATHB-8) is not reported, its expression is restricted to provascular cells  and promotes differentiation in vascular meristems .
The presence of the START domain in HD transcription factors suggests the possibility of lipid/sterol regulation of gene transcription for HD-START proteins, as previously hypothesized . One advantage of such a mechanism is that the metabolic state of the cell in terms of lipid/sterol synthesis could be linked to developmental events such as regulation of transcription during differentiation. Changes in the activity of a HD-START transcription factor could be controlled via a lipid/sterol-binding induced conformational change. For instance, a protein-lipid/sterol interaction involving the START domain may regulate the activity of the transcription factor directly by affecting its DNA-binding affinity or interaction with accessory proteins at the promoter. Alternatively, or in addition, protein-lipid/sterol binding may positively or negatively affect transport or sequestration of the transcription factor to the nucleus.
PH-START proteins differ in plants and animals
A subset of animal and plant START proteins contain an amino-terminal PH domain, which is found in a wide variety of eukaryotic proteins implicated in signaling. PH domains are characterized by their ability to bind phosphoinositides, thereby influencing membrane and/or protein interactions . In some cases, phosphoinositide interactions alone may not be sufficient for membrane association, but may require cooperation with other cis-acting anchoring motifs, such as the START domain, to drive membrane attachment.
Although both plant and animal genomes encode START domains in association with an amino-terminal PH domain, the sequences of the PH-START proteins are not conserved between kingdoms (Figure 1; data not shown). In plants, the START domain is adjacent to the PH domain, whereas in animals the PH and START domains are separated by two serine-rich domains . The PH-START protein from humans, GPBP, has serine/threonine kinase activity and Goodpasture (GP) antigen binding affinity, two functions that involve the serine-rich domains. In contrast, the plant PH-START proteins contain a plant-specific carboxy-terminal domain (of around 230 amino acids) of unknown function, DUF1336 (Protein families database (Pfam)) . In addition, amino-terminal sequence analysis (TargetP; see Materials and methods) predicts that three PH-START proteins from Arabidopsis (At3g54800, At4g19040 and At5g45560) and two PH-START proteins from rice (BAD07818 and BAC22213) localize to mitochondria. This suggests a common lipid/sterol-regulated function of these proteins that is related to their subcellular localization.
Membrane localization of START-domain proteins
Transmembrane segments may act to tether START-domain proteins to intracellular membranes. One START protein, MLN64 from humans, is anchored to late endosomes via four putative transmembrane segments that are required for its localization [38, 39]. The carboxy-terminal START domain of MLN64 faces the outside of this compartment, consistent with a role in cholesterol trafficking to a cytoplasmic acceptor . Using sequence scanning and alignments with MLN64, putative transmembrane segments in the MLN64-related proteins from D. melanogaster (CG3522/Start1)  and C. elegans (3F991) are predicted (Table 3). No other START proteins from animals are predicted to contain transmembrane segments. A single putative amino-terminal transmembrane segment is predicted for one bacterial START protein (Xanthomonas axonopodis XAC0537), while one START protein from the unicellular protist Cryptosporidium parvum (1Mx.08) contains seven putative transmembrane segments amino-terminal to the START domain (Table 4). Transmembrane segments are also predicted for several Arabidopsis and rice START proteins (Tables 1,2). Among Arabidopsis START-domain proteins, four (At1g55960, At1g64720, At3g13062, and At3g23080) are predicted to encode a single transmembrane segment at the amino-terminal end of the protein. Two START proteins of very similar sequence (At4g14500 and At5g54170) contain two putative transmembrane segments, one at the amino terminus of the protein and the second very close to the carboxy terminus. This feature is also found in one rice protein (ID215312), and appears to be plant specific. Possible functions of these putative transmembrane proteins include lipid/sterol transport to and from membrane-bound organelles and/or maintenance of membrane integrity.
Putative ligands for START domains in plants
Cholesterol, phosphatidylcholine, carotenoid and ceramide are examples of lipids that are known to bind START domains from animals [3–6]. We explored the potential for predicting lipid/sterol ligands for plant-specific START domains by aligning plant START domains with related animal START domains having defined ligands. However, the mammalian MLN64 and StAR START domains, which have been shown to bind cholesterol, do not share convincing sequence conservation with any particular subfamily of Arabidopsis or rice START domains. Similarly, the START domains of CBP1, which binds the carotenoid lutein, and CERT, which binds ceramide, both show insufficient amino-acid conservation when compared in alignment to the most closely related Arabidopsis proteins. Knowledge of the precise ligand contacts within the cavity of these START domains is required to make more accurate predictions about their specific binding properties.
Although most START domains appear to have diverged in the evolution of plants and animals, one notable exception is the preservation of PCTP-like sequences (Figures 1,4a). The major function of mammalian PCTP, a START-domain minimal protein, is to replenish the plasma membrane with phosphatidylcholine and to allow its efflux. Mammalian PCTP binds phosphatidylcholine and its crystal structure in the ligand-bound form has recently been solved . Phosphatidylcholine is the most abundant phospholipid in animals and plants. It is a key building block of membrane bilayers, comprising a high proportion of the outer leaflet of the plasma membrane.
To identify plant START domains that are most similar to mammalian PCTP, we constructed a phylogeny with PCTP against the complete set of 64 START domains from Arabidopsis and rice. Figure 4a shows the branch from this phylogeny grouping mammalian PCTP with START domains from 10 plant proteins. Unlike PCTP, a protein of 214 amino acids, the plant PCTP-like proteins are larger (around 400 amino acids) and contain divergent amino-terminal and carboxy-terminal sequences in addition to the START domain. Most (8/10) of these proteins are predicted to encode an amino-terminal transmembrane segment while a subset (3/8) also encode a putative carboxy-terminal transmembrane segment (Tables 1,2). A sequence alignment of PCTP against the START domains from these 10 PCTP-like proteins is illustrated in Figure 4b. Among 27 amino-acid residues in PCTP that contact the phosphatidylcholine ligand  (Figure 4b), 63% (17/27) are similar to residues in one or more of the plant proteins.
Homology modeling was used to generate three-dimensional structures of the 10 PCTP-like START domains from Arabidopsis and rice, incorporating templates from PCTP and two other mammalian START domain structures (MLN64  and StarD4 ). The models proved most accurate for the At3g13062 and ID215312 START domains from Arabidopsis and rice, respectively. Figure 4c-f depicts the model of the 212-amino-acid START domain from At3g13062. Overall, the PCTP and At3g13062 START domain sequences are 24% identical and 45% similar.
The backbone of the At3g13062 START domain overlaps with the crystal structure of PCTP START . The model exhibits α helices at the amino and carboxy termini separated by nine β strands and two shorter α helices (Figure 4c). The curved β sheets form a deep pocket with the carboxy-terminal α helix (α4) acting as a lid, resulting in an internal hydrophobic cavity (Figures 4e-g). The amino-terminal α helix (α1), which is in a region of the START domain that is least conserved, is peripheral to the hydrophobic cavity, and is not predicted to contact the ligand. Lipid entry or egress would require a major conformational change, most likely opening or unfolding of the carboxy-terminal α helix lid. In the case of PCTP, the α4 helix has been implicated in membrane binding and phosphatidylcholine extraction . The overall conformation is distinct from other hydrophobic cavity lipid-binding proteins, such as the intracellular sterol carrier protein 2 (SCP2)  from animals, and the orthologous nonspecific lipid-transfer protein from plants . It is not clear how specific START is for its particular ligand, as the START domain of human StAR has been shown in vitro to have affinity for both cholesterol and sitosterol , the latter being the major plant sterol. However, modeling of highly similar START domains based on the crystal structure of PCTP implicates a similar molecular ligand, and suggests that these proteins have retained a common function in evolution.
START domain-containing genes in unicellular organisms
Of about 100 bacterial genomes sequenced, only eight species encode proteins with similarity to START, whereas one-half (4/8) of unicellular protist genomes thus far sequenced appear to contain a START domain (Figure 1, Table 4). Expression data regarding the putative START-encoding genes from bacteria and unicellular protists is lacking to date. Most of the bacterial sequences represent START-domain minimal proteins of less than 245 amino acids. A START gene from the opportunistic pathogen Pseudomonas aeruginosa (PA1579) was previously hypothesized to have arisen by horizontal gene transfer [7, 8]. Here we show that similar START-domain genes also occur in one other pathogenic Pseudomonas species as well as in two related nonpathogens. The citrus-tree pathogen X. axonopodis and human pathogen Vibrio vulnificus also contain START genes. Nonpathogenic species that contain START genes include the anaerobic green sulfur bacterium Chlorobium tepidum and the anaerobic dehalogenating bacterium Desulfitobacterium hafniense. Among unicellular protists, the mammalian pathogens C. parvum, Giardia lamblia, Plasmodium falciparum, and Plasmodium yoelii yoelii contain START genes. By contrast, the genome of the unicellular green alga Chlamydomonas reinhardtii lacks START domains (Shinhan Shiu, personal communication).
START domains appear to be absent from the Archaea. Moreover the yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe, as well as another fungal species, Neurospora crassa, lack START domains. Perhaps START-domain genes are rarely found in prokaryotes and protists because their function is not required for normal unicellular proliferation. Their presence in a number of pathogenic species is consistent with the hypothesis of horizontal acquisition. This raises the possibility that in pathogens the START-domain proteins have a role in the interaction with host cells. For example, pathogens may recognize, transport and utilize lipids/sterols from the host organism via their START domains.
In this study we have explored putative START domains from plants in comparison to animal, protist and bacterial START domains. The function of START in binding lipid/sterol ligands has been demonstrated in only a few cases in the context of mammalian cells. One START domain configuration appears to be conserved from animals to plants: mammalian PCTP, which specifically binds phosphatidylcholine and exhibits sequence similarity to six Arabidopsis and four rice START domains. Using three-dimensional molecular modeling we predict that the plant proteins will bind orthologous lipids/sterols. Thus our finding provides a rationale for experimentally testing potential lipid/sterol interactions with PCTP-like START proteins. The majority of START domains appear to have evolved divergently in animals and plants. This is exemplified by RhoGAP START proteins, which are found in animals but are absent from plants. Strikingly, plant genomes contain two subfamilies of HD-START proteins that are not found in animals. Our results illustrate that these plant-specific START domains are amplified and conserved in dicot and monocot plant genomes, and that they originated in an ancient ancestor of the plant lineage. The modular coupling detected raises the possibility of regulation of signal transduction by lipids/sterols in plants. Although HD-START proteins are unique to plants, the presence of the lipid/sterol-binding domain in a transcription factor has a parallel in the animal kingdom (Figure 5). In animals, steroid-binding transcription factors of the nuclear receptor superfamily control numerous metabolic and developmental processes. The structure of a classic nuclear receptor consists of amino-terminal transactivation and zinc-finger DNA-binding domains and a carboxy-terminal steroid hormone receptor domain of around 240 amino acids . Steroid nuclear receptors are complexed with the proteins Hsp90, Hsp70, immunophilins or cyclophilins in the cytoplasm and upon steroid hormone binding move into the nucleus to bind specific DNA response elements for transcriptional control of target genes. The presence of HD-START transcription factors in plants leads us to postulate that these have an analogous role by binding regulatory lipids/sterols. The next hurdle in elucidating the role of the START domain in plant signal transduction will be the identification of in vivo binding partners.
Materials and methods
Searching databases for putative START proteins and transcribed sequences
Putative START proteins were identified by BLASTP against databases at the National Center for Biotechnology Information (NCBI) , The Institute for Genomic Research (TIGR) , and the Munich Information Center for Protein Sequences (MIPS) . Putative START proteins from rice Oryza sativa ssp. indica were identified by BLASTP of sequences predicted on whole-genome shotgun deposited in MIPS Oryza sativa database (MOsDB)  using START-domain proteins from Arabidopsis (GL2, ATML1), rice (ROC1, AAP54082) and humans (PCTP) as query. START-domain proteins were verified using InterPro  and NCBI Conserved Domain (CD)-Search [51, 52], using cutoff E-values set to 1e-05. Arabidopsis proteins were checked against Arabidopsis EST databases using WU-BLAST2  to identify corresponding transcribed sequences. Supporting mRNA expression data was obtained from the microarray elements search tool  available through The Arabidopsis Information Resource (TAIR). Rice ESTs and cDNAs were identified using BLASTN at E-values of < 1e-60 (with most sequences showing an E-value = 0) and other plant ESTs were analyzed using TBLASTN at an E-value threshold of < 1e-80. Putative START protein-encoding plant cDNAs were identified from Sputnik EST database  and NCBI using Arabidopsis and rice START proteins as query. START domains were identified from the following plant species: Arabidopsis thaliana, Beta vulgaris, Capsicum annum, Cycas rumphii, Glycine max, Gossypium arboretum, Gossypium hirsutum, Helianthus annuus, Hordeum vulgare, Lactuca sativa, Lotus japonicus, Lycopersicon esculentum, Malus domesticus, Medicago sativa, Medicago truncatula, Oryza sativa, Phalaenopsis sp., Physcomitrella patens, Picea abies, Prunus persica, Solanum tuberosum, Sorghum bicolor, Triticum aestivum, Zea mays and Zinnia elegans.
Signal peptide and transmembrane analysis of the predicted START proteins was accomplished with TargetP  and TMHMM  using default settings. For the human MLN64 protein, and the MLN64-like proteins from D. melanogaster (CG3522/Start1) and C. elegans (3F991), sequence alignments were utilized to confirm the number and coordinates of the transmembrane segments. Experimental data for MLN64  was incorporated to predict the orientation of the transmembrane segments.
Phylogenetic analysis and sequence alignments
Alignments were made using ClustalW  with BLOSUM62 settings. Alignments of START domains were converted into the Molecular Evolutionary Genetics Analysis (MEGA) format in MEGA2.1 [59, 60]. For sequences originating from the same species, alignments were made using PileUp from Genetics Computer Group (GCG), screening for similarities and identities manually. After eliminating redundancies, datasets were realigned with ClustalW and a phylogenetic tree was derived. Phylogenies containing sequences from different organisms were constructed by the neighbor-joining method, assuming that sequences from different lineages evolve at different rates, followed by bootstrapping with 2,000 replicates [60, 61]. The Poisson correction was used to account for multiple substitutions at the same site. Gaps were handled using complete deletion, assuming that different amino-acid sequences evolve differently. Alignments were constructed by GCG's PileUp to examine trends seen in the phylogenies. For analysis of proteins from different organisms, pairwise deletion was used instead of complete deletion. Domain alignments and consensus sequences shown in Figures 3 and 4 were constructed using GCG's Pretty .
Tertiary structure analysis of START domains
Homology modeling of Arabidopsis and rice PCTP-like proteins was accomplished using the SWISS-MODEL web server [63, 64] and DeepView (Swiss-Pdb Viewer) . PDB templates of structures used for homology modeling were PCTP (PDB ID 1LN1A), MLN64 (1EM2A), and StarD4 (1JSSA). The WhatCheck summary reports for the highest likelihood models for Arabidopsis and rice PCTP-like proteins, At3g13062 and ID215312, respectively, were as follows: structure Z-scores: 1st generation packing quality: -1.779, -2.132; 2nd generation packing quality: -4.818, -4.425; Ramachandran plot appearance: -2.674, -2.393; chi-1/chi-2 rotamer normality: 0.138, 1.450; backbone conformation: -1.861, -1.446; and root mean square (RMS) Z-scores: bond lengths: 0.785, 0.859; bond angles: 1.389, 1.243; omega angle restraints: 1.007, 0.960; side chain planarity: 2.374, 1.346; improper dihedral distribution: 1.877, 1.544; inside/outside distribution: 1.324, 1.269. Coordinates and complete WhatCheck Reports are available on request. Cartoon images of the models were produced in RasMol Version 18.104.22.168.
Additional data files
The following additional data are available with the online version of this article: a pdb file (Additional data file 1), the WhatCheck Report (Additional data file 2) and the log (Additional data file 3) for the 3D model of the At3g13062 START domain; a pdb file (Additional data file 4), the WhatCheck Report (Additional data file 5) and the log (Additional data file 6) for the 3D model of the ID215312 START domain; a text file (Additional data file 7) and a MEGA file (Additional data file 8) of the alignment corresponding to the tree in Figure 1; a text file (Additional data file 9) and the MEGA file (Additional data file 10) of the alignment corresponding to the tree in Figure 2a; a text file (Additional data file 11) and the MEGA file (Additional data file 12) of the alignment corresponding to the tree in Figure 2b; a text file (Additional data file 13) and a MEGA file (Additional data file 14) of the alignment corresponding to a tree comparing PCTP START with Arabidopsis and rice START domains, a branch of which is shown in Figure 4a.
Stocco DM: StAR protein and the regulation of steroid hormone biosynthesis. Annu Rev Physiol. 2001, 63: 193-213. 10.1146/annurev.physiol.63.1.193.
Ponting CP, Aravind L: START: a lipid-binding domain in StAR, HD-ZIP and signaling proteins. Trends Biochem Sci. 1999, 24: 130-132. 10.1016/S0968-0004(99)01362-6.
Tsujishita Y, Hurley JH: Structure and lipid transport mechanism of a StAR-related domain. Nat Struct Biol. 2000, 7: 408-414. 10.1038/75192.
Roderick SL, Chan WW, Agate DS, Olsen LR, Vetting MW, Rajashankar KR, Cohen DE: Structure of the human phosphatidylcholine transfer protein in complex with its ligand. Nat Struct Biol. 2002, 9: 507-511.
Tabunoki H, Sugiyama H, Tanaka Y, Fujii H, Banno Y, Jouni ZE, Kobayashi M, Sato R, Maekawa H, Tsuchida K: Isolation, characterization, and cDNA sequence of a carotenoid binding protein from the silk gland of Bombyx mori larvae. J Biol Chem. 2002, 277: 32133-32140. 10.1074/jbc.M204507200.
Hanada K, Kumagai K, Yasuda S, Miura Y, Kawano M, Fukasawa M, Nishijima M: Molecular machinery for non-vesicular trafficking of ceramide. Nature. 2003, 426: 803-809. 10.1038/nature02188.
Romanowski MJ, Soccio RE, Breslow JL, Burley SK: Crystal structure of the Mus musculus cholesterol-regulated START protein 4 (StarD4) containing a StAR-related lipid transfer domain. Proc Natl Acad Sci USA. 2002, 99: 6949-6954. 10.1073/pnas.052140699.
Iyer L, Koonin EV, Aravind L: Adaptations of the helix-grip fold for ligand binding and catalysis in the START domain superfamily. Proteins. 2001, 43: 134-144. 10.1002/1097-0134(20010501)43:2<134::AID-PROT1025>3.3.CO;2-9.
The RCSB Protein Data Bank. [http://www.rcsb.org/pdb]
SCOP: structural classification of proteins. [http://scop.mrc-lmb.cam.ac.uk/scop]
Soccio RE, Breslow JL: StAR-related lipid transfer (START) proteins: mediators of intracellular lipid metabolism. J Biol Chem. 2003, 278: 22183-22186. 10.1074/jbc.R300003200.
Ennis HL, Dao DN, Pukatzki SU, Kessin RH: Dictyostelium amoebae lacking an F-box protein form spores rather than stalk in chimeras with wild-type. Proc Natl Acad Sci USA. 2000, 97: 3292-3297. 10.1073/pnas.050005097.
Nelson MK, Clark A, Abe T, Nomura A, Yadava N, Funair CJ, Jermyn KA, Mohanty S, Firtel RA, Williams JG: An F-Box/WD-40 repeat-containing protein important for Dictyostelium cell-type proportioning, slug behaviour, and culmination. Dev Biol. 2000, 224: 42-59. 10.1006/dbio.2000.9793.
Tekinay T, Ennis HL, Wu MY, Nelson M, Kessin RH, Ratner DI: Genetic interactions of the E3 ubiquitin ligase component FbxA with cyclic AMP metabolism and a histidine kinase signaling pathway during Dictyostelium discoidum development. Eukaryotic Cell. 2003, 2: 618-626. 10.1128/EC.2.3.618-626.2003.
Riechmann JL: Transcriptional regulation: a genomic overview. In The Arabidopsis Book. Edited by: Somerville CR, Meyerowitz EM. 2002, Rockville, MD: American Society of Plant Biologists, [http://www.aspb.org/downloads/arabidopsis/reichm.pdf]
Szymanski DB, Jilk RA, Pollack SM, Marks MD: Control of GL2 expression in Arabidopsis leaves and trichomes. Development. 1998, 125: 1161-1171.
Zhong R, Ye Z-H: IFL1, a gene regulating interfascicular fiber differentiation in Arabidopsis, encodes a homeodomain-leucine zipper protein. Plant Cell. 1999, 11: 2139-2152. 10.1105/tpc.11.11.2139.
Ohashi Y, Oka A, Rodrigues-Pousada R, Possenti M, Ruberti I, Morelli G, Aoyama T: Modulation of phospholipid signaling by GLABRA2 in root-hair pattern formation. Science. 2003, 300: 1427-1430. 10.1126/science.1083695.
Abe M, Takahashi T, Komeda Y: Identification of a cis-regulatory element for L1 layer-specific gene expression, which is targeted by an L1-specific homeodomain protein. Plant J. 2001, 26: 487-494. 10.1046/j.1365-313x.2001.01047.x.
Abe M, Katsumata H, Komeda Y, Takahashi T: Regulation of shoot epidermal cell differentiation by a pair of homeodomain proteins in Arabidopsis. Development. 2003, 130: 635-643. 10.1242/dev.00292.
Kikuchi S, Satoh K, Nagata T, Kawagashira N, Doi K, Kishimoto N, Yazaki J, Ishikawa M, Yamada H, Ooka H, et al: Collection, mapping and annotation of over 28,000 cDNA clones from japonica rice. Science. 2003, 301: 376-379. 10.1126/science.1081288.
Schena M, Davis RW: Structure of homeobox-leucine zipper genes suggests a model for the evolution of gene families. Proc Natl Acad Sci USA. 1994, 91: 8393-8397.
Yang J-Y, Chung M-C, Tu C-Y, Leu W-M: OSTF1 : A HD-GL2 family homeobox gene is developmentally regulated during early embryogenesis in rice. Plant Cell Physiol. 2002, 43: 628-638. 10.1093/pcp/pcf076.
Di Cristina M, Sessa G, Dolan L, Linstead P, Baima S, Ruberti I, Morelli G: The Arabidopsis Athb-10 (GLABRA2) is an HD-Zip protein required for regulation of root hair development. Plant J. 1996, 10: 393-402. 10.1046/j.1365-313X.1996.10030393.x.
Landschulz WH, Johnson PF, McKnight SL: The leucine zipper: a hypothetical structure common to a new class of DNA binding proteins. Science. 1988, 240: 1759-1764.
O'Shea EK, Rutkowski R, Kim PS: Evidence that the leucine zipper is a coiled coil. Science. 1989, 243: 538-542.
Ito M, Sentoku N, Nishimura A, Hong S-K, Sato Y, Matsuoka M: Position dependent expression of GL2-type homeobox gene, Roc1: significance for protoderm differentiation and radial pattern formation in early rice embryogenesis. Plant J. 2002, 29: 497-507. 10.1046/j.1365-313x.2002.01234.x.
Kubo H, Peeters AJM, Aarts MGM, Pereira A, Koornnneef M: ANTHOCYANINLESS2, a homeobox gene affecting anthocyanin distribution and root development. Plant Cell. 1999, 11: 1217-1226. 10.1105/tpc.11.7.1217.
Rerie WG, Feldmann KA, Marks MD: The GLABRA2 gene encodes a homeo domain protein required for normal trichome development in Arabidopsis. Genes Dev. 1994, 8: 1388-1399.
McConnell JR, Emery J, Eshed Y, Bao N, Bowman J, Barton K: Role of PHABULOSA and PHAVOLUTA in determining radial patterning in shoots. Nature. 2001, 411: 709-713. 10.1038/35079635.
Otsuga D, DeGuzman B, Prigge MJ, Drews GN, Clark SE: REVOLUTA regulates meristem initiation at lateral positions. Plant J. 2001, 25: 223-236. 10.1046/j.1365-313x.2001.00959.x.
Talbert PB, Adler HT, Parks DW, Comai L: The REVOLUTA gene is necessary for apical meristem development and for limiting cell divisions in the leaves and stems of Arabidopsis thaliana. Development. 1995, 121: 2723-2735.
Ratcliffe OJ, Riechmann JL, Zhang JZ: INTERFASCICULAR FIBERLESS1 is the same gene as REVOLUTA. Plant Cell. 2000, 12: 315-317. 10.1105/tpc.12.3.315.
Baima S, Nobili F, Sessa G, Lucchetti S, Ruberti I, Morelli G: The expression of the Athb-8 homeobox gene is restricted to provascular cells in Arabidopsis thaliana. Development. 1995, 121: 4171-4182.
Baima S, Possenti M, Matteucci A, Wisman E, Altamura MM, Ruberti I, Morelli G: The Arabidopsis ATHB-8 HD-Zip protein acts as a differentiation-promoting transcription factor of the vascular meristems. Plant Physiol. 2001, 126: 643-655. 10.1104/pp.126.2.643.
Lemmon MA, Ferguson KM: Molecular determinants in pleckstrin homology domains that allow specific recognition of phosphoinositides. Biochem Soc Trans. 2001, 29: 377-384.
Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, et al: The Pfam contribution to the annual NAR database issue. Nucleic Acids Res. 2004, 32 Database issue: D138-D141. 10.1093/nar/gkh121.
Alpy F, Stoeckel ME, Dierich A, Escola JM, Wendling C, Chenard MP, Vanier MT, Gruenberg J, Tomasetto C, Rio MC: The steroidogenic acute regulatory protein homolog MLN64, a late endosomal cholesterol-binding protein. J Biol Chem. 2001, 276: 4261-4269. 10.1074/jbc.M006279200.
Moog-Lutz C, Tomasetto C, Regnier CH, Wendling C, Lutz Y, Muller D, Chenard MP, Basset P, Rio MC: MLN64 exhibits homology with the steroidogenic acute regulatory protein (StAR) and is over-expressed in human breast carcinomas. Int J Cancer. 1997, 71: 183-191. 10.1002/(SICI)1097-0215(19970410)71:2<183::AID-IJC10>3.0.CO;2-J.
Roth GE, Gierl MS, Vollborn L, Meise M, Lintermann R, Korge G: The Drosophila gene Start1: a putative cholesterol transporter and key regulator of ecdysteroid synthesis. Proc Natl Acad Sci USA. 2004, 101: 1601-1606. 10.1073/pnas.0308212100.
Feng L, Chan WW, Roderick SL, Cohen DE: High-level expression and mutagenesis of recombinant human phosphatidyl transfer protein using a synthetic gene: evidence for a C-terminal membrane binding domain. Biochemistry. 2000, 39: 15399-15409. 10.1021/bi001076a.
Choinowski T, Hauser H, Piontek K: Structure of sterol carrier protein 2 at 1.8 Å resolution reveals a hydrophobic tunnel suitable for lipid binding. Biochemistry. 2000, 39: 1897-1902. 10.1021/bi992742e.
Lee JY, Min K, Cha H, Shin DH, Hwang KY, Suh SW: Rice non-specific lipid transfer protein: The 1.6 Å crystal structure in the unliganded state reveals a small hydrophobic cavity. J Mol Biol. 1998, 276: 437-448. 10.1006/jmbi.1997.1550.
Kallen C, Billheimer JT, Summers SA, Stayrook SE, Lewis M, Strauss JF: Steroidogenic acute regulatory protein (StAR) is a sterol transfer protein. J Biol Chem. 1998, 273: 26285-26288. 10.1074/jbc.273.41.26285.
Laudet V, Gronemeyer H: The Nuclear Receptors FactsBook. 2002, London: Academic Press
NCBI BLAST. [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/blast]
Arabidopsis thaliana sequence BLAST search. [http://tigrblast.tigr.org/euk-blast]
MATDB BLAST query. [http://mips.gsf.de/proj/thal/db/search/blast_arabi.html]
Karlowski WM, Schoof H, Janakiraman V, Stuempflen V, Mayer KF: MOsDB: an integrated information resource for rice genomics. Nucleic Acids Res. 2003, 31: 190-192. 10.1093/nar/gkg073.
InterPro: home. [http://www.ebi.ac.uk/interpro]
NCBI CD-search. [http://www.ncbi.nih.gov/Structure/cdd/wrpsb.cgi]
Marchler-Bauer A, Anderson JB, DeWeese-Scott C, Fedorova ND, Geer LY, He S, Hurwitz DI, Jackson JD, Jacobs AR, Lanczycki CJ, et al: CDD: a curated Entrez database of conserved domain alignments. Nucleic Acids Res. 2003, 31: 383-387. 10.1093/nar/gkg087.
TAIR WU-BLAST 2.0 form. [http://www.arabidopsis.org/wublast/index2.jsp]
TAIR: microarray elements search and download. [http://www.arabidopsis.org/tools/bulk/microarray/index.html]
The Sputnik resource for the annotation of clustered plant ESTs. [http://mips.gsf.de/proj/sputnik]
TargetP server v1.01. [http://www.cbs.dtu.dk/services/TargetP]
TMHMM server v.2.0. [http://www.cbs.dtu.dk/services/TMHMM]
Molecular Evolutionary Genetics Analysis MEGA. [http://www.megasoftware.net]
Kumar S, Tamura K, Jakobsen IB, Nei M: MEGA2: molecular evolutionary genetics analysis software. Bioinformatics. 2001, 17: 1244-1245. 10.1093/bioinformatics/17.12.1244.
Nei M, Kumar S: Molecular Evolution and Phylogenetics. 2000, Oxford: Oxford University Press
Devereux J, Haeberli P, Smithies O: A comprehensive set of sequence analysis programs for the VAX. Nucleic Acids Res. 1984, 12: 387-395.
Schwede T, Kopp J, Guex N, Peitsch MC: SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 2003, 31: 3381-3385. 10.1093/nar/gkg520.
Geux N, Peitsch MC: SWISS-MODEL and the Swiss-Pdb Viewer: an environment for comparative protein modeling. Electrophoresis. 1997, 18: 2714-2723.
The authors are grateful to Animesh Ray for the generous use of his facilities at the Keck Graduate Institute. We thank David Wild for help with protein modeling and Shinhan Shiu of the Department of Ecology and Evolution at the University of Chicago for searching the Chlamydomonas reinhardtii genome database at the Department of Energy Joint Genome Institute http://www.jgi.doe.gov for START domains. We thank Elliot Meyerowitz, Animesh Ray, and David Wild for valuable discussions and critical reading of the manuscript, and Raymond Soccio for many helpful suggestions. This work was supported by the National Science Foundation (Animesh Ray: NSF EIA-0205061 and NSF EIA-0130059), a grant from the Deutsche Forschungsgemeinschaft (SCHR 744/1-2), and a research fellowship to K.S. (SCHR 744/1-1,1-3).
Electronic supplementary material
About this article
Cite this article
Schrick, K., Nguyen, D., Karlowski, W.M. et al. START lipid/sterol-binding domains are amplified in plants and are predominantly associated with homeodomain transcription factors. Genome Biol 5, R41 (2004) doi:10.1186/gb-2004-5-6-r41
- Leucine Zipper
- Transmembrane Segment
- Carboxy Terminus
- Start Gene
- Sterile Alpha Motif