BACKGROUND
the basic/helix-loop-helix  proteins have dna-binding and dimerization capabilities.
they are a superfamily of transcription factor  that have been found to have many different functions in essential physiological and developmental process in animals and plants  <cit> .
the bhlh domain contains approximately  <dig> amino acids with two functionally distinct regions, the basic region and the hlh regions  <cit> .
the basic region was  <dig> amino acids long and typically included six basic residues.
it was located at the n terminus of the domain and functions as a dna binding motif  <cit> .
the hlh region contains two amphipathic α helices separated by a loop region of variable length.
the hlh region acts as a dimerization domain and allows the formation of homodimers or heterodimers  <cit> .
among all bhlh motifs,  <dig> amino acids have been found to be highly conserved in organisms ranging from yeast to mammals  <cit> .
outside of these conserved bhlh domains, the proteins exhibited considerable sequence divergence.
some bhlh proteins have been shown to bind to the sequences containing the core element known as the e box , with the most common form of g-box .
the nucleotides flanking the core element may also have a role in binding specificity  <cit> .
based on the phylogenetic relationships, dna-binding motifs, and functional properties, the bhlh tfs family has been divided into six main groups in metazoans  <cit> .
in brief, group a bhlh proteins can bind to the cagctg core sequences of e-boxes.
group b includes a large number of functionally proteins  and bind to the g-box sequence cacgtg  <cit> .
group c contains an additional protein-protein interaction region  and binds to acgtg or gcgtg sequences.
group d proteins have the hlh region but lack the basic dna binding domain  <cit> .
group e proteins have pro or gly residues within the basic region and can bind preferentially to a typical sequence, cacgng  <cit> .
group f consists of the coe domain; they have diverse sequences compared with other groups and another domain for dimerization and dna binding  <cit> .
only a small number of plant bhlh proteins have been characterized functionally, far fewer than in animals.
in arabidopsis,  <dig> bhlh-encoding genes which were divided into  <dig> subfamilies according to their phylogenetic relationships have been identified from the analysis of genome sequences  <cit> .
a total of  <dig> and  <dig> bhlh tfs have been identified in the rice  and chinese cabbage  genomes, respectively  <cit> .
these have been divided into  <dig> and  <dig> subfamilies, respectively.
phylogenetic analysis showed that the plant bhlh proteins comprised  <dig> subfamilies,  <dig> of which were present in the common ancestors of extant mosses and vascular plants  <cit> .
most bhlh proteins identified so far have been functionally characterized in arabidopsis, and their roles have been shown to include regulation of fruit dehiscence, anther and epidermal cell development, hormone signaling, and stress responses  <cit> .
tomato  is an economically important vegetable worldwide.
the annual global production of tomato in  <dig> was more than  <dig> million tons including  <dig> million tons in china .
the tomato genome has been sequenced and assembled by the international tomato genome sequencing project , because tomatoes are economically important and it is model species for the study of fruit ripening  <cit> .
a high-quality genome sequence for domesticated tomato and more than  <dig>  proteins have been obtained.
tomato yellow leaf curl virus  is the most widespread and currently ranks 3rd among the most economically and scientifically most important plant viruses worldwide  <cit> .
the symptoms of tylcv infection in young plants include stunted growth, upward curling of leaf margins, marked reduction in leaf size, mottling and yellowing of young leaves, and flower abscission, leading to severe yield loss  <cit> .
currently five major loci resistant to tylcv have been identified from different wild tomato relatives, ty- <dig>  ty- <dig> and ty- <dig> from s.chilense, ty- <dig> from s.habrochaites, and ty- <dig> from s.peruvianum  <cit> .
among them, ty- <dig> and ty- <dig> were found to be allelic and have been cloned.
ty- <dig> and ty- <dig> were found to be allelic and have been cloned.
they are rna-dependent rna polymerases  and may be involved in rna silencing  <cit> .
in addition, ty- <dig>  ty- <dig> and ty- <dig> have been mapped to chromosomes  <dig>   <dig>  and  <dig> respectively, using molecular markers  <cit> .
cdna library comparisons of susceptible and resistant tomato lines before and after tylcv infection showed approximately  <dig> genes that are preferentially expressed in a tomato line with a resistance introgressed from s. habrochaites  <cit> .
using whole transcriptome sequencing of the tylcv-resistant tomato breeding line cln2777a  and tylcv-susceptible tomato breeding line tmxa48-4- <dig> ,  <dig> and  <dig> genes were found to be differentially expressed in the r and s tomato lines, respectively  <cit> .
in tomatoes, lefer, a bhlh protein encoded by solyc06g <dig> .
<dig>  slybhlh <dig>  was the first identified regulator of iron nutrition in plants.
lefer plays an important role in the fe-deficiency response of tomatoes  <cit> .
style <dig> , encoded by solyc02g <dig> .
<dig>  slybhlh <dig>  is the major quantitative trait locus responsible for style length; this important floral attribute has been shown to be associated with the evolution of self-pollination and was cloned in cultivated tomatoes  <cit> .
however, the tomato bhlh protein family has not been analyzed at a genome-wide level, and the phylogenetic relationship of this protein family remains poorly understood.
in this study, a total of  <dig> slybhlh genes were identified in the tomato genomic sequence and phylogenetic analyses were carried out to evaluate the relationships among these genes.
changes in global expression pattern of slybhlh genes in r and s lines infected by tylcv were analyzed to provide insight into the regulation of response to tylcv.
the expression of slybhlh exhibited a variety of expression patterns, suggesting a novel layer of regulation for the response to tylcv in tomato.
methods
database search for bhlh genes
the pfam database   <cit>  was used to screen the genome of tomato  and potato .
proteins with helix-loop-helix dna-binding domains  were used to identify the putative bhlh proteins in tomato and potato using the hidden markov model .
the hmmsearch tool, with an expected value  cut-off of  <dig>  was used to identify the proteins.
these sequences were then verified using the smart tool   <cit> .
the arabidopsis thaliana bhlh proteins were retrieved from the tair database  using a previous report  <cit> .
phylogenetic analysis and identification of conserved motifs and gene structure
the complete amino acid sequences were screened against the pfam database to identify the domains of bhlh transcription factors.
mega <dig> software was used to construct neighbor-joining  distance trees using tomato bhlh protein domain sequences  <cit> .
the bootstrap was set as  <dig>  replicates, which provided information regarding their statistical reliability.
meanwhile, the nj method of the phylip software  was also used with bootstrap of  <dig> replicates to create another phylogenetic tree to validate the results from the nj method by mega  <dig> software.
a phylogenetic tree of all the identified bhlh protein domains was also constructed.
the identified bhlh domains were aligned using a clustalx  <dig>  program with default settings  <cit> .
to identify the conserved motifs in tomato bhlh proteins, the multiple expectation-maximization for motif elicitation  program version  <dig> .
<dig>  <cit>  was used with default parameters, except for the following parameters:  optimum motif width was set to ≥ <dig> and ≤100;  the maximum number of motifs was set to identify ten motifs.
meme software  was used to search for conserved motifs in the complete amino acid sequences of bhlh proteins.
the coding domain sequences  and dna sequences of tomato bhlh genes were used to assess gene structure using gsds   <cit> .
collinear correlations of bhlh genes in the tomato, potato, and arabidopsis genomes
orthomcl program   <cit>  was used to identify the orthologous and paralogous genes in tomatoes, potatoes and arabidopsis.
briefly, the tools blastp, with an e-value ≤ 1e− <dig>  and orthomclpairs were used to find orthologs, inparalogs and coorthologs in these three species.
the circos tool was used to link these genes to chromosomes  <cit> .
in addition, the relationships of orthologous and paralogous genes in these three species were also shown using the circos tool  <cit> .
the bhlh genes in tomato were searched for duplication events .
chromosome distribution and gene duplications
to determine the physical locations of bhlh genes, the starting and ending positions of all bhlh genes on each chromosome were obtained from the tomato database.
the mapinspect software was used to draw the images of the locations of the tomato bhlh genes .
we used the plant genome duplication database  to retrieve the duplicate chromosomal blocks and then identify the bhlh genes in the duplication block which allowed us to identify duplicate tomato bhlh genes  <cit> .
the pgdd is a public database used to identify and catalogue plant genes in terms of intra-genomic or cross-genomic syntenic relationships.
rna data collection and data mining
transcriptomic data of tylcv-resistant breeding line, cln2777a  and susceptible breeding line, tmxa48-4- <dig>  with uninfected  and mixed infection samples of  <dig>   <dig>  and 7 days post infection  were downloaded from ncbi sra database  and analyzed as described in a previous study  <cit> .
enrichment of gene ontology  categories was performed with an agrigo analysis toolkit   <cit>  using the topgo ‘elim’ algorithm  <cit>  for the aspects ‘biological process’ and ‘subcellular localization’.
the selected categories were sorted from the lowest to the highest p value .
validation of differentially expressed genes by quantitative rt-pcr
four bhlh genes with differentially expressed in r or s lines were selected and subjected to quantitative rt-pcr validation .
primers for quantitative rt-pcr were designed using primer <dig> software and primer specificity was evaluated by blasting primer sequences against the ncbi database.
pcr amplifications were performed in a real-time thermal cycler qtower  <dig> / <dig>   with 15 μl of final volumes containing  <dig>  μl of cdna,  <dig>  μl each primer , 6 μl of sterile water, and  <dig>  μl  sybr premix extaq™ ii kit .
the conditions for amplification were as follows: 5 min of denaturation at 95°c followed by 40 cycles of 95°c for 10 s, 60°c for 20 s, and 72°c for 10 s. the expression levels of selected genes were normalized to α-tubulin  expression  <cit> .
relative gene expression was calculated using the 2-ΔΔct method  <cit> .
three biological replicates were performed for each of the selected genes.
validation of candidate genes with virus-induced gene silencing  and cell death analysis
the tobacco rattle virus  mediated vigs system was used to silence a bhlh gene   <cit> .
briefly, ptrv-containing agrobacterium eha <dig> was cultured in liquid lb medium and resuspended in infiltration medium at an o.d.
value of  <dig>  and cultured at room temperature for 4 h. three week old seedlings were infiltrated by pressure inoculation in the leaves with a needleless syringe.
for the vigs experiments, agro infiltration was performed two weeks after trv inoculation.
death cells were identified by staining with lacto-phenol trypan blue and dab as previous described  <cit> .
to visualize cell death, stem were stained by boiling in lacto-phenol trypan blue , followed by destaining with chloral hydrate .
then the death cell was examined with ni-u microscope .
RESULTS
identification and classification of bhlh proteins in tomato
to identify the putative bhlh proteins in the tomato genome, a hidden markov model search resulted in the identification of  <dig> bhlh proteins .
to verify the reliability of our criteria, we performed simple modular architecture research tool  analysis of  <dig> putative slybhlh protein sequences and found all of them had a typical bhlh domain.
the number of bhlh tfs in tomatoes exceeded that of many metazoans and fungi, but was less than that found in some plants, such as rice , chinese cabbage , potatoes , soybeans  and maize   <cit> .
in addition, the density of bhlh proteins in the entire tomato  and potato  genomes was found to be less than that in most plant species, such as arabidopsis  and rice   <cit> .
most angiosperm plant lineages have experienced one of more rounds of ancient polyploidy  <cit> .
and the genomes of tomato and potato have undergone recent triplication events, whereas few individual tomato/potato genes remain triplicated  <cit> .
therefore, this event might lead to relatively fewer bhlh genes in tomato compared to other plant species.
multiple sequence alignments, predicted dna-binding ability and conserved residues
to examine sequence features of these tomato bhlh domains, multiple sequence alignment of the  <dig> bhlh amino acid sequences were performed.
there were four conserved regions in the bhlh domain sequences, including one basic region, two helix regions and one loop region .
the basic regions have five basic residues, but five of these proteins did not have the basic region .
the loop was found to be the most divergent region in terms of both length and amino acid composition.
from the alignment,  <dig> residues were identified that were identical in at least 50% of the  <dig> tomato bhlh domains .
among these  <dig> residues, nine residues were present in more than 75% sequences .figure  <dig>
the characterization and distribution of the bhlh domains.
a. sequence logo of the slybhlh domain by meme.
the h <dig>  e <dig> and r <dig> amino acids in the basic domain that are important for dna binding are indicated by stars.
amino acids important for dimerization of the helix-loop-helix domain are indicated by arrows.
b. distribution of amino acids in the bhlh consensus motif among tomatoes.
the numbers of horizontal ordinate refer to the positions of the residues in the alignments of the studies.
c. predicted dna-binding characteristics of the bhlh domain of slybhlh proteins.
five residues , five residues , two residues  and seven residues  made up the basic region, the first helix region, the loop region and the second helix region, respectively.
all of these conserved residues were consistent with previous studies  <cit> .
the leu- <dig> in the basic region was conserved in all  <dig> bhlh proteins, suggesting that this residue is extremely important for promoting the formation of dimerization among bhlh proteins  <cit> .
the basic region of the bhlh domain can bind to dna and is critical for function  <cit> .
using the criteria described by massari and murre, the slybhlh proteins were divided into several categories based on sequence information in the n-terminal region of the bhlh domains   <cit> .
as was done with arabidopsis and chinese cabbage, the slybhlh proteins of tomato were also divided into two major groups according to 17 n-terminal amino acids within bhlh protein domain, including  <dig> dna-binding and  <dig> non-dna binding proteins.
the dna-binding bhlhs were further divided into two groups with different predicted target sequences depending on the presence or absence of residues glu- <dig> and arg- <dig> in the basic region.
group  proteins had  <dig> putative e-box-binding proteins with conserved glu-9/arg- <dig> residues and group  proteins had  <dig> non-e-box-binding proteins lacking these residues .
the three residues in the basic region of the bhlh domain, his/lys- <dig>  glu- <dig> and arg- <dig>  were found to constitute the classic g-box-binding region  <cit> .
group  can therefore be subdivided further into two subgroups: 1a <dig>  whose  <dig> proteins are predicted to bind g-boxes, and 1a <dig>  whose three members are predicted to bind other types of e boxes .
phylogenetic analysis of the bhlh transcription factor family
to assess the evolutionary relationships of the slybhlh genes, an nj phylogenetic tree was generated using the multiple sequences alignments of the conserved bhlh tf domains in tomato and arabidopsis with a bootstrap value of  <dig> .
twenty-six subfamilies were identified according to the clades support values, topology of the tree, and classification of the arabidopsis  <cit> .
no slybhlh proteins in the xiii, ii, subfamilies relative to those of arabidopsis, therefore, tomato contained  <dig> bhlh subfamilies in our analysis .
to further validate the reliability of the nj tree with mega  <dig> , nj and maximum parsimony analysis was also used to generate phylogenetic trees using phylip software .
<dig> %  of the slybhlh proteins with nj model using phylip software were placed into the same subfamilies as those in the nj tree with mega  <dig> , indicating that both methods are in very good agreement.figure  <dig>
phylogenetic tree constructed from the neighbor-joining method using the bhlh transcription factor domain in tomato.
the phylogenetic tree was constructed using mega <dig> software.
the numbers are bootstrap values based on  <dig> iterations.
only bootstrap values larger than 50% support are indicated.
in order to assess differences in protein structure, meme was used to identify conserved motifs in the tomato bhlh proteins.
ten conserved motifs were identified and named motif  <dig> through motif  <dig>  in general, the bhlh proteins were clustered in the same subfamilies and shared similar motif compositions, which indicated functional similarities among members of the same subfamilies  <cit> .
the tomato bhlh proteins were found to have a similar structure for every subfamily .
the pattern of intron position can also provide important evidence to support phylogenetic relationships in a gene family.
here, gsds tools were used to show the gene structures for slybhlh genes .
among the  <dig> tomato bhlh genes, the number of introns ranged from  <dig> to  <dig>  and most members of the same subfamilies had similar intron/exon structures.
for example, the members of subfamilies iii and xv each have only one intron.
these results demonstrated that proteins within the same subfamily share close evolutionary relationships.
collinear correlations of bhlh genes in tomatoes, potatoes, and arabidopsis
the solanum lineage has experienced two consecutive genome triplications: one is ancient and shared with rosids, and the other more recent  <cit> .
in this study, the correlation between tomato, potato, and arabidopsis bhlh genes was analyzed using the orthomcl program.
here,  <dig> gene pairs were found to be orthologous between tomatoes and potatoes, but only  <dig> orthologous pairs were found between tomato and arabidopsis .
these results were consistent with the close relationship of tomatoes and potatoes.
among the orthologous gene pairs shared by tomatoes and potatoes, each tomato bhlh gene had one to four potato bhlh genes.
these results demonstrated that bhlh tf genes in potato were duplicated accompanied with evolution processes.
in addition, paralogous bhlh gene pairs were also analyzed.
a total of  <dig>   <dig>  and  <dig> bhlh gene pairs were identified in arabidopsis, tomatoes, and potatoes, respectively .
visualization of the relationships of paralogous and paralogous bhlh genes among these three species was performed using the circos software .figure  <dig>
comparative analysis of synteny and expansion of bhlh genes.
a. twelve tomato chromosomes  and  <dig> of potato chromosomes  maps were based on the orthologous and paralogous pair positions and demonstrate highly conserved synteny.
the red lines represent the orthologous bhlh genes between tomato and potato.
the blue and yellow lines represent the paralogous bhlh genes in potato and tomato, respectively.
b. twelve tomato  and arabidopsis chromosome  maps were based on the orthologous and paralogous pair positions and demonstrate highly conserved synteny.
the red lines represent the orthologous bhlh genes between tomato and arabidopsis.
the blue and yellow lines represent the paralogous bhlh genes in tomato and arabidopsis, respectively.
chromosome distribution and gene duplication of the bhlh tf family
the physical map positions of the bhlh genes on tomato chromosomes were identified .
among the  <dig> bhlh tfs,  <dig> were mapped onto the twelve tomato chromosomes except slybhlh <dig>  most bhlh tfs were found on chromosome  <dig>  and  <dig> .
in contrast, there are only  <dig>  and  <dig>  bhlh tf genes on chromosome  <dig> and  <dig>  respectively.
furthermore, bhlh genes were found to be mapped on the chromosomes with an obviously uneven distribution, and some bhlh genes gathered on part of the chromosome.
relative high densities of bhlh genes were observed in some chromosomal regions, including the bottom of chromosomes  <dig>   <dig>   <dig>  and  <dig>  for example,  <dig> genes clustered in the end of chromosome  <dig> and  <dig> with density of  <dig>  and  <dig>  genes per mb, respectively.
in contrast, several large chromosomal regions lacked bhlh genes, such as the top half of chromosomes  <dig> and  <dig> and the central section of chromosomes  <dig>   <dig>   <dig>  and  <dig> .figure  <dig>
distribution of  <dig> bhlh genes on the twelve tomato chromosomes using mapinspect software.
only slybhlh <dig> could not be anchored onto a specific chromosome.
the scale is in megabases .
the duplicate bhlh genes are connected with blue lines.
previous reports have analyzed duplication events in rice and chinese cabbage  <cit> .
in the current analysis, we first retrieved the genome chromosome blocks in the tomato with pgdd database and identified the  <dig> duplicate blocks.
a total of  <dig> duplicate bhlh genes pairs were located in these blocks .
these duplication bhlh genes are derived from the same subfamily, indicating that members of some bhlh subfamilies originated from the duplication events.
differential expression of slybhlh genes in response to tylcv infection
the rna-seq technology has been shown to provide precise digital information related on gene expression and can discriminate genes of high sequence identity  <cit> .
using this technology, global gene expression changes of the leaves of two tomato breeding lines, tylcv-resistant cln2777a and tylcv-susceptible tmxa48-4- <dig>  have been analyzed before  and after  tylcv infection with viruliferous whiteflies  <cit> .
a total of  <dig>  transcripts were detected from r and s lines by alignment to the tomato genome including the  <dig>  novel transcripts predicted in tomatoes.
the expression levels of mapped genes were normalized with a value of fragments per kilobase of exon per million fragments mapped .
if the fpkm value of gene was above zero, the gene was considered expressed.
on the basis of this criterion, the expression level of  <dig> slybhlh genes was confirmed in the s and r line with  <dig> and  <dig> dpi.
a subset of  <dig> slybhlh genes  was expressed under both conditions .
go enrichment analysis revealed that the products of most of slybhlh genes were localized in the nucleus .
in addition, some biological processes such as ‘regulation of macromolecule metabolic process’ , ‘regulation of gene expression’  and ‘regulation of metabolic process’  were overrepresented among all slybhlh genes, indicating that the bhlh genes were involved in transcription and metabolic regulation .figure  <dig>
the expression levels and related functions of bhlh genes.
a. venn’s diagram of the cross comparison of expressed genes based on rna-seq data.
b. gene ontology  enrichment analysis of expressed bhlh genes involved in the biological processes.
c. subcellular location of bhlh genes.
out of  <dig> slybhlh genes in the current study, only four were differentially expressed before and after tylcv infection  .
among the four differentially expressed slybhlh genes, slybhlh <dig>  was up-regulated in the r line and down-regulated in the s line after tylcv infection.
the expression level of four differentially expressed genes was determined using quantitative rt-pcr to validate the gene expression data from rna-seq .
the results demonstrated that all tested genes revealed a similar trend of transcript accumulation as in rna-seq analysis .table  <dig>
differentially expressed
slybhlh
genes between before and after tylcv infection
ty-2
4840
<dig> dpi
<dig> dpi
fold change
<dig> dpi
<dig> dpi
fold changequantitative rt-pcr was used to measure the relative expression levels of five pathogen resistance related genes in the r  and s  lines, with tomato α-tubulin  as an internal reference.
a. slybhlh <dig>  b. slybhlh <dig>  c. slybhlh <dig>  d. slybhlh <dig>
vigs validation of slybhlh <dig> gene
to investigate the role of tylcv as it related to resistant in tomato, the slybhlh <dig> gene was further challenged with tylcv after vigs at the cotyledon stage.
one month after agroinfiltration, the success of the trv silencing system was confirmed by the appearance of cell death in the leaves of s plantlets treated with ptrv <dig> and ptrv2- slybhlh <dig> .
we also observed the cell death development under ni-u microscope and found that the vigs lines with slybhlh <dig> gene triggered a rapid cell death response in comparison with empty vector lines by trypan blue and dab staining .
a quantitative rt-pcr analysis showed there were significantly fewer slybhlh <dig> transcripts in slybhlh131-silenced tomato leaves during trv infection , indicating that slybhlh <dig> was effectively silenced in tomatoes.figure  <dig>
validation of the bhlh gene  with virus-induced gene silencing.
a. cotyledon agroinfiltration of trv vectors was carried out in the s plantlets at the cotyledon stage.
r plants treated with phytoene desaturase  gene silencing constructs ptrv <dig> and ptrv <dig> -  showed bleached areas in leaflets .
s plants treated with ptrv <dig> and ptrv <dig> vectors showed the normal phenotype .
s plantlets treated with slybhlh <dig> gene silencing constructs ptrv <dig> and ptrv <dig> -  showed a cell death phenotype .
arrow indicates cell death.
b. trypan blue staining of vigs and empty vector lines leaves.
c. dab staining of vigs and empty vector lines leaves.
d. relative expression of slybhlh <dig> transcript using real-time rt-pcr analysis in the vigs-treated lines 20 days after agroinfiltration with trv vectors.
tomato α-tubulin  was used as an internal reference.
error bars represented se of three biological replicates and asterisks indicate significant differences by student’s t test for p <  <dig> .
discussion
the solanaceae is one of the largest angiosperm genera and includes many different vegetables consumed by humans.
in recent years, many plants in this genera have been sequenced, including tomatoes  <cit> , potatoes  <cit> , peppers  <cit> , and tobacco  <cit> .
with the rapid development in bioinformatic analyses, the information stored in various genomes may be explored to elucidate the mechanisms that regulate the development and response to biotic and abiotic stresses.
tomatoes are major crop plant and a model system for fruit development.
as of  <dig>  tomato production had doubled from  <dig> million tons in  <dig>  <cit> .
however, tomato yellow leaf curl virus disease causes huge losses in tomato production worldwide.
this disease is caused by different related be gomovirus species.
a previous study extended our basci understanding of the response of tomato to tylcv infection by comparing whole transcriptome expression changes between a tylcv-resistant line and a tylcv-susceptible line.
in present study,  <dig> bhlh transcription factor genes were identified in tomatoes and their differential expression was analyzed in r and s lines before and after tylcv infection.
four differentially expressed genes were identified in r and s lines, in which slybhlh <dig> and slybhlh <dig> were derived from the chromosome  <dig> and  <dig>  and slybhlh <dig> and slybhlh <dig> were mapped on the chromosome  <dig>  these genes are not located in the region of five known loci , so they are not the most important genes but might be involved in the regulation network of tylcv resistance.
phylogenetic analysis of the bhlh domain allows division of the slybhlh family into  <dig> subfamilies.
the clustering of the members within these subfamilies was further supported by additional analysis with regard to other criteria, such as predicted dna binding capacity and sequence specificity, exon/intron distribution pattern with the domain.
these data support the general conclusion that members within subfamilies may have recent common evolutionary origins resulting from various genomic duplication events.
they may have related molecular functions.
the bhlh subfamily iiid transcription  function function redundantly to negatively regulate jasmonate  responses in arabidopsis  <cit> .
however, the strong sequence diversity outside of the bhlh domain across the members of the slybhlh family suggests that the expansion of this family in tomatoes involved extensive domain shuffling, as in other organisms.
however, the non-bhlh amino acid motifs are conserved in each of the bhlh subfamilies , suggesting that the conservation of these extra domains during plant evolution may have been essential to the function of the bhlh proteins in their respective subfamilies  <cit> .
the core dna binding domain of the bhlh proteins contained the basic region of the bhlh domain and these residues were found to recognize and bind to the core hexanucleotide  <cit> .
the amino acid sequence in this region provides the major subdivision of the bhlh family, dividing these proteins into those that are predicted to bind dna and those that are not .
key residues in this region confer the capacity to discriminate the variants of the hexanucleotide motif with the canonical e-box  and non-e-box motifs.
additional residues within the basic region confer further dna binding site sequence selectivity.
these include g-box and non-g-box core motifs.
the current analysis showed there are  <dig> e-box binding bhlh proteins in tomatoes, and this ration was lower than in arabidopsis and chinese cabbage, indicating that many different binding motifs exist in the tomato bhlh protein family.
most bhlh proteins identified so far were functionally characterized in arabidopsis, and their roles include plant development, fruit dehiscence, phytochrome signaling, hormone signaling and stress responses, such as cold, heat, abscisic acid, jasmonic acid, and the light signaling pathway  <cit> .
only two bhlh genes have been functionally characterized in tomato.
one is the fer gene , which is involved in the response to iron acquisition and supply  <cit> .
another is style <dig>  , that is associated with the evolution of self-pollination; its natural genetic variation in the promoter is responsible for evolution from allogamy to autogamy in the cultivated tomatoes  <cit> .
in the current rna-seq data set, the fer gene was not expressed in any of the four samples, and the style <dig>  gene was expressed but not differentially expressed in r or s lines after tylcv infection.
members of the same plant bhlh subfamily are frequently involved in the same biological processes.
usually the functions of these proteins overlap, causing them to be partially or totally redundant  <cit> .
the characterized function of bhlh in other species can help the user to predict the function of tomato bhlh in the same subfamily.
jasmonates  are lipid-derived hormones that regulate plant responses to stresses such as wounding and pathogens invasion  <cit> .
ja negatively regulates plant growth and is considered to modulate the distribution of energy to defense responses  <cit> .
in arabidopsis, the bhlh subfamily iiid members  act as transcription repressors and function redundantly to negatively regulate ja response.
in tomatoes, eight members of the iii subfamily have been identified, so these bhlh tfs might be involved in the ja response network and defense against tylcv infection.
slybhlh <dig>  a member of the ib subfamily was up-regulated in the r line, and down-regulated in the s line.
in barley and rice, this subfamily is involved in fe uptake and homeostasis  <cit> .
this suggests that slybhlh <dig> might perform some previously unknown function in tylcv infection.
therefore, our results will pave the way for studies of new functions of bhlh genes in tylcv infection and will further our understanding of this gene family under other biotic and abiotic stresses in tomato.
CONCLUSIONS
an extensive analysis of the tomato bhlh genes was performed, identifying  <dig> bhlh tfs in the entire tomato genome.
these genes can be divided into  <dig> subfamilies using phylogeny, protein motifs, and gene structures.
this phylogenetic analysis is in consistent with previous results.
the members of subfamilies may share conserved functions not shared by other species.
the pattern of expression of slybhlh genes was observed in r and s lines infected with tylcv.
results showed that slybhlh <dig> might be involved in the tylcv infection by vigs.
in summary, this is the first comprehensive and systemic analysis of bhlh transcription factors in tomatoes, and the results of this study revealed the importance of bhlh genes during tylcv infection.
they may also provide new opportunities for the investigation of previously unknown mechanisms by which tomatoes tolerate tylcv infection.
furthermore, our results have established a solid foundation for future studies of other biotic and abiotic stresses using biochemical and physiological approaches that will probably reveal the functional significance of this family in tomato.
additional files
additional file 1: table s <dig>  the primer sequences used for quantitative real-time pcr amplification of actin and slybhlh <dig>  table s <dig>  summary of bhlh genes in tomato genome.
table s <dig>  the bhlh domain sequences of tomato.
table s <dig>  the orthologous bhlh genes in tomato, potato and arabidopsis genome.
table s <dig>  the paralogous bhlh genes in tomato, potato and arabidopsis genome.
table s <dig>  the expression of  <dig> slybhlh genes in ty- <dig> and  <dig> after tylcv infection.
additional file 2: figure s <dig>  alignment of all the bhlh domain of tomato proteins.
shown at the top are the boundaries used in this study to distinguish the dna-binding basic region, the two a-helices and the variable loop region.
additional file 3: figure s <dig>  the nj phylogenetic tree of bhlh proteins motif using phylip software.
additional file 4: figure s <dig>  the nj phylogenetic tree and conserved motif compositions of tomato.
the neighbor-joining tree of tomato bhlh genes and their motif locations.
additional file 5: figure s <dig>  the nj phylogenetic tree and slybhlh gene structure of tomato.
the yellow and green blocks indicate introns and exons, respectively.
jinyan wang and zhongze hu contributed equally to this work.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
jyw, blz, wgy performed the data analysis and drafted the manuscript.
tmz, blz, ywy, and jyw participated in the analysis of the data.
jyw, lxt, zzh, mly, and tzc performed the experiments.
all authors approved the final version of the manuscript.

 the basic/helix-loop-helix proteins have dna-binding and dimerization capabilities. they are a superfamily of transcription factor that have been found to have many different functions in essential physiological and developmental process in animals and plants. the bhlh domain contains approximately amino acids with two functionally distinct regions, the basic region and the hlh regions. the basic region was amino acids long and typically included six basic residues. it was located at the n terminus of the domain and functions as a dna binding motif. the hlh region contains two amphipathic α helices separated by a loop region of variable length. the hlh region acts as a dimerization domain and allows the formation of homodimers or heterodimers. among all bhlh motifs, amino acids have been found to be highly conserved in organisms ranging from yeast to mammals. outside of these conserved bhlh domains, the proteins exhibited considerable sequence divergence. some bhlh proteins have been shown to bind to the sequences containing the core element known as the e box, with the most common form of g-box. the nucleotides flanking the core element may also have a role in binding specificity. based on the phylogenetic relationships, dna-binding motifs, and functional properties, the bhlh tfs family has been divided into six main groups in metazoans. in brief, group a bhlh proteins can bind to the cagctg core sequences of e-boxes. group b includes a large number of functionally proteins and bind to the g-box sequence cacgtg. group c contains an additional protein-protein interaction region and binds to acgtg or gcgtg sequences. group d proteins have the hlh region but lack the basic dna binding domain. group e proteins have pro or gly residues within the basic region and can bind preferentially to a typical sequence, cacgng. group f consists of the coe domain; they have diverse sequences compared with other groups and another domain for dimerization and dna binding. only a small number of plant bhlh proteins have been characterized functionally, far fewer than in animals. in arabidopsis, bhlh-encoding genes which were divided into subfamilies according to their phylogenetic relationships have been identified from the analysis of genome sequences. a total of and bhlh tfs have been identified in the rice and chinese cabbage genomes, respectively. these have been divided into and subfamilies, respectively. phylogenetic analysis showed that the plant bhlh proteins comprised subfamilies, of which were present in the common ancestors of extant mosses and vascular plants. most bhlh proteins identified so far have been functionally characterized in arabidopsis, and their roles have been shown to include regulation of fruit dehiscence, anther and epidermal cell development, hormone signaling, and stress responses. tomato is an economically important vegetable worldwide. the annual global production of tomato in was more than million tons including million tons in china. the tomato genome has been sequenced and assembled by the international tomato genome sequencing project, because tomatoes are economically important and it is model species for the study of fruit ripening. a high-quality genome sequence for domesticated tomato and more than proteins have been obtained. tomato yellow leaf curl virus is the most widespread and currently ranks 3rd among the most economically and scientifically most important plant viruses worldwide. the symptoms of tylcv infection in young plants include stunted growth, upward curling of leaf margins, marked reduction in leaf size, mottling and yellowing of young leaves, and flower abscission, leading to severe yield loss. currently five major loci resistant to tylcv have been identified from different wild tomato relatives, ty- ty- and ty- from s.chilense, ty- from s.habrochaites, and ty- from s.peruvianum. among them, ty- and ty- were found to be allelic and have been cloned. ty- and ty- were found to be allelic and have been cloned. they are rna-dependent rna polymerases and may be involved in rna silencing. in addition, ty- ty- and ty- have been mapped to chromosomes and respectively, using molecular markers. cdna library comparisons of susceptible and resistant tomato lines before and after tylcv infection showed approximately genes that are preferentially expressed in a tomato line with a resistance introgressed from s. habrochaites. using whole transcriptome sequencing of the tylcv-resistant tomato breeding line cln2777a and tylcv-susceptible tomato breeding line tmxa48-4-, and genes were found to be differentially expressed in the r and s tomato lines, respectively. in tomatoes, lefer, a bhlh protein encoded by solyc06g. slybhlh was the first identified regulator of iron nutrition in plants. lefer plays an important role in the fe-deficiency response of tomatoes. style, encoded by solyc02g. slybhlh is the major quantitative trait locus responsible for style length; this important floral attribute has been shown to be associated with the evolution of self-pollination and was cloned in cultivated tomatoes. however, the tomato bhlh protein family has not been analyzed at a genome-wide level, and the phylogenetic relationship of this protein family remains poorly understood. in this study, a total of slybhlh genes were identified in the tomato genomic sequence and phylogenetic analyses were carried out to evaluate the relationships among these genes. changes in global expression pattern of slybhlh genes in r and s lines infected by tylcv were analyzed to provide insight into the regulation of response to tylcv. the expression of slybhlh exhibited a variety of expression patterns, suggesting a novel layer of regulation for the response to tylcv in tomato. methods database search for bhlh genes the pfam database was used to screen the genome of tomato and potato. proteins with helix-loop-helix dna-binding domains were used to identify the putative bhlh proteins in tomato and potato using the hidden markov model. the hmmsearch tool, with an expected value cut-off of was used to identify the proteins. these sequences were then verified using the smart tool. the arabidopsis thaliana bhlh proteins were retrieved from the tair database using a previous report. phylogenetic analysis and identification of conserved motifs and gene structure the complete amino acid sequences were screened against the pfam database to identify the domains of bhlh transcription factors. mega software was used to construct neighbor-joining distance trees using tomato bhlh protein domain sequences. the bootstrap was set as replicates, which provided information regarding their statistical reliability. meanwhile, the nj method of the phylip software was also used with bootstrap of replicates to create another phylogenetic tree to validate the results from the nj method by mega software. a phylogenetic tree of all the identified bhlh protein domains was also constructed. the identified bhlh domains were aligned using a clustalx program with default settings. to identify the conserved motifs in tomato bhlh proteins, the multiple expectation-maximization for motif elicitation program version. was used with default parameters, except for the following parameters: optimum motif width was set to ≥ and ≤100; the maximum number of motifs was set to identify ten motifs. meme software was used to search for conserved motifs in the complete amino acid sequences of bhlh proteins. the coding domain sequences and dna sequences of tomato bhlh genes were used to assess gene structure using gsds. collinear correlations of bhlh genes in the tomato, potato, and arabidopsis genomes orthomcl program was used to identify the orthologous and paralogous genes in tomatoes, potatoes and arabidopsis. briefly, the tools blastp, with an e-value ≤ 1e− and orthomclpairs were used to find orthologs, inparalogs and coorthologs in these three species. the circos tool was used to link these genes to chromosomes. in addition, the relationships of orthologous and paralogous genes in these three species were also shown using the circos tool. the bhlh genes in tomato were searched for duplication events. chromosome distribution and gene duplications to determine the physical locations of bhlh genes, the starting and ending positions of all bhlh genes on each chromosome were obtained from the tomato database. the mapinspect software was used to draw the images of the locations of the tomato bhlh genes. we used the plant genome duplication database to retrieve the duplicate chromosomal blocks and then identify the bhlh genes in the duplication block which allowed us to identify duplicate tomato bhlh genes. the pgdd is a public database used to identify and catalogue plant genes in terms of intra-genomic or cross-genomic syntenic relationships. rna data collection and data mining transcriptomic data of tylcv-resistant breeding line, cln2777a and susceptible breeding line, tmxa48-4- with uninfected and mixed infection samples of and 7 days post infection were downloaded from ncbi sra database and analyzed as described in a previous study. enrichment of gene ontology categories was performed with an agrigo analysis toolkit using the topgo ‘elim’ algorithm for the aspects ‘biological process’ and ‘subcellular localization’. the selected categories were sorted from the lowest to the highest p value. validation of differentially expressed genes by quantitative rt-pcr four bhlh genes with differentially expressed in r or s lines were selected and subjected to quantitative rt-pcr validation. primers for quantitative rt-pcr were designed using primer software and primer specificity was evaluated by blasting primer sequences against the ncbi database. pcr amplifications were performed in a real-time thermal cycler qtower / with 15 μl of final volumes containing  μl of cdna,  μl each primer, 6 μl of sterile water, and  μl sybr premix extaq™ ii kit. the conditions for amplification were as follows: 5 min of denaturation at 95°c followed by 40 cycles of 95°c for 10 s, 60°c for 20 s, and 72°c for 10 s. the expression levels of selected genes were normalized to α-tubulin expression. relative gene expression was calculated using the 2-ΔΔct method. three biological replicates were performed for each of the selected genes. validation of candidate genes with virus-induced gene silencing and cell death analysis the tobacco rattle virus mediated vigs system was used to silence a bhlh gene. briefly, ptrv-containing agrobacterium eha was cultured in liquid lb medium and resuspended in infiltration medium at an o.d. value of and cultured at room temperature for 4 h. three week old seedlings were infiltrated by pressure inoculation in the leaves with a needleless syringe. for the vigs experiments, agro infiltration was performed two weeks after trv inoculation. death cells were identified by staining with lacto-phenol trypan blue and dab as previous described. to visualize cell death, stem were stained by boiling in lacto-phenol trypan blue, followed by destaining with chloral hydrate. then the death cell was examined with ni-u microscope.  identification and classification of bhlh proteins in tomato to identify the putative bhlh proteins in the tomato genome, a hidden markov model search resulted in the identification of bhlh proteins. to verify the reliability of our criteria, we performed simple modular architecture research tool analysis of putative slybhlh protein sequences and found all of them had a typical bhlh domain. the number of bhlh tfs in tomatoes exceeded that of many metazoans and fungi, but was less than that found in some plants, such as rice, chinese cabbage, potatoes, soybeans and maize. in addition, the density of bhlh proteins in the entire tomato and potato genomes was found to be less than that in most plant species, such as arabidopsis and rice. most angiosperm plant lineages have experienced one of more rounds of ancient polyploidy. and the genomes of tomato and potato have undergone recent triplication events, whereas few individual tomato/potato genes remain triplicated. therefore, this event might lead to relatively fewer bhlh genes in tomato compared to other plant species. multiple sequence alignments, predicted dna-binding ability and conserved residues to examine sequence features of these tomato bhlh domains, multiple sequence alignment of the bhlh amino acid sequences were performed. there were four conserved regions in the bhlh domain sequences, including one basic region, two helix regions and one loop region. the basic regions have five basic residues, but five of these proteins did not have the basic region. the loop was found to be the most divergent region in terms of both length and amino acid composition. from the alignment, residues were identified that were identical in at least 50% of the tomato bhlh domains. among these residues, nine residues were present in more than 75% sequences.figure the characterization and distribution of the bhlh domains. a. sequence logo of the slybhlh domain by meme. the h e and r amino acids in the basic domain that are important for dna binding are indicated by stars. amino acids important for dimerization of the helix-loop-helix domain are indicated by arrows. b. distribution of amino acids in the bhlh consensus motif among tomatoes. the numbers of horizontal ordinate refer to the positions of the residues in the alignments of the studies. c. predicted dna-binding characteristics of the bhlh domain of slybhlh proteins. five residues, five residues, two residues and seven residues made up the basic region, the first helix region, the loop region and the second helix region, respectively. all of these conserved residues were consistent with previous studies. the leu- in the basic region was conserved in all bhlh proteins, suggesting that this residue is extremely important for promoting the formation of dimerization among bhlh proteins. the basic region of the bhlh domain can bind to dna and is critical for function. using the criteria described by massari and murre, the slybhlh proteins were divided into several categories based on sequence information in the n-terminal region of the bhlh domains. as was done with arabidopsis and chinese cabbage, the slybhlh proteins of tomato were also divided into two major groups according to 17 n-terminal amino acids within bhlh protein domain, including dna-binding and non-dna binding proteins. the dna-binding bhlhs were further divided into two groups with different predicted target sequences depending on the presence or absence of residues glu- and arg- in the basic region. group proteins had putative e-box-binding proteins with conserved glu-9/arg- residues and group proteins had non-e-box-binding proteins lacking these residues. the three residues in the basic region of the bhlh domain, his/lys- glu- and arg- were found to constitute the classic g-box-binding region. group can therefore be subdivided further into two subgroups: 1a whose proteins are predicted to bind g-boxes, and 1a whose three members are predicted to bind other types of e boxes. phylogenetic analysis of the bhlh transcription factor family to assess the evolutionary relationships of the slybhlh genes, an nj phylogenetic tree was generated using the multiple sequences alignments of the conserved bhlh tf domains in tomato and arabidopsis with a bootstrap value of. twenty-six subfamilies were identified according to the clades support values, topology of the tree, and classification of the arabidopsis. no slybhlh proteins in the xiii, ii, subfamilies relative to those of arabidopsis, therefore, tomato contained bhlh subfamilies in our analysis. to further validate the reliability of the nj tree with mega, nj and maximum parsimony analysis was also used to generate phylogenetic trees using phylip software. % of the slybhlh proteins with nj model using phylip software were placed into the same subfamilies as those in the nj tree with mega, indicating that both methods are in very good agreement.figure phylogenetic tree constructed from the neighbor-joining method using the bhlh transcription factor domain in tomato. the phylogenetic tree was constructed using mega software. the numbers are bootstrap values based on iterations. only bootstrap values larger than 50% support are indicated. in order to assess differences in protein structure, meme was used to identify conserved motifs in the tomato bhlh proteins. ten conserved motifs were identified and named motif through motif in general, the bhlh proteins were clustered in the same subfamilies and shared similar motif compositions, which indicated functional similarities among members of the same subfamilies. the tomato bhlh proteins were found to have a similar structure for every subfamily. the pattern of intron position can also provide important evidence to support phylogenetic relationships in a gene family. here, gsds tools were used to show the gene structures for slybhlh genes. among the tomato bhlh genes, the number of introns ranged from to and most members of the same subfamilies had similar intron/exon structures. for example, the members of subfamilies iii and xv each have only one intron. these results demonstrated that proteins within the same subfamily share close evolutionary relationships. collinear correlations of bhlh genes in tomatoes, potatoes, and arabidopsis the solanum lineage has experienced two consecutive genome triplications: one is ancient and shared with rosids, and the other more recent. in this study, the correlation between tomato, potato, and arabidopsis bhlh genes was analyzed using the orthomcl program. here, gene pairs were found to be orthologous between tomatoes and potatoes, but only orthologous pairs were found between tomato and arabidopsis. these results were consistent with the close relationship of tomatoes and potatoes. among the orthologous gene pairs shared by tomatoes and potatoes, each tomato bhlh gene had one to four potato bhlh genes. these results demonstrated that bhlh tf genes in potato were duplicated accompanied with evolution processes. in addition, paralogous bhlh gene pairs were also analyzed. a total of and bhlh gene pairs were identified in arabidopsis, tomatoes, and potatoes, respectively. visualization of the relationships of paralogous and paralogous bhlh genes among these three species was performed using the circos software.figure comparative analysis of synteny and expansion of bhlh genes. a. twelve tomato chromosomes and of potato chromosomes maps were based on the orthologous and paralogous pair positions and demonstrate highly conserved synteny. the red lines represent the orthologous bhlh genes between tomato and potato. the blue and yellow lines represent the paralogous bhlh genes in potato and tomato, respectively. b. twelve tomato and arabidopsis chromosome maps were based on the orthologous and paralogous pair positions and demonstrate highly conserved synteny. the red lines represent the orthologous bhlh genes between tomato and arabidopsis. the blue and yellow lines represent the paralogous bhlh genes in tomato and arabidopsis, respectively. chromosome distribution and gene duplication of the bhlh tf family the physical map positions of the bhlh genes on tomato chromosomes were identified. among the bhlh tfs, were mapped onto the twelve tomato chromosomes except slybhlh most bhlh tfs were found on chromosome and. in contrast, there are only and bhlh tf genes on chromosome and respectively. furthermore, bhlh genes were found to be mapped on the chromosomes with an obviously uneven distribution, and some bhlh genes gathered on part of the chromosome. relative high densities of bhlh genes were observed in some chromosomal regions, including the bottom of chromosomes and for example, genes clustered in the end of chromosome and with density of and genes per mb, respectively. in contrast, several large chromosomal regions lacked bhlh genes, such as the top half of chromosomes and and the central section of chromosomes and.figure distribution of bhlh genes on the twelve tomato chromosomes using mapinspect software. only slybhlh could not be anchored onto a specific chromosome. the scale is in megabases. the duplicate bhlh genes are connected with blue lines. previous reports have analyzed duplication events in rice and chinese cabbage. in the current analysis, we first retrieved the genome chromosome blocks in the tomato with pgdd database and identified the duplicate blocks. a total of duplicate bhlh genes pairs were located in these blocks. these duplication bhlh genes are derived from the same subfamily, indicating that members of some bhlh subfamilies originated from the duplication events. differential expression of slybhlh genes in response to tylcv infection the rna-seq technology has been shown to provide precise digital information related on gene expression and can discriminate genes of high sequence identity. using this technology, global gene expression changes of the leaves of two tomato breeding lines, tylcv-resistant cln2777a and tylcv-susceptible tmxa48-4- have been analyzed before and after tylcv infection with viruliferous whiteflies. a total of transcripts were detected from r and s lines by alignment to the tomato genome including the novel transcripts predicted in tomatoes. the expression levels of mapped genes were normalized with a value of fragments per kilobase of exon per million fragments mapped. if the fpkm value of gene was above zero, the gene was considered expressed. on the basis of this criterion, the expression level of slybhlh genes was confirmed in the s and r line with and dpi. a subset of slybhlh genes was expressed under both conditions. go enrichment analysis revealed that the products of most of slybhlh genes were localized in the nucleus. in addition, some biological processes such as ‘regulation of macromolecule metabolic process’, ‘regulation of gene expression’ and ‘regulation of metabolic process’ were overrepresented among all slybhlh genes, indicating that the bhlh genes were involved in transcription and metabolic regulation.figure the expression levels and related functions of bhlh genes. a. venn’s diagram of the cross comparison of expressed genes based on rna-seq data. b. gene ontology enrichment analysis of expressed bhlh genes involved in the biological processes. c. subcellular location of bhlh genes. out of slybhlh genes in the current study, only four were differentially expressed before and after tylcv infection. among the four differentially expressed slybhlh genes, slybhlh was up-regulated in the r line and down-regulated in the s line after tylcv infection. the expression level of four differentially expressed genes was determined using quantitative rt-pcr to validate the gene expression data from rna-seq. the results demonstrated that all tested genes revealed a similar trend of transcript accumulation as in rna-seq analysis.table differentially expressed slybhlh genes between before and after tylcv infection ty-2 4840 dpi dpi fold change dpi dpi fold changequantitative rt-pcr was used to measure the relative expression levels of five pathogen resistance related genes in the r and s lines, with tomato α-tubulin as an internal reference. a. slybhlh b. slybhlh c. slybhlh d. slybhlh vigs validation of slybhlh gene to investigate the role of tylcv as it related to resistant in tomato, the slybhlh gene was further challenged with tylcv after vigs at the cotyledon stage. one month after agroinfiltration, the success of the trv silencing system was confirmed by the appearance of cell death in the leaves of s plantlets treated with ptrv and ptrv2- slybhlh. we also observed the cell death development under ni-u microscope and found that the vigs lines with slybhlh gene triggered a rapid cell death response in comparison with empty vector lines by trypan blue and dab staining. a quantitative rt-pcr analysis showed there were significantly fewer slybhlh transcripts in slybhlh131-silenced tomato leaves during trv infection, indicating that slybhlh was effectively silenced in tomatoes.figure validation of the bhlh gene with virus-induced gene silencing. a. cotyledon agroinfiltration of trv vectors was carried out in the s plantlets at the cotyledon stage. r plants treated with phytoene desaturase gene silencing constructs ptrv and ptrv - showed bleached areas in leaflets. s plants treated with ptrv and ptrv vectors showed the normal phenotype. s plantlets treated with slybhlh gene silencing constructs ptrv and ptrv - showed a cell death phenotype. arrow indicates cell death. b. trypan blue staining of vigs and empty vector lines leaves. c. dab staining of vigs and empty vector lines leaves. d. relative expression of slybhlh transcript using real-time rt-pcr analysis in the vigs-treated lines 20 days after agroinfiltration with trv vectors. tomato α-tubulin was used as an internal reference. error bars represented se of three biological replicates and asterisks indicate significant differences by student’s t test for p < . discussion the solanaceae is one of the largest angiosperm genera and includes many different vegetables consumed by humans. in recent years, many plants in this genera have been sequenced, including tomatoes, potatoes, peppers, and tobacco. with the rapid development in bioinformatic analyses, the information stored in various genomes may be explored to elucidate the mechanisms that regulate the development and response to biotic and abiotic stresses. tomatoes are major crop plant and a model system for fruit development. as of tomato production had doubled from million tons in. however, tomato yellow leaf curl virus disease causes huge losses in tomato production worldwide. this disease is caused by different related be gomovirus species. a previous study extended our basci understanding of the response of tomato to tylcv infection by comparing whole transcriptome expression changes between a tylcv-resistant line and a tylcv-susceptible line. in present study, bhlh transcription factor genes were identified in tomatoes and their differential expression was analyzed in r and s lines before and after tylcv infection. four differentially expressed genes were identified in r and s lines, in which slybhlh and slybhlh were derived from the chromosome and and slybhlh and slybhlh were mapped on the chromosome these genes are not located in the region of five known loci, so they are not the most important genes but might be involved in the regulation network of tylcv resistance. phylogenetic analysis of the bhlh domain allows division of the slybhlh family into subfamilies. the clustering of the members within these subfamilies was further supported by additional analysis with regard to other criteria, such as predicted dna binding capacity and sequence specificity, exon/intron distribution pattern with the domain. these data support the general conclusion that members within subfamilies may have recent common evolutionary origins resulting from various genomic duplication events. they may have related molecular functions. the bhlh subfamily iiid transcription function function redundantly to negatively regulate jasmonate responses in arabidopsis. however, the strong sequence diversity outside of the bhlh domain across the members of the slybhlh family suggests that the expansion of this family in tomatoes involved extensive domain shuffling, as in other organisms. however, the non-bhlh amino acid motifs are conserved in each of the bhlh subfamilies, suggesting that the conservation of these extra domains during plant evolution may have been essential to the function of the bhlh proteins in their respective subfamilies. the core dna binding domain of the bhlh proteins contained the basic region of the bhlh domain and these residues were found to recognize and bind to the core hexanucleotide. the amino acid sequence in this region provides the major subdivision of the bhlh family, dividing these proteins into those that are predicted to bind dna and those that are not. key residues in this region confer the capacity to discriminate the variants of the hexanucleotide motif with the canonical e-box and non-e-box motifs. additional residues within the basic region confer further dna binding site sequence selectivity. these include g-box and non-g-box core motifs. the current analysis showed there are e-box binding bhlh proteins in tomatoes, and this ration was lower than in arabidopsis and chinese cabbage, indicating that many different binding motifs exist in the tomato bhlh protein family. most bhlh proteins identified so far were functionally characterized in arabidopsis, and their roles include plant development, fruit dehiscence, phytochrome signaling, hormone signaling and stress responses, such as cold, heat, abscisic acid, jasmonic acid, and the light signaling pathway. only two bhlh genes have been functionally characterized in tomato. one is the fer gene, which is involved in the response to iron acquisition and supply. another is style, that is associated with the evolution of self-pollination; its natural genetic variation in the promoter is responsible for evolution from allogamy to autogamy in the cultivated tomatoes. in the current rna-seq data set, the fer gene was not expressed in any of the four samples, and the style gene was expressed but not differentially expressed in r or s lines after tylcv infection. members of the same plant bhlh subfamily are frequently involved in the same biological processes. usually the functions of these proteins overlap, causing them to be partially or totally redundant. the characterized function of bhlh in other species can help the user to predict the function of tomato bhlh in the same subfamily. jasmonates are lipid-derived hormones that regulate plant responses to stresses such as wounding and pathogens invasion. ja negatively regulates plant growth and is considered to modulate the distribution of energy to defense responses. in arabidopsis, the bhlh subfamily iiid members act as transcription repressors and function redundantly to negatively regulate ja response. in tomatoes, eight members of the iii subfamily have been identified, so these bhlh tfs might be involved in the ja response network and defense against tylcv infection. slybhlh a member of the ib subfamily was up-regulated in the r line, and down-regulated in the s line. in barley and rice, this subfamily is involved in fe uptake and homeostasis. this suggests that slybhlh might perform some previously unknown function in tylcv infection. therefore, our results will pave the way for studies of new functions of bhlh genes in tylcv infection and will further our understanding of this gene family under other biotic and abiotic stresses in tomato.  an extensive analysis of the tomato bhlh genes was performed, identifying bhlh tfs in the entire tomato genome. these genes can be divided into subfamilies using phylogeny, protein motifs, and gene structures. this phylogenetic analysis is in consistent with previous results. the members of subfamilies may share conserved functions not shared by other species. the pattern of expression of slybhlh genes was observed in r and s lines infected with tylcv. results showed that slybhlh might be involved in the tylcv infection by vigs. in summary, this is the first comprehensive and systemic analysis of bhlh transcription factors in tomatoes, and the results of this study revealed the importance of bhlh genes during tylcv infection. they may also provide new opportunities for the investigation of previously unknown mechanisms by which tomatoes tolerate tylcv infection. furthermore, our results have established a solid foundation for future studies of other biotic and abiotic stresses using biochemical and physiological approaches that will probably reveal the functional significance of this family in tomato. additional files additional file 1: table s the primer sequences used for quantitative real-time pcr amplification of actin and slybhlh table s summary of bhlh genes in tomato genome. table s the bhlh domain sequences of tomato. table s the orthologous bhlh genes in tomato, potato and arabidopsis genome. table s the paralogous bhlh genes in tomato, potato and arabidopsis genome. table s the expression of slybhlh genes in ty- and after tylcv infection. additional file 2: figure s alignment of all the bhlh domain of tomato proteins. shown at the top are the boundaries used in this study to distinguish the dna-binding basic region, the two a-helices and the variable loop region. additional file 3: figure s the nj phylogenetic tree of bhlh proteins motif using phylip software. additional file 4: figure s the nj phylogenetic tree and conserved motif compositions of tomato. the neighbor-joining tree of tomato bhlh genes and their motif locations. additional file 5: figure s the nj phylogenetic tree and slybhlh gene structure of tomato. the yellow and green blocks indicate introns and exons, respectively. jinyan wang and zhongze hu contributed equally to this work. competing interests the authors declare that they have no competing interests. authors’ contributions jyw, blz, wgy performed the data analysis and drafted the manuscript. tmz, blz, ywy, and jyw participated in the analysis of the data. jyw, lxt, zzh, mly, and tzc performed the experiments. all authors approved the final version of the manuscript. 
however, no systemic analysis of the bhlh transcription factor family has yet been reported in tomatoes. distribution mapping results showed bhlh genes were localized on the tomato chromosomes. among the bhlh genes from the tomato genome, bhlh genes were detected in the tylcv-susceptible and resistant tomato breeding line before and after tylcv infection. as anticipated, gene ontology analysis indicated that most bhlh genes are related to the regulation of macromolecule metabolic processes and gene expression. only four bhlh genes were differentially expressed between and dpi. virus-induced gene silencing of one bhlh genes slybhlh in resistant lines can lead to the cell death. one of which bhlh genes, slybhlh was found to be involved in the tylcv infection through qrt-pcr expression analysis and vigs validation. a total of bhlh transcription factors were identified from the entire tomato genome. the isolation and identification of these bhlh transcription factors facilitated clarification of the molecular genetic basis for the genetic improvement of tomatoes and the development of functional gene resources for transgenic research. in addition, these findings may aid in uncovering an unexplored mechanism during the tylcv infection in tomatoes. phylogenetic analysis of bhlh domain sequences from arabidopsis and tomato facilitated classification of these genes into subfamilies. in the present study, bhlh transcription factor genes were identified. tomato yellow leaf curl virus threatens tomato production worldwide by causing leaf yellowing, leaf curling, plant stunting and flower abscission.
 the basic helix-loop-helix proteins are a superfamily of transcription factors that can bind to specific dna target sites. they have been well characterized in model plants such as arabidopsis and rice and have been shown to be important regulatory components in many different biological processes. however, no systemic analysis of the bhlh transcription factor family has yet been reported in tomatoes. tomato yellow leaf curl virus threatens tomato production worldwide by causing leaf yellowing, leaf curling, plant stunting and flower abscission.  a total of bhlh transcription factors were identified from the entire tomato genome. phylogenetic analysis of bhlh domain sequences from arabidopsis and tomato facilitated classification of these genes into subfamilies. the evolutionary and possible functional relationships revealed during this analysis are supported by other criteria, including the chromosomal distribution of these genes, the conservation of motifs and exon/intron structural patterns, and the predicted dna binding activities within subfamilies. distribution mapping results showed bhlh genes were localized on the tomato chromosomes. among the bhlh genes from the tomato genome, bhlh genes were detected in the tylcv-susceptible and resistant tomato breeding line before and after tylcv infection. as anticipated, gene ontology analysis indicated that most bhlh genes are related to the regulation of macromolecule metabolic processes and gene expression. only four bhlh genes were differentially expressed between and dpi. virus-induced gene silencing of one bhlh genes slybhlh in resistant lines can lead to the cell death.  in the present study, bhlh transcription factor genes were identified. one of which bhlh genes, slybhlh was found to be involved in the tylcv infection through qrt-pcr expression analysis and vigs validation. the isolation and identification of these bhlh transcription factors facilitated clarification of the molecular genetic basis for the genetic improvement of tomatoes and the development of functional gene resources for transgenic research. in addition, these findings may aid in uncovering an unexplored mechanism during the tylcv infection in tomatoes. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords genome-wide analysisbhlh transcription factortylcvvigstomatoissue-copyright-statement© the author 2015 