Identification and characterization of LIM gene family in Brassica rapa
BMC Genomics volume 15, Article number: 641 (2014)
LIM (Lin-11, Isl-1 and Mec-3 domains) genes have been reported to trigger the formation of actin bundles, a major higher-order cytoskeletal assembly, in higher plants; however, the stress resistance related functions of these genes are still not well known. In this study, we collected 22 LIM genes designated as Brassica rapa LIM (BrLIM) from the Brassica database, analyzed the sequences, compared them with LIM genes of other plants and analyzed their expression after applying biotic and abiotic stresses in Chinese cabbage.
Upon sequence analysis these genes were confirmed as LIM genes and found to have a high degree of homology with LIM genes of other species. These genes showed distinct clusters when compared to other recognized LIM proteins upon phylogenetic analysis. Additionally, organ specific expression of these genes was observed in Chinese cabbage plants, with BrPLIM2a, b, c, BrDAR1, BrPLIM2e, f and g only being expressed in flower buds. Furthermore, the expression of these genes (except for BrDAR1 and BrPLIM2e) was high in the early flowering stages. The remaining genes were expressed in almost all organs tested. All BrDAR genes showed higher expression in flower buds compared to other organs. These organ specific expressions were clearly correlated with the phylogenetic grouping. In addition, BrWLIM2c and BrDAR4 responded to Fusarium oxysporum f. sp. conglutinans infection, while commonly two BrDARs and eight BrLIMs responded to cold, ABA and pH (pH5, pH7 and pH9) stress treatments in Chinese cabbage plants.
Taken together, the results of this study indicate that BrLIM and BrDAR genes may be involved in resistance against biotic and abiotic stresses in Brassica.
LIM proteins have mainly been described in animals in which they play key roles in a variety of important developmental processes [1, 2]. These proteins contain one or several (up to five) double zinc finger motifs, known as LIM domains, which function by mediating protein-protein interactions [3, 4]. The three transcription factors, LIN11 from Caenorhabditis elegans, ISL-1 from rat  and MEC-3 from C. elegans, from which the acronym LIM is derived, operate during transcription in association with a DNA-binding homeodomain. LIM domains have since been found in a wide variety of eukaryotic proteins of diverse functions. LIM domain-containing transcription factors without the homeodomain have also been described. Specifically, the LIM-only protein LMO2 was found to act as a bridging molecule in assembly of the erythroid DNA-binding complex , while other LIM-domain-containing proteins such as LIM kinases are known to participate in regulation of actin dynamics through phosphorylation of cofilin [9, 10]. In contrast, plants possess two distinct sets of LIM proteins, one that is plant–specific and has been partially functional characterized [11, 12]. Another, cysteine rich protein (CRP)-like that comprises CRPs exhibiting the same overall structure found in animals (i.e., two very similar LIM domains separated by an ≈ 50 amino acid–long interLIM domain and a relatively short and variable C–terminal domain) .
Plants have evolved mechanisms to exploit, survive, or minimize the negative impact of a diverse range of environmental factors including potential pathogens, and in many cases the plant cytoskeleton is instrumental in mediating the plant’s response. Changes in organization of the plant cytoskeleton during plant interactions with these factors are complex and varied, and much still remains to be elucidated, especially in terms of the molecules that signal and induce the dramatic reorganizations that are often observed. Physical and chemical barriers resulting from actin-dependent cytoplasmic aggregation, secretion, and papilla formation appear to constitute an important and probably ancient form of basal resistance to pathogen attack . Reorganization of cytoskeletal elements may also be a mechanism responsible for abiotic stress responses. For example, cytoskeletal elements translocate chloroplasts under high light conditions , facilitate gravity sensing , and direct cellular response to wounding [17, 18]. A full understanding of the cytoskeletal basis of cytoplasmic aggregation will require identification and characterization of proteins involved in the signaling pathways those induce cytoskeletal rearrangement and responsible proteins for bringing about cytoskeletal reorganization and their function. This mechanism utilizes regulatory proteins that generate intracellular motility, such as myosin  and thick microfilament bundles focusing on the site of infection, such as the actin-bundling proteins villin and fimbrin [20, 21].
In the cytoplasm, actin monomers polymerize into actin filaments (AFs), which constitute the core elements of the actin cytoskeleton, providing mechanical support to the cytoplasm and serving as tracks for myosin-dependent intracellular transport [22, 23]. To fulfill its various roles, the actin cytoskeleton requires a sophisticated regulatory system to control its organization and dynamics at both spatial and temporal levels. Primary components of this system are the actin binding proteins (ABPs) that directly interact with monomeric and/or polymerized actin to promote AF nucleation, polymerization, depolymerization, stabilization, severing, capping, and cross-linking . Recently, a number of vertebrate LIM domain proteins belonging to the Cys-rich protein (CRP) family and several structurally related plant proteins (hereafter referred to as plant LIMs) have been shown to function as ABPs [25–28]. In Arabidopsis, two distinct LIM gene subfamilies have been reported. The first subfamily contains six genes, which encode proteins with two LIM domains that share similarity with animal CRPs . These LIM domain-containing proteins bind actin filaments and influence actin cytoskeleton organization . In vitro, chicken CRP1 and tobacco (Nicotiana tabacum) WLIM1 bind directly to AFs to trigger the formation of thick actin bundles. Importantly, over-expression of CRP1 and WLIM1 proteins has been shown to be sufficient to significantly increase the bundling of AFs in rat fibroblasts and tobacco cells, respectively [13, 26, 27]. The genes in the second subfamily encode DA1 and DA1 related (DAR) proteins with a single LIM domain and are specific to plants . DA1 encodes a putative ubiquitin receptor with two ubiquitin interacting motifs (UIMs) and a single LIM domain, while e.g. DAR2 contains a single LIM domain and a putative zinc-binding domain but lacks UIMs suggesting a possible functional divergence with DA1 . The members of this family are involved in diverse functions e.g. seed and organ size control, root meristem size control, resistance signaling and cold response [11, 30, 31].
In this study, we retrieved 22 Brassica rapa LIM genes from the Brassica database, analyzed the sequences and studied their homology with LIM genes of other species. We also analyzed their phylogenetic classification and discussed their correlation with organ specific expression. Expression of these genes was also analyzed following application of different biotic and abiotic stresses and their association was discussed with biotic and abiotic stress resistance in Brassica.
Results and discussion
Sequence analysis and phylogenetic classification
We identified 22 LIM proteins, designated as B. rapa LIM (BrLIM), after which the sequences were analyzed (Table 1). The predicted size of the 22 BrLIMs ranged from 92 to 1315 amino acids (10.36 to 150.61 kDa), and the predicted isoelectric points varied from 5.28 to 9.77. Analysis of the protein domain organization showed that proteins contained the characteristic LIM domain in the conserved region. Additionally, genomic DNA sequence of the 22 BrLIM genes was determined from the B. rapa chromosome sequences and the introns and exons were identified by sequence analysis. Structural information pertaining to the 22 BrLIM genes is presented in Additional file 1: Figure S1. This analysis confirmed that the identified genes code for LIM domain containing proteins. A BLAST search of the NCBI database conducted for comparison of these BrLIMs with LIM genes of other species revealed that the deduced amino acid sequences of 22 BrLIMs shared high homology, primarily with LIMs of A. thaliana and some other homologous species (Table 2). Pairwise amino acid sequence comparisons among these 22 BrLIM proteins were also calculated using BLAST of the NCBI database and shown in Additional file 2: Table S2. The similarity within these BrLIM protein sequences ranged from 23 to 99% and 8 BrLIMs showed above 90% similarity within the species indicating their probable duplication.
We also retrieved all LIM family protein sequences of A. thaliana, Populus trichocarpa and Oryza sativa from NCBI and constructed a phylogenetic tree with 22 deduced amino acid sequences of BrLIM using the NJ method (Figure 1, Additional file 3: Table S3). In the phylogenetic tree, previously identified LIM1 and LIM2 groups  were clearly separated and the plant LIM family was further divided into four groups, αLIM1, βLIM1, γLIM2, and δLIM2, resulting from division of the LIM1 and LIM2 groups. This phylogenetic analysis confirms the existence of the PLIM1, WLIM1, PLIM2, and WLIM2 subgroups as previously described . The αLIM1 group includes the PLIM1 and WLIM1 subgroups, while βLIM1 is a newly identified group . The WLIM2 and PLIM2 subgroups belong to the γLIM2 and δLIM2 groups, respectively. Among the 22 BrLIMs, BrLIM10 and 11 were clustered with WLIM1 subgroup under the αLIM1 group. BrLIM5, 6, 7, 8 and 9 were separated with the WLIM2 subgroup under the γLIM2 group. Moreover, BrLIM1, 2, 3, 4, 20, 21 and 22 separated with the PLIM2 subgroup under the δLIM2 group, showing close relationships with the LIM sequences of other species. Multiple alignments among members of the WLIM1, WLIM2, PLIM2 and DAR groups also showed clearly distinct amino acid residues conserved in the domains and other regions of their protein sequences (Figure 2). In Arabidopsis, LIM proteins are encoded by genes from two distinct subfamilies, one encoding proteins with two LIM domains that are homologous to animal CRP proteins . Consequently, highly homologous WLIM1, PLIM2, and WLIM2 proteins of B. rapa of this subfamily also contain two LIM domains except BrLIM6. Interestingly, with two LIM domains N-terminal extensions were observed in BrLIM7 and 9 protein sequences and BrLIM9 had additional 11 transmembrane domains in the N-terminal extension (Figure 2, Table 1). These N-terminal extensions suggest probable misannotation or divergent functions of these two genes which require further experimental evidence.
To investigate the existence of N-terminal extensions of BrLIM7 and 9, we analyzed expression of both genes over different organs specifying forward and reverse primers at several locations of the sequences (Additional file 1: Figure S2). Both genes showed expression, when forward (#5 of BrLIM7 and #7 of BrLIM9) and reverse (#6 of BrLIM7 and #8 of BrLIM9) primers were selected at the conserved parts. Conversely, no expression was found, when the forward primers (#1 of BrLIM7 and #3 of BrLIM9) were selected at the extended parts and the reverse primers (#3 and #4 of BrLIM7; #5 and #6 BrLIM9) at several locations of the conserved parts. Subsequently, 3’-RACE PCR using primer #1 of BrLIM7 and primer #3 of BrLIM9 also showed discontinuation of the transcription towards conserved parts and found to show polyA tails indicating the end of the transcript. Alongside, genomic DNA amplification specifying these transition portions of both genes revealed no sequence distortion. Thus, it is clear that the extended portions are not part of BrLIM7 and 9. It was a wrong annotation of the splice sites in the original prediction. The correct exons of BrLIM7 and 9 were then drawn in the postulated gene structures (Additional file 1: Figure S2) and the primer combinations, #5 and #6 of BrLIM7 and #7 and #8 of BrLIM9 located at the corrected exons, were used in latter expression studies. These corrected sequences were deposited in the genbank under the accession numbers, KJ686594 (BrLIM7) and KJ686595 (BrLIM9). These corrected sequences are the basis for all further analysis and all figures and tables were corrected accordingly.
The genes in the second subfamily encode DA1 and DA1-related (DAR) proteins specific to plants contain a single conserved LIM domain . Interestingly, BrLIM12, 13, 14, 15, 16, 17, 18 and 19 were clustered into a new group, showing close relationships with eight DA1 and DAR sequences of A. thaliana and four DAR sequences of O. sativa. Additionally, all members of this new group contain only one LIM domain, while some also contain one or more UIM(s). Therefore, BrLIM proteins were renamed according to their phylogenetic relationship as B. rapa DAR (BrDAR) (Table 2). Alignment of the motif sequences of the DA1 and DAR group showed very low similarity with WLIM1, WLIM2 and PLIM2 sequences, indicating their distinctness from these groups (Figure 2). However, the BrDARs were highly similar to DA1 and DAR sequences of A. thaliana and O. sativa, and these all contained similar patterns of motifs in their conserved regions. DA1 of A. thaliana contains two UIMs, a single LIM domain and the highly conserved C-terminal region . In contrast, DAR2 contains a single LIM domain and a putative zinc-binding domain but lacks UIMs compared with DA1 and DAR1 . Similarly, BrDAR4 and 5 contain two UIMs and a single LIM domain, while BrDAR8 contains four UIM and a single LIM domain with highly conserved C-terminal region in their protein sequences. Other BrDARs contain only a single LIM domain and highly conserved C-terminal regions (Table 1). Thus, the BrLIM genes showed high similarity with corresponding groups of LIM genes of different species, indicating their close evolutionary relationship.
Organ specific expression analysis under unstressed condition
Expression analysis was performed using specific primers with equal amounts of cDNA templates prepared from the mRNA of roots, stems, leaves and flower buds of healthy Chinese cabbage plants without applying any stresses by RT-PCR and organ specific expression was observed. Specifically, BrPLIM2d, BrWLIM2a, b, c, d, e, BrWLIM1a and b showed expression in all the organs tested. BrPLIM2a, b, c, BrPLIM2e, f and g were only expressed in flower buds (Figure 3A). Expression of the genes that were only expressed in flower buds was further analyzed during six developmental stages (from the very young bud stage to the mature bud stage), and all genes except BrPLIM2e showed higher expression in young buds (Figure 3B). This expression pattern and phylogenetic classification revealed two clear groups of BrLIM genes, with one group (BrPLIMs except BrPLIM2d) only being expressed in flower buds and another group (WLIM) being expressed in all vegetative parts of the Chinese cabbage plants tested in this study. Thus, organ specific expressions of BrLIMs showed clear correlations with phylogenetic grouping as well.
DAR subfamily genes of B. rapa showed a different expression pattern compared to other LIM domain containing genes. Specifically, BrDAR3, 4 and 5 showed expression in all the organs tested. BrDAR2, 7 and 8 showed expression in all organs except roots, while BrDAR1 was only expressed in flower buds. Noticeably, all BrDAR genes expressed highly in flower buds compared to other organs (Figure 3A). Expression of BrDAR1 that was only expressed in flower buds was further analyzed during six developmental stages (from the very young bud stage to the mature bud stage), and found higher expressions in young and middle aged buds (Figure 3B). Microarray expression of DA1 and DAR genes of A. thaliana were also lower at different vegetative parts compared to floral organs and comparatively higher at young flower bud stages (Figure 3C).
The PLIM-1 protein of sunflower appeared in the microspore stage in a limited number of cytoplasmic bodies, became undetectable in bicellular pollen, and reappeared in tricellular pollen grains in cortical patches, being particularly concentrated in the F-actin-enriched germination cones of the vegetative cell. Taken together, these findings suggest a dual function of PLIM-1 gene during pollen development . Conversely, expression of NtPLIM1 of tobacco in germinating pollen suggests a possible role in pollen germination and/or pollen tube growth . LiLIM1, an effective actin bundling protein, also plays an important role in pollen tube growth of Lilium longiflorum. Likewise, the BrPLIMs expressed in young flower bud stage might play a role in the development of pollen tube. Earlier studies have also suggested the existence of two main LIM gene subfamilies that differ in their expression patterns . Again, based on phylogenetic analysis of 149 LIMs and comparison of the available expression data, a complex classification of LIM genes has been proposed by Arnaud et al., . According to this classification, the Arabidopsis LIM gene family comprises three vegetative (WLIM1 and WLIM2a and b) and three reproductive (PLIM2a-c) isoforms, of which, all PLIMs have been shown to be expressed in pollen [41, 42]. Ye et al.  also suggested the critical roles of Arabidopsis PLIMs in actin configuration during pollen germination and tube growth. Specific control of PLIM actin regulatory activities by pH is particularly relevant with regard to the potential biological functions of these proteins in growing tip region of pollen tubes .
Expression analysis under stress conditions
Biotic stress responses are a crucial issue for the Brassica crops and functions of LIM genes in biotic stress responses are not well studied yet. We investigated the responses of the 22 BrLIM genes identified in this study to infection by F. oxysporum f.sp. conglutinans in Chinese cabbage after various times. BrWLIM2c and BrDAR4 were found to be up-regulated with the time course of infection. Specifically, BrDAR4 showed very high responses in the late stages of fungal infection, while BrWLIM2c responded during earlier stages (Figure 4). AtWLIM1, AtWLIM2b, AtDAR7, AtDAR4 and AtDAR1 genes also showed responses during microarray expression analysis after infection with F. oxysporum in A. thaliana, which were collected from public database (Additional file 4: Figure S3). We also analyzed the expression of all BrLIMs after infection with the soft rot disease causing necrotroph bacteria, Pectobacterium carotovorum subsp. carotovorum; however, no differential expression was observed in Chinese cabbage plants (data not shown). These results suggest probable involvement of BrWLIM2c and BrDAR4 with F. oxysporum f.sp. conglutinans resistance and other BrLIM and BrDAR genes might also have association with resistance against other biotic stress agents in Brassica crops.
Plants encounter and must deal with a range of other organisms that may be potential partners or pathogens, and the plant cytoskeleton plays a key role in their response to such organisms. Cell wall appositions, or papillae, are important barriers formed by plants in defense against attempted penetration by fungal and oomycete pathogens . Prior to the development of papillae, plant cytosol and subcellular components are rapidly translocated to the site of pathogen penetration via a mechanism known as cytoplasmic aggregation that is dependent on the action of the actin component of the cytoskeleton [44, 45]. This cytoplasmic aggregation has been observed in many plant-microbe interactions  and is a common resistance response to pathogens by both dicotyledonous and monocotyledonous plants to invading filamentous pathogens. Microfilaments are the main cytoskeletal element responsible for penetration resistance based on cytoplasmic aggregation and papilla formation at the site of infection. For example, focusing of actin microfilaments at the penetration site is more extensive in resistant plants than susceptible hosts in flax (Linum usitatissimum)-Melampsora and barley (Hordeum vulgare)-Erysiphe interactions [47, 48], suggesting that accumulation of material at the site of infection is related to the degree of host resistance. Some other previous observations suggest that inhibition of penetration through cytoplasmic aggregation and papilla formation is an early, if not the first, tactic in plant resistance, and that it may be backed up by the hypersensitive response . Thus, the physical and chemical barrier resulting from actin-dependent cytoplasmic aggregation, secretion, and papilla formation appears to constitute an important and probably ancient form of basal resistance to pathogen attack . For a full understanding of the cytoskeletal basis of cytoplasmic aggregation, it is very essential to elucidate the role of responsible regulatory proteins that participate in other forms of cytoplasmic streaming, as well as additional factors specific for the localized defense response. Proteins those generate intracellular motility (e.g., myosin) are likely to be important components of such regulations . In addition, several CRP-like plant LIM proteins have been shown to function as ABPs [25–28], and LiLIM1 has been proposed as an effective ABP, that plays an important role in pollen tube growth of L. longiflorum. In vitro, tobacco (N. tabacum) WLIM1 binds directly to AFs and triggers the formation of thick actin bundles. Importantly, overexpression of WLIM1 proteins was sufficient to significantly increase the bundling of AFs in tobacco cells [13, 27]. Thus, actin binding and bundling function of LIM genes could also be an important factor for actin-dependent resistance mechanism against pathogen infection. In this study, high expression of BrWLIM2c gene against pathogen infection is tempting to speculate that BrLIM protein may also be implicated in the perception of pathogen resistance mechanism.
Genes in the second subfamily encode DA1 proteins, which are specific to plants. The Arabidopsis LIM domain protein DA1 was characterized to function as an ubiquitin receptor . Ubiquitin is a highly conserved and wide-spread small protein modifier that is engaged in a wide range of cellular  and biological processes, such as resistance to disease and abiotic stresses [51–53]. The consequences of ubiquitination are accusing the target protein to proteolysis or relocalization or endocytosis . DA1 is inevitably involved in ubiquitination resulting extension of the cellular proliferation period, thus increasing cell numbers and ultimately plant organ size [11, 12]. Our highly expressed BrDAR4 in response to pathogen infection might be involved in biotic stress resistance of B. rapa following aforesaid mechanism. Another member of this family CHS3/AtDAR4 containing TIR-NBS-LRR domain is associated with resistance signaling and cold response  and AtDAR2 with controlling of root meristem size . Alongside, no BrDAR protein contained TIR-NBS-LRR domain like as CHS3/AtDAR4. However, reports on the responses of LIM genes to biotic stresses are still lacking, and this is the first report of such possible role of LIM domain containing genes against biotic stress caused by F. oxysporum f.sp. conglutinans in Chinese cabbage. Thus, BrLIM and BrDAR genes might be involved in biotic stress resistance related functions in B. rapa.
In addition to biotic stresses, expression of BrLIMs was also analyzed in response to abiotic stress conditions such as cold, ABA and pH stresses. Commonly ten BrLIMs among them, showed differential expressions in response to the abiotic stresses applied in this study. Specifically, BrPLIM2d, BrWLIM2a, b, c, d, e, BrWLIM1a, b, BrDAR3 and 4 showed differential expression after cold, ABA and pH stress treatments compared to mock treated samples (control) of Chinese cabbage (Figure 5). BrPLIM2d, BrWLIM2a, b and BrDAR4 were down-regulated in response to cold stress throughout the treatment period, while BrWLIM2c, d, e, BrWLIM1a, b and BrDAR3 were up-regulated during the first 8 h of treatment, after which they were down-regulated. Following ABA stress treatments, BrPLIM2d, BrWLIM2b, e and BrDAR4 were also down-regulated throughout the treatment period, while BrWLIM2a, c, d, BrWLIM1a, b and BrDAR3 were up-regulated during the first 8 h of treatment, after which they were down-regulated. Moreover, BrPLIM2d, BrWLIM2a and BrDAR4 were down-regulated and BrWLIM1a, b BrWLIM2b and c were up-regulated throughout the treatment period in response to pH 5, 7 and 9 stress treatments compared to mock treated samples (control) which were maintained after transferring into a fresh medium with pH 5.8. BrWLIM2d and BrWLIM2e were down-regulated in response to pH 5 and 9, while was up-regulated in response to pH 7. BrDAR3 was down-regulated in response to pH 7 and up-regulated in response to pH 5 and 9. Interestingly, BrPLIM2d and BrDAR4 were down-regulated and BrWLIM1a, b and BrWLIM2c were up-regulated in response to all abiotic stress treatments applied in this study. AtWLIM1, AtWLIM2a, AtWLIM2b, DAR7, DAR4, DAR2 and DAR1 genes also showed responses after cold stress treatments in A. thaliana during microarray expressions study collected from public database (Additional file 4: Figure S4).
Plants actin cytoskeleton is a major signaling target and in response to numerous abiotic and biotic stimuli, it changes dramatically in various ways, ranging from filament bundling, to massive actin depolymerization, to assembly of new filament arrays [55–57]. Plant cells can detect mechanical stress on the epidermal cell surface and respond by reorganizing subcellular components like, actin microfilaments, endoplasmic reticulum (ER) and peroxisomes in a manner similar to that induced during attack by potential fungal or oomycete pathogens . Expression analysis of poplar LIMs in tension wood also suggested a connection between plant LIMs and mechanical stress . Thus, rearrangement of cytoskeletal elements through actin binding and bundling may also be the mechanism of abiotic stress responses of the CRP-like LIM genes. Alongside, aforementioned ubiquitination mechanism of DA1 and related proteins might also be implicated in abiotic stress resistance function of BrDAR proteins. In a very recent study, expression of GmaDA1 gene of soybeans (Glycine max) showed responses to salt, drought, acid and alkali stresses and ABA treatment . AtDA1 expression occurs in response to ABA trteatment  and CHS3/AtDAR4 is associated with resistance signaling and cold response in A. thaliana. Supporting this report Bi et al.  hypothesized that disease resistance and environmental stress response pathways are intricately tied in plants, although the degree of overlap and the precise signaling events required for each pathway remain unclear. Thus, the results of this study indicate the likely involvement of the analyzed genes towards abiotic stress tolerance in Brassica.
This present study identified 22 LIM domain containing genes in B. rapa and classified into two distinct functional groups. One group is plant–specific and designated as BrDAR, and another is CRP-like designated as BrLIM. These sequences were analyzed, compared them with LIM domain containing genes of other plants and analyzed their expression after applying biotic and abiotic stresses in Chinese cabbage. Overall, BrPLIM2d, BrWLIM2a, b, c, d, e, BrWLIM1a, b, BrDAR3 and 4 responded after cold, ABA and pH stress treatments in B. rapa, of which BrWLIM2c and BrDAR4 showed responses to F. oxysporum f.sp. conglutinans infection as well. Thus, these eight BrLIMs and two BrDARs might be promising candidates for multiple abiotic stresses resistance, while BrWLIM2c and BrDAR4 might be the most promising among them for both biotic and abiotic stresses resistance in Brassica crops. The remaining BrLIMs and BrDARs might also be useful resources for induction of resistance to other biotic and abiotic stresses and/or various developmental processes of Brassica.
Brassica rapa morphotype Chinese cabbage cv. SUN-3061 plants were grown at the Department of Horticulture, Sunchon National University, Korea. Fresh roots, stems, leaves and flower buds of the Chinese cabbage plants were harvested by freezing them in liquid nitrogen and stored at -80°C to carry out subsequent organ specific expression studies.
Chinese cabbage plants were infected with F. oxysporum f.sp. conglutinans at the Screening Center for Disease Resistant Vegetable Crops, Korea. A modified version of the root-dip inoculation (RDI) method was used to inoculate the Chinese cabbage with the fungus . Briefly, three week old seedlings were removed from the soil and immersed in the conidial suspension. Samples were then collected from infected and mock-infected plants at 0 h, 3 h, 6 h, 24 h, 6 d and 11 d. Upon collection, the samples were immediately frozen in liquid nitrogen, after which they were stored at -80°C until RNA isolation. The local (fourth) and systemic (fifth) leaves were harvested as samples.
For biotic stress treatments, Chinese cabbage plants were grown in soil for six weeks under culture room conditions with a 16 h light 8 h dark cycle at 25°C temperature prior to treatment. P. carotovorum subsp. carotovorum stock (10 μl) was cultured in 25 ml of liquid YEP medium until the OD600 was 1.4, which was equivalent to 1,170,000 colony forming units (cfu) ml-1, after which it was diluted to an OD600 of 1.19, which was equal to 1 × 106 CFU ml-1, by adding double distilled water. For pathogen inoculation, 10 μl of P. carotovorum subsp. carotovorum culture solution (1 × 106 cfu ml-1) was added to a freshly needle wounded site at the lower 1/3 of the midrib of the upper third leaves, after which the sample was covered with polyvinyl bags to maintain 80-90% humidity and incubated at 25°C. All inoculations were performed three times, and the infection was confirmed by observing disease lesions in the leaves of the Chinese cabbage plants. Approximately one-third of the top of the infected leaves were harvested for RNA extraction at 0 h, 6 h, 12 h, 24 h, and 72 h after inoculation, frozen in liquid nitrogen, and stored at -80°C.
For abiotic stress treatments, Chinese cabbage seeds were grown aseptically on half-strength MS (HMS) agar medium (without sucrose) with pH of 5.8 in a culture room under a 16 h light photoperiod at 25°C. After three weeks of growth, the seedlings were transferred to fresh liquid HMS medium containing 100 μM of ABA for ABA stress treatments. To induce cold stress, the seedlings were maintained at 4°C in fresh liquid HMS medium for 24 h. In addition, pH stress treatment was applied by transferring the seedlings into fresh liquid HMS medium with pH levels of 5, 7 or 9 for 24 h. For each stress, the leaf samples were collected after 0 h, 8 h and 24 h of treatment and each sample was collected three times from two plants every time. Mock treated samples (control) were collected from the seedlings after transferring into same fresh liquid HMS medium with pH levels of 5.8 during same time interval. The collected samples were frozen immediately in liquid nitrogen and stored at -80°C until RNA isolation.
Total RNA was extracted from the roots, stems, leaves, flower buds of six different stages (early to mature flower bud) and stress treated frozen samples using an RNeasy mini kit (Qiagen, USA), after which it was treated with RNase-free DNase (Promega, USA) to remove genomic DNA contaminants.
Sequence analysis of BrLIMgenes
We recovered 22 LIM genes of B. rapa from the Brassica database (http://brassicadb.org/brad/index.php) using the key word “LIM protein”. The primary structure of genes was then analyzed using protParam (http://expasy.org/tools/protparam.html). Exons were drawn using gene structure display server 2.0 (http://gsds2.cbi.pku.edu.cn/) and open reading frames (ORFs) were determined using ORF finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html). Additionally, an alignment search was conducted using NCBI BLAST (http://www.ncbi.nlm.nih.gov/BLAST/) and the BLASTp program with the nr database. Typical domains were analyzed using the EMBL web tool (http://smart.embl.de/smart/set_mode.cgi?GENOMIC=1). Multiple protein sequences were aligned using PIR (http://pir.georgetown.edu/pirwww/search/multialn.shtml) and phylogenetic tree was constructed according to the neighbor-joining method using the MEGA5.1 software (http://www.megasoftware.net) [63, 64].
AMV one step RT-PCR kit (Takara, Japan) was used for RT-PCR. Specific primers for all genes were used for RT-PCR, and actin primers of B. rapa (FJ969844) were used as a control (Additional file 2: Table S1). We used primer 3 program for designing primers and Tm values for forward and reverse primers of each gene was determined at 55 ± 1°C. PCR was performed using 50 ng of cDNA from the roots, leaves, stems and flower buds as templates in master mixes composed of 20 pmol of each primer, 150 μM of each dNTP, 1.2 U of Taq polymerase, 1× Taq polymerase buffer, and double-distilled H2O diluted to a total volume of 20 μl in 0.5 ml PCR tubes. The samples were then subjected to the following conditions: pre-denaturing at 94°C for 5 min, followed by 30 cycles of denaturation at 94°C for 30 s, annealing at 55°C for 30 s and extension at 72°C for 1 min, with a final extension for 5 min at 72°C. Microarray expression of A. thaliana genes were created by GENEVESTIGATOR (https://www.genevestigator.com/gv/plant.jsp).
The cDNA synthesized from B. rapa sample was reverse transcribed and the trapped sequence adjacent to the vector was amplified by 3’-RACE using a TAKARA 3’-Full RACE Core set according to the manufacturer’s instructions. For BrWLIM2c and e-specific PCR the primers F-1 (5’-GTG AAC ACA ACA GAG GAG GT-3’) and F-2 (5’-AGT GAG CAC AAT CCC TCT TA-3’) were used for amplification of the target portions, respectively. The PCR products were then cloned using a TOPO TA 2.1 cloning kit (Invitrogen, USA) according to the manufacturer’s instructions and sequenced with an ABI3700 sequencer (SolGent Co., Ltd., Korea).
Real-time quantitative PCR analysis
The cDNA was synthesized using a Superscript® III First-Strand synthesis kit (Invitrogen, USA) according to the manufacturer’s instructions using oligo dT primer. Real-time quantitative PCR was performed using 1 μl of cDNA in a 25 μl reaction volume employing iTaq™ SYBR® Green PCR kit (Qiagen, California, USA). The same specific primers were used for real-time PCR and actin primers of B. rapa (FJ969844) were used as a control and the conditions for real-time PCR were as follows: 10 min at 95°C, followed by 40 cycles of 94°C for 30 s (denaturation), 58°C for 30 s (extension step), and 72°C for 45 s (signal aquisition). The fluorescence was measured following the last step of each cycle, and three replicates were used for each sample. Amplification, detection, and data analysis were conducted using a Rotor-Gene 6000 real-time rotary analyzer (Corbett Life Science, Australia). Threshold cycle (Ct) represents the number of cycles at which the fluorescence intensity was significantly higher than the background at the initial exponential phase of PCR amplification. The Br-Actin was used as the internal reference in all analyses and relative gene expression level was calculated on the basis of the 2-ΔΔCt method . Finally, the relative gene expressions in treated samples were presented as fold change compared to mock treated samples and at each time point gene expression values in mock treated samples are set to 1.
- LIM :
Lin-11, Isl-1 and Mec-3 domains
- BrLIM :
Brassica rapa LIM
Cysteine rich proteins
Actin binding proteins
Colony forming units
Half-strength Murashige and Skoog
Ubiquitin interacting motif
Rapid amplification of cDNA ends.
Dawid IB, Toyama R, Taira M: LIM domain proteins. CR Acad Sci [III]. 1995, 318: 295-306.
Taira M, Evrard JL, Steinmetz A, Dawid IB: Classification of LIM proteins. Trends Genet. 1995, 11: 431-432. 10.1016/S0168-9525(00)89139-8.
Schmeichel KL, Beckerle MC: The LIM domain is a modular protein-binding interface. Cell. 1994, 79: 211-219. 10.1016/0092-8674(94)90191-0.
Jurata LW, Pfaff SL, Gill GN: The nuclear LIM domain interactor NLI mediates homo- and heterodimerization of LIM domain transcription factors. J Biol Chem. 1998, 273: 3152-3157. 10.1074/jbc.273.6.3152.
Freydt G, Kim SK, Horvitz HR: Novel cysteine-rich motif and homeodomain in the product of the Caenorhabditis elegans cell lineage gene lin-11. Nature. 1990, 344: 876-879. 10.1038/344876a0.
Karlsson O, Thor S, Norberg T, Ohlsson H, Edlund T: Insulin gene enhancer binding protein Isl-1 is a member of a novel class of proteins containing both a homeo- and a Cys-His domain. Nature. 1990, 344: 879-882. 10.1038/344879a0.
Way JC, Chalfie M: mec-3, a homeobox-containing gene that specifies differentiation of the touch receptor neurons in C. elegans. Cell. 1988, 54: 5-16. 10.1016/0092-8674(88)90174-2.
Wadman IA, Osada H, Grutz GG, Agulnick AD, Westphal H, Forster A, Rabbitts TH: The LIM-only protein Lmo2 is a bridging molecule assembling an erythroid, DNA binding complex which includes the TAL1, E47, GATA-1 and Ldb1/NLI proteins. EMBO J. 1997, 16: 3145-3157. 10.1093/emboj/16.11.3145.
Arber S, Barbayannis FA, Hanser H, Schneider C, Stanyon CA, Bernard O, Caroni P: Regulation of actin dynamics through phosphorylation of cofilin by LIM-kinase. Nature. 1998, 393: 805-809. 10.1038/31729.
Yang N, Higuchi O, Ohashi K, Nagata K, Wada A, Kangawa K, Nishida E, Mizuno K: Cofilin phosphorylation by LIM-kinase 1 and its role in Rac-mediated actin reorganization. Nature. 1998, 393: 809-812. 10.1038/31735.
Li Y, Zheng L, Corke F, Smith C, Bevan MW: Control of final seed and organ size by the DA1 gene family in Arabidopsis thaliana. Genes Dev. 2008, 22: 1331-1336. 10.1101/gad.463608.
Xia T, Li N, Dumenil J, Li J, Kamenski A, Bevan MW, Gao F, Li Y: The Ubiquitin Receptor DA1 Interacts with the E3 Ubiquitin Ligase DA2 to Regulate Seed and Organ Size in Arabidopsis. Plant Cell. 2013, 25 (9): 3347-3359. 10.1105/tpc.113.115063.
Thomas C, Dieterle M, Gatti S, Hoffmann C, Moreau F, Papuga J, Steinmetz A: Actin bundling via LIM domains. Plant Signal Behav. 2008, 3: 320-321. 10.4161/psb.3.5.5310.
Thordal-Christensen H: Fresh insights into processes of nonhost resistance. Curr Opin Plant Biol. 2003, 6: 351-357. 10.1016/S1369-5266(03)00063-3.
Takagi S: Roles for actin filaments in chloroplast motility and anchoring. Actin: A Dynamic Framework for Multiple Plant Cell Functions. Edited by: Staiger C. 2000, Dordrecht, The Netherlands: Kluwer Academic Publishers, 203-212.
Blancaflor EB: The cytoskeleton and gravitropism in higher plants. J Plant Growth Regul. 2002, 21: 120-136. 10.1007/s003440010041.
Foissner I, Lichtscheidl IK, Wasteneys GO: Actin-based vesicle dynamics and exocytosis during wound wall formation in characean internodal cells. Cell Motil Cytoskeleton. 1996, 35: 35-48. 10.1002/(SICI)1097-0169(1996)35:1<35::AID-CM3>3.0.CO;2-H.
Hush JM, Overall RL: Cortical microtubule reorientation in higher plants: dynamics and regulation. J Microsc. 1996, 181: 129-139. 10.1046/j.1365-2818.1996.106386.x.
Shimmen T, Yokota E: Cytoplasmic streaming in plants. Curr Opin Cell Biol. 2004, 16: 68-72. 10.1016/j.ceb.2003.11.009.
Vidali L, Yokota E, Cheung AY, Shimmen T, Hepler PK: The 135 kDa actin-binding protein from Lilium longiflorum pollen is the plant homologue of villin. Protoplasma. 1999, 209: 283-291. 10.1007/BF01453456.
Kovar DR, Staiger CJ, Weaver EA, McCurdy DW: AtFim1 is an actin filament crosslinking protein from Arabidopsis thaliana. Plant J. 2000, 24: 625-636. 10.1046/j.1365-313x.2000.00907.x.
Hepler PK, Vidali L, Cheung AY: Polarized cell growth in higher plants. Annu Rev Cell Dev Biol. 2001, 17: 159-187. 10.1146/annurev.cellbio.17.1.159.
Shimmen T: The sliding theory of cytoplasmic streaming: Fifty years of progress. J Plant Res. 2007, 120: 31-43. 10.1007/s10265-006-0061-0.
Winder SJ, Ayscough KR: Actin-binding proteins. J Cell Sci. 2005, 118: 651-654. 10.1242/jcs.01670.
Grubinger M, Gimona M: CRP2 is an autonomous actinbinding protein. FEBS Lett. 2004, 557: 88-92. 10.1016/S0014-5793(03)01451-0.
Tran TC, Singleton C, Fraley TS, Greenwood JA: Cysteine-rich protein 1 (CRP1) regulates actin filament bundling. BMC Cell Biol. 2005, 6: 45-10.1186/1471-2121-6-45.
Thomas C, Hoffmann C, Dieterle M, Van Troys M, Ampe C, Steinmetz A: Tobacco WLIM1 is a novel F-actin binding protein involved in actin cytoskeleton remodeling. Plant Cell. 2006, 18: 2194-2206. 10.1105/tpc.106.040956.
Wang HJ, Wan AR, Jauh GY: An actin-binding protein, LlLIM1, mediates calcium and hydrogen regulation of actin dynamics in pollen tubes. Plant Physiol. 2008, 147: 1619-1636. 10.1104/pp.108.118604.
Papuga J, Hoffman C, Dieterle M, Moes D, Moreau F, Tholl S, Steinmetz A, Thomas C: Arabidopsis LIM proteins: a family of actin bundlers with distinct expression patterns and modes of regulation. Plant Cell. 2010, 22: 3034-3052. 10.1105/tpc.110.075960.
Peng Y, Ma W, Chen L, Yang L, Li S, Zhao H, Zhao Y, Jin W, Li N, Bevan MW, Li X, Tong Y, Li Y: Control of Root Meristem Size by DA1-RELATED PROTEIN2 in Arabidopsis. Plant Physiol. 2013, 161: 1542-1556. 10.1104/pp.112.210237.
Yang H, Shi Y, Liu J, Guo L, Zhang X, Yang S: A mutant CHS3 protein with TIR- NB-LRR-LIM domains modulates growth, cell death and freezing tolerance in a temperature-dependent manner in Arabidopsis. Plant J. 2010, 63: 283-296. 10.1111/j.1365-313X.2010.04241.x.
Lin X, Kaul S, Rounsley S, Shea TP, Benito MI, Town CD, Fujii CY, Mason T, Bowman CL, Barnstead M: Sequence and analysis of chromosome 2 of the plant Arabidopsis thaliana. Nature. 1999, 402 (6763): 761-768. 10.1038/45471.
Cheung AY, Duan QH, Costa SS, de Graaf BH, Di Stilio VS, Feijo J, Wu HM: The dynamic pollen tube cytoskeleton: live cell studies using actin-binding and microtubule-binding reporter proteins. Mol Plant. 2008, 1 (4): 686-702. 10.1093/mp/ssn026.
Theologis A, Ecker JR, Palm CJ, Federspiel NA, Kaul S, White O, Alonso J, Altafi H, Araujo R, Bowman CL, Brooks SY, Buehler E, Chan A, Chao Q, Chen H: Sequence and analysis of chromosome 1 of the plant Arabidopsis thaliana. Nature. 2000, 408 (6814): 816-820. 10.1038/35048500.
Tabata S, Kaneko T, Nakamura Y, Kotani H, Kato T, Asamizu E, Miyajima N, Sasamoto S: Sequence and analysis of chromosome 5 of the plant Arabidopsis thaliana. Nature. 2000, 408 (6814): 823-826. 10.1038/35048507.
Mayer K, Schuller C, Wambutt R, Murphy G, Volckaert G, Pohl T, Dusterhoft A, Stiekema W, Entian KD: Sequence and analysis of chromosome 4 of the plant Arabidopsis thaliana. Nature. 1999, 402 (6763): 769-777. 10.1038/47134.
Eliasson A, Gass N, Mundel C, Baltz R, Kräuter R, Evrard JL, Steinmetz A: Molecular and expression analysis of a LIM protein gene family from flowering plants. Mol Gen Genet. 2000, 264: 257-267. 10.1007/s004380000312.
Arnaud D, Dejardin A, Leple JC, Lesage-Descauses MC, Pilate G: Genome-wide analysis of LIM gene family in Populus trichocarpa, Arabidopsis thaliana, and Oryza sativa. DNA Res. 2007, 14: 103-116. 10.1093/dnares/dsm013.
Baltz R, Schmit AC, Kohnen M, Hentges F, Steinmetz A: Differential localization of the LIM domain protein PLIM-1 in microspores and mature pollen grains from sunflower. Sex Plant Reprod. 1999, 12: 60-65. 10.1007/s004970050172.
Sweetman J, Spurr C, Eliasson A, Gass N, Steinmetz A, Twell D: Isolation and characterization of two pollen-specific LIM domain protein cDNAs from Nicotiana tabacum. Sex Plant Reprod. 2000, 12: 339-345. 10.1007/PL00009836.
Alves-Ferreira M, Wellmer F, Banhara A, Kumar V, Riechmann JL, Meyerowitz EM: Global expression profiling applied to the analysis of Arabidopsis stamen development. Plant Physiol. 2007, 145: 747-762. 10.1104/pp.107.104422.
Ye JR, Zhou LM, Xu ML: Arabidopsis LIM proteins PLIM2a and PLIM2b regulate actin configuration during pollen tube growth. Biol Plant. 2013, 57 (3): 433-441. 10.1007/s10535-013-0323-3.
Aist JR: Papillae and related wound plugs of plant cells. Annu Rev Phytopathol. 1976, 14: 145-163. 10.1146/annurev.py.14.090176.001045.
Tomiyama K, Sato K, Doke N: Effect of cytochalasin B and colchicines on hypersensitive death of potato cells infected by incompatible race of Phytophthora infestans. Ann Phytopathol Soc Jpn. 1982, 48: 228-230. 10.3186/jjphytopath.48.228.
Hazen BE, Bushnell WR: Inhibition of the hypersensitive reaction in barley to powdery mildew by heat shock and cytochalasin B. Physiol Plant Pathol. 1983, 23: 421-438. 10.1016/0048-4059(83)90026-7.
Takemoto D, Jones DA, Hardham AR: GFP-tagging of cell components reveals the dynamics of subcellular re-organization in response to infection of Arabidopsis by oomycete pathogens. Plant J. 2003, 33: 775-777. 10.1046/j.1365-313X.2003.01673.x.
Kobayashi I, Kobayashi Y, Hardham AR: Dynamic reorganization of microtubules and microfilaments in flax cells during the resistance response to flax rust infection. Planta. 1994, 195: 237-247.
Kobayashi I, Kobayashi Y, Yamaoka N, Kunoh H: Recognition of a pathogen and a nonpathogen by barley coleoptile cells. III. Responses of microtubules and actin filaments in barley coleoptile cells to penetration attempts. Can J Bot. 1992, 70: 1815-1823. 10.1139/b92-225.
Collins NC, Thordal-Christensen H, Lipka V, Bau S, Kombrink E, Qiu JL, Huckelhoven R, Stein M, Freialdenhoven A, Somerville SC, Schulze-Lefert P: SNARE-protein-mediated disease resistance at the plant cell wall. Nature. 2003, 425: 973-977. 10.1038/nature02076.
Vierstra RD: The ubiquitin-26S proteasome system at the nexus of plant biology. Nat Rev Mol Cell Biol. 2009, 10: 385-397. 10.1038/nrm2688.
Dreher K, Callis J: Ubiquitin, hormones and biotic stress in plants. Ann Bot. 2007, 99: 787-822. 10.1093/aob/mcl255.
Trujillo M, Shirasu K: Ubiquitination in plant immunity. Curr Opin Plant Biol. 2010, 13: 402-408. 10.1016/j.pbi.2010.04.002.
Liu YC, Wu YR, Huang XH, Sun J, Xie Q: AtPUB19, a U-box E3 ubiquitin ligase, negatively regulates abscisic acid and drought responses in Arabidopsis thaliana. Mol Plant. 2011, 4 (6): 938-946. 10.1093/mp/ssr030.
Ikeda F, Dikic I: A typical ubiquitin chains: new molecular signals, ‘Protein modifications: Beyond the Usual Suspects’ review series. EMBO Rep. 2008, 9: 536-542. 10.1038/embor.2008.93.
Smertenko A, Franklin-Tong VE: Organisation and regulation of the cytoskeleton in plant programmed cell death. Cell Death Differ. 2011, 18: 1263-1270. 10.1038/cdd.2011.39.
Higaki T, Kurusu T, Hasezawa S, Kuchitsu K: Dynamic intracellular reorganization of cytoskeletons and the vacuole in defense responses and hypersensitive cell death in plants. J Plant Res. 2011, 124: 315-324. 10.1007/s10265-011-0408-z.
Day B, Henty JL, Porter KJ, Staiger CJ: The pathogen-actin connection: A platform for defense signaling in plants. Annu Rev Phytopath. 2011, 49: 489-506.
Hardham AR, Takemoto D, White RG: Rapid and dynamic subcellular reorganization following mechanical stimulation of Arabidopsis epidermal cells mimics responses to fungal and oomycete attack. BMC Plant Biol. 2008, 8: 63-10.1186/1471-2229-8-63.
Arnaud D, Dejardin A, Leple JC, Lesage-Descauses MC, Boizot N, Villar M, Benedetti H, Pilate G: Expression analysis of LIM gene family in poplar, toward an updated phylogenetic classification. BMC Res Notes. 2012, 5: 102-10.1186/1756-0500-5-102.
Zhao M, He L, Gu Y, Wang Y, Chen Q, He C: Genome-wide analyses of a plant-specific LIM-domain gene family implicate its evolutionary role in plant diversification. Genome Biol Evol. 2014, 6 (4): 1000-1012. 10.1093/gbe/evu076.
Bi D, Johnson KCM, Zhu Z, Huang Y, Chen F, Zhang Y, Li X: Mutations in an atypical TIR-NB-LRR-LIM resistance protein confer autoimmunity. Front Plant Sci. 2011, 2: 71-
Ospina-Giraldo MD, Mullins E, Kang S: Loss of function of the Fusarium oxysporum SNF1 gene reduces virulence on cabbage and Arabidopsis. Curr Genet. 2003, 44: 49-57. 10.1007/s00294-003-0419-y.
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 5.1. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
This research was supported by basic Science Research Program though the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF-2010-355-C00054, 2012R1A1A2044500). and Golden Seed Project (Center for Horticultural Seed Development, No. 213003-04-2-SB420), Ministry of Agriculture, Food and Rural Affairs (MAFRA), Ministry of Oceans and Fisheries (MOF), Rural Development Administration (RDA) and Korea Forest Service (KFS).
The authors declare that they have no competing interests.
Conceived and designed the experiments: JIP, NUA, MW, ISN. Performed the experiments: NUA, JIP, HJJ. Analyzed the data/gave suggestions: NUA, JIP, HJJ, SKTA, MYC, YGC, MW, ISN. Contributed reagents/materials/analysis tools: JIP, NUA, HJJ, SKTA, MYC, YGC, MW, ISN. Wrote the paper: NUA and JIP. All authors read and approved the final manuscript.
Jong-In Park, Nasar Uddin Ahmed contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1: Gene structure of 22 LIM genes of Brassica rapa. Solid boxes and lines indicate exons and introns respectively. Length of exons and introns are mentioned below in base pairs (bps). Figure S2. Schematic representation of the existing and the postulated genomic structures of A) BrLIM7 and B) BrLIM9. Solid boxes and lines indicate exons and introns respectively and their lengths are mentioned above in base pairs (bps). Specific primers were indicated by arrows, forward (→) and reverse (←), with corresponding numbers for RT-PCR expression analysis over different organs, roots (R), stems (S), leaves (L) and flower buds (Fb), and genomic DNA amplifications, which are represented beneath each genes. Primer #1 and #3, indicated by green arrows (→), were used for 3’-RACE PCR of BrLIM7 and BrLIM9, respectively. (PDF 449 KB)
Additional file 2: Table S1: Primers specific for 22 BrLIM genes used for RT and Real time PCR analysis. Table S2. Sequence relatedness among 22 LIM proteins of Brassica rapa.(PDF 80 KB)
Additional file 3: Table S3: List of genes selected for phylogenetic analysis with sequences of deduced proteins of B. rapa.(PDF 164 KB)
Additional file 4: Figure S3: Microarray expressions of six LIM, one DA1 and seven DAR genes of Arabidopsis thaliana after infection with Fusarium oxysporum (created by GENEVESTIGATOR). Figure S4. Microarray expressions of six LIM, one DA1 and seven DAR genes of Arabidopsis thaliana after cold stress treatments (created by GENEVESTIGATOR). (XLS 44 KB)
Authors’ original submitted files for images
About this article
Cite this article
Park, JI., Ahmed, N.U., Jung, HJ. et al. Identification and characterization of LIM gene family in Brassica rapa. BMC Genomics 15, 641 (2014). https://doi.org/10.1186/1471-2164-15-641
- Brassica rapa
- Expression analysis
- LIM gene
- Stress resistance