- Research article
- Open Access
Comparative analysis of the repertoire of G protein-coupled receptors of three species of the fungal genus Trichoderma
BMC Microbiology volume 13, Article number: 108 (2013)
Eukaryotic organisms employ cell surface receptors such as the seven-transmembrane G protein-coupled receptors (GPCRs) as sensors to connect to the environment. GPCRs react to a variety of extracellular cues and are considered to play central roles in the signal transduction in fungi. Several species of the filamentous ascomycete Trichoderma are potent mycoparasites, i.e. can attack and parasitize other fungi, which turns them into successful bio-fungicides for the protection of plants against fungal phytopathogens. The identification and characterization of GPCRs will provide insights into how Trichoderma communicates with its environment and senses the presence of host fungi.
We mined the recently published genomes of the two mycoparasitic biocontrol agents Trichoderma atroviride and Trichoderma virens and compared the identified GPCR-like proteins to those of the saprophyte Trichoderma reesei. Phylogenetic analyses resulted in 14 classes and revealed differences not only among the three Trichoderma species but also between Trichoderma and other fungi. The class comprising proteins of the PAQR family was significantly expanded both in Trichoderma compared to other fungi as well as in the two mycoparasites compared to T. reesei. Expression analysis of the PAQR-encoding genes of the three Trichoderma species revealed that all except one were actually transcribed. Furthermore, the class of receptors with a DUF300 domain was expanded in T. atroviride, and T. virens showed an expansion of PTH11-like receptors compared to T. atroviride and T. reesei.
Comparative genome analyses of three Trichoderma species revealed a great diversity of putative GPCRs with genus- and species- specific differences. The expansion of certain classes in the mycoparasites T. atroviride and T. virens is likely to reflect the capability of these fungi to establish various ecological niches and interactions with other organisms such as fungi and plants. These GPCRs consequently represent interesting candidates for future research on the mechanisms underlying mycoparasitism and biocontrol.
Fungi are eukaryotes and include organisms with important ecological and economic roles. The relatively simple structure and the ease of cultivation and genetic manipulation make fungi interesting eukaryotic models for studying fundamental biological processes. They share important features with even mammalian cells such as conserved signal transduction pathways that regulate cell function [1, 2]; thus studying fungal signaling and environmental sensing contributes to our knowledge on conserved basic molecular principles of life.
Communication of cells with each other and with their environment is crucial for survival of organisms. Consequently, ingenious mechanisms of sensing environmental signals and elaborated ways of adaption to the environment evolved . Cell surface receptors connect the cell to the environment by functioning as sensors. Among these receptors, G protein-coupled receptors (GPCRs) comprise the largest class with roles in virtually every physiological function . GPCRs have a common domain structure containing seven stretches of hydrophobic amino acids spanning the cytoplasmic membrane connected by intra- and extracellular loops with the N-terminus located outside of the cell and the C-terminus within the cytoplasm . The classic paradigm is based on a physical interaction of the GPCR with an intracellular GÎ± subunit once the receptor is activated by ligand binding which leads to dissociation of GÎ± from GÎ²Î³ subunits . Both signalling units then regulate activities of downstream effectors [7–9]. In eukaryotic organisms a plenty of different GPCRs is facing a small amount of G proteins. If G proteins were the only transmitters of GPCR-mediated signaling, this unequal ratio seems to limit the specificity of signal transduction. In recent years several intracellular partners other than G proteins were identified that are capable of mediating signals originating from these receptors. These include arrestins, G protein-coupled receptor kinases, small GTP-binding proteins, and many more [10–13]. Accordingly, GPCRs are extremely diverse in sequence and function and missing genome sequence information and constraints in structure prediction for a long time impaired research on these proteins. Although pheromone- and nutrient- sensing GPCRs have been studied extensively in yeast and some filamentous fungi [14–26] far more GPCRs remain to be identified and characterized.
The fungal genus Trichoderma comprises saprophytic and mycoparasitic species, and species interacting with plants and animals . Because of these versatile lifestyles and the variety of interactions with other organisms, Trichoderma fungi are valuable models for studying organismic cross-talk and signaling. Studies on heterotrimeric G proteins revealed a multitude of processes being regulated by these signal transduction compounds in Trichoderma. The class I adenylate cyclase-inhibiting as well as the class III adenylate cyclase-activating GÎ± subunits regulate vegetative growth and conidiation of the fungus and affect processes relevant for mycoparasitism , i.e. a lifestyle where Trichoderma parasitizes other fungi. Trichoderma atroviride Tga1 as well as Tga3 govern the production of extracellular chitinases and antifungal metabolites, and Tga3 is essential for transmitting signals that regulate the recognition of the host fungus and attachment to its hyphae. Both, T. atroviride âˆ†tga1 as well as âˆ†tga3 mutants, are unable to overgrow and lyse host fungi [29–31], while Trichoderma virens TgaA regulates mycoparasitism in a host-specific manner . For T. virens âˆ†tgaB mutants missing the class II GÎ±-encoding gene, unaltered growth, conidiation, and mycoparasitic activity have been reported . In the saprophyte Trichoderma reesei, the heterotrimeric G protein pathway is crucial for the interconnection of nutrient signaling and light response. Besides the GÎ± subunits GNA1 and GNA3, which transmit signals positively impacting cellulase gene expression, GNB1 (GÎ²), GNG1 (GÎ³) and the phosducin PhLP1 influence light responsiveness, glycoside hydrolase expression and sexual development [33, 34].
Here we present an exploration of the genomes of the two mycoparasites T. atroviride and T. virens and identify members of the G protein-coupled receptor family from the entire deduced proteomes. The identified proteins are classified and compared to those encoded in the saprophyte T. reesei and several other fungi. In contrast to the presence of only three GÎ± subunits, one beta and one gamma subunit in each of the genomes of the three Trichoderma species, our analyses revealed a great diversity of GPCRs and differences both between the three Trichoderma species and between Trichoderma and other fungi.
Results and discussion
Identification of G protein-coupled receptor-like proteins in the genomes of three Trichoderma species
The T. atroviride, T. virens and T. reesei genome databases were searched for putative GPCRs using a homology (BLAST)-based strategy. Together with the putative GPCRs identified in the genome of Neurospora crassa and Phytophtora sojae GPR11 , the 18 GPCRs previously identified in Aspergillus spp.  and the three new GPCRs predicted in the Verticillium genome  were used in a BLASTP search against the predicted proteomes of the following species of the Sordariomycetes (Magnaporthe grisea, Podospora anserina, Chaetomium globosum, Fusarium graminearum, Nectria haematococca, T. reesei, T. atroviride and T. virens), a subgroup within the Ascomycota. In an analogous manner, the PTH11 receptor of M. grisea[14, 37] was used as a query. All consequently identified GPCR-like proteins were next used as a query in similar BLAST searches of the proteomes of the other species. In the end each possible combination was tested. By additionally applying a HMM-based approach, which is suitable for detecting candidates lacking significant sequence similarity to known GPCR-like proteins and therefore may escape detection by BLAST-based homology searches, two additional proteins of the PTH11-like class could be identified (Triat86665, Trive78137).
All identified Trichoderma proteins were evaluated for the typical topology of seven transmembrane regions and, if conducive, a manual editing of candidate GPCR sequences was performed including movement of exon-intron boundaries and sequence extension or truncation. This total set of analyses resulted in the identification of 65 and 76 putative GPCRs in T. atroviride and T. virens, including 38 and 52 PTH11-like receptors, respectively, which are facing 58 predicted GPCRs in the T. reesei genome (Table 1). Among the PTH11-like receptors, a protein exhibiting 15 transmembrane domains was found in all three Trichoderma species. An orthologue of this putative GPCR has previously been identified in M. grisea and A. nidulans suggesting conservation of this particular receptor.
Phylogenetic analysis of the identified Trichoderma GPCR-like proteins
Previous studies led to the categorization of fungal GPCRs into the following classes: pheromone receptors, carbon sensors, putative nitrogen sensors, cAMP receptor-like proteins, GPCRs with an RGS domain, GPCRs related to rat growth hormone releasing factor, mPR-like GPCRs, microbial opsins and those related to M. grisea PTH11 [1, 2, 14]. Recently, this classification has been extended by three novel classes whose members show similarity to PTM proteins (putative tumor necrosis factor receptors), to GPR89A of higher eukaryotes, and to family C-like GPCRs (metabotropic glutamate/pheromone receptors of Gallus gallus), respectively .
A phylogenetic analysis of all putative GPCRs identified in this study including those previously described for T. reesei[38, 39] revealed that the Trichoderma proteins were distributed over 14 classes including PTH11-like GPCRs and putative receptors similar to P. sojae GPR11 (Figure 1, Table 1). Phylogeny also showed that the orthologous proteins from T. atroviride, T. virens and T. reesei mainly formed the topologies ((Tr, Tv) Ta) and ((Ta, Tv) Tr) with 14 and 9 cases, respectively, whereas the ((Ta, Tr) Tv) topology resulted only once (Figure 1). This suggests that some of the GPCRs of T. virens are more related to those of T. atroviride and some are more related to those of T. reesei. This is in accordance to the phylogeny of these species based on other genes showing that T. atroviride resembles the more ancient state of Trichoderma and that both T. virens and T. reesei evolved later . Accordingly, comparative genome analysis showed that the lineage to T. reesei appears to have lost a significant number of genes present in T. atroviride and maintained in T. virens.
Trichoderma members of classes I to VII of fungal GPCRs
Two putative pheromone receptors are encoded in the genomes of the three Trichoderma species analyzed. Similar to other fungi, these proteins group to classes I and II of fungal GPCRs (Figure 1, Additional file 1), respectively, and harbor the typical STE2 (pfam02116; Triat36032, Trive147400, Trire64018) and STE3 (pfam02076; Triat147894, Trive40681, Trire57526) domains. Functional analysis of the pheromone receptors of T. reesei (H. jecorina) showed that HPR1 and HPR2 confer female fertility in their cognate mating types, mediate induction of fruiting body development, and are involved in ascosporogenesis . While sexual crossing remains to be experimentally shown for T. atroviride and T. virens, a respective MAT1-2 mating type locus is present in their genomes and the corresponding teleomorphs, Hypocrea atroviridis and Hypocrea virens, have already been described [41, 42].
There is no direct sequence homologue of the class III carbon-sensing GPCRs Gpr1 of Saccharomyces cerevisiae and GPR-4 of N. crassa[21, 43, 44] in Trichoderma. Nevertheless, we could identify a 7-transmembrane domain protein in T. atroviride (Triat246916), T. virens (Trive29548) and T. reesei (Trire59778) sharing sequence and structural similarity with Aspergillus nidulans GprC, GprD and GprE, and GprC and GprD of Aspergillus fumigatus and Aspergillus oryzae, which have previously been described as class III GPCRs . GprD negatively regulates sexual development in A. nidulans and A. fumigatus and GprC and GprD of A. fumigatus are furthermore involved in integrating and processing stress signals via modulation of the calcineurin pathway [45, 46]. Recently, GprD was further shown to be involved in the sensing of oxylipins in A. nidulans and A. flavus. Due to the absence of a locus similar to that of N. crassa GPR-4 in the T. reesei genome, it has been postulated that T. reesei does not possess a class III GPCR. Trire59778 was instead grouped to the cAMP receptor-like class . However, structural analyses of receptors of classes III and V revealed distinct topologies: whereas class III members display seven transmembrane regions at their amino-terminal end and a long carboxy-terminal cytoplasmic domain, class V receptors exhibit five domains at the N-terminal end, a long intracellular loop and two helices next to the C-terminus . Consistent with a clustering of Triat246916, Trive29548 and Trire59778 with A. nidulans GprC, GprD and GprE in the phylogenetic analysis (Additional file 1), the Trichoderma proteins clearly share the topology of class III members and contain a Git3 (pfam11710; G protein-coupled glucose receptor) domain. Whether these proteins actually are implicated in glucose sensing, remains to be elucidated.
Fungal GPCRs with similarity to Schizzosaccharomyces pombe Stm1 have been designated as class IV. The Stm1 receptor has been previously shown to be required for proper recognition of nitrogen starvation signals and to couple to the Gpa2 GÎ± subunit in S. pombe. This class of GPCRs, all containing PQ-loop repeats, is well conserved in filamentous fungi , although their function remains elusive. Two PQ-loop containing 7-transmembrane proteins grouping to class IV are encoded in the mycoparasites T. atroviride and T. virens (Figure 1, Table 1) which is consistent with previous reports on T. reesei[38, 39]. Interestingly, one of the two class IV members of T. atroviride, Triat300620, has been found in an EST-based study to be expressed exclusively under mycoparasitic conditions (i.e. in direct confrontation with the host fungus Rhizoctonia solani) . This transcriptome analysis further revealed that T. atroviride faces stress from nitrogen limitation when it is confronted with a fungal host accompanied by an up-regulation of genes encoding proteolytic enzymes. Consequently, oligopeptides emerging from an initial degradation of the host by secreted proteases have been suggested as signals for nitrogen deficiency by binding to the Stm1-receptor in a ligand-receptor-specific manner . A possible role of Triat300620 in nitrogen signaling during mycoparasitism is further supported by the fact that T. atroviride knock-out mutants missing the Tga3 GÎ± protein (orthologue of S. pombe Gpa2) are completely deficient in mycoparasitism, e.g. unable to attack and parasitize host fungi .
The class V of fungal GPCRs comprises cAMP receptor-like (CRL) proteins that are distantly related to the four cAMP receptors of Dictyostelium discoideum[1, 2]. Similar to T. reesei, four CRL proteins harboring a Dicty_CAR (pfam05462) domain were identified in the genomes of the two mycoparasitic Trichoderma species T. atroviride and T. virens (Figure 1, Table 1). Two of these (Gpr1/ Triat160995 and Gpr2/ Triat 50902) have been functionally characterized in T. atroviride. While mutants silenced in the gpr2 gene did not show any phenotypic alterations [28, 38], gpr1 mutants were unable to attach to host hyphae and to respond to host fungi with the production of cell wall-degrading enzymes. Besides these defects in mycoparasitism-relevant activities, Gpr1 further affects vegetative growth and conidiation of T. atroviride. As Gpr1 did not interact with any of the three T. atroviride GÎ± proteins (Tga1, Tga2, or Tga3) in a split-ubiquitin yeast-two-hybrid assay , signal transduction in a G protein-independent manner cannot be ruled out at the moment.
Members of class VI of fungal GPCRs are characterized by the presence of both 7-transmembrane regions and an RGS (regulator of G protein signaling) domain in the cytoplasmic part of the proteins. They show similarity to Arabidopsis thaliana AtRGS1 which modulates plant cell proliferation via the Gpa1 GÎ± subunit . In contrast to other filamentous ascomycetes like F. graminearum, N. crassa, A. nidulans, A. fumigatus, A. oryzae, Verticillium spp. and M. grisea, which possess only one or two members of class VI [1, 2], three putative RGS domain-containing GPCRs could be identified in both T. reesei[38, 39] and the two mycoparasitic species T. atroviride and T. virens (Table 1).
A putative receptor distantly related to mammalian GPCRs like the rat growth hormone-releasing factor receptor has been initially identified in the M. grisea genome . Similar to closely related fungi like N. crassa and F. graminearum one orthologue with more than 50% amino acid identity to MG00532 is encoded in the genomes of T. atroviride, T. virens and T. reesei which accordingly was assigned to class VII (Table 1).
The PAQR family is expanded in mycoparasitic Trichoderma species
Receptors responding to progesterone and adiponectin as ligands have previously been classified as progestin-adipoQ receptors (PAQR , a group of 7-transmembrane proteins lacking significant sequence similarity to any previously described GPCRs but with ancient evolutionary roots. The PAQR family also includes prokaryotic hemolysin-type proteins and members have been identified throughout the eukaryotic kingdom including 11 paralogues in mammals . In S. cerevisiae the PAQR family members Izh1p, Izh2p, Izh3p, and Izh4p are involved in the regulation of intracellular zinc levels. Izh2p has further been reported to play a role in lipid and phosphate metabolism [53, 54], and to function as a receptor for the plant defense protein osmotin which induces programmed cell death in yeast .
In the genomes of filamentous fungi such as N. crassa, A. nidulans, F. graminearum, and M. grisea two to three PAQR-type proteins are encoded and have been designated as class VIII of fungal GPCRs [1, 2]. Our mining of the genomes of T. virens and T. atroviride revealed the presence of six and seven PAQR members (Table 1, Figure 1), respectively, all of which bear the hemolysin III motif (pfam03006, HlyIII) and which face five members identified in T. reesei[38, 39]. Phylogenetic analysis showed the Trichoderma orthologues Triat136196, Trive180426, Trire56426 in a clade together with yeast Izh3 (Figure 2). Izh3 possesses a long N-terminal tail with unknown function as a distinctive characteristic . Similar extracellular N-terminal extensions of ~280 amino acids were found in the Trichoderma Izh3-like proteins Triat136196, Trive180426 and Trire56426. It is worth mentioning that some of the Trichoderma class VIII members do not share the typical GPCR topology but have an extracellular C-terminus and the N-terminal domain within the cytoplasm. Triat210209, Triat46847, Triat142943, Trire82246, Trive92622 are in the same, although not well supported, cluster with the human adinopectin receptors adipor1-human and adipor2-human, which share the same topology .
To analyze whether the class VIII genes identified in the Trichoderma genomes are actually transcribed, their expression was assessed by RT-qPCR. Respective transcripts were detected for all five and six genes of T. reesei and T. virens, respectively, as well as for six of the seven genes identified in the T. atroviride genome (Figure 3). Triat46847 was not transcribed under the growth condition tested (PDA). Analysis of mRNA levels after co-cultivation of Trichoderma with Rhizoctonia solani revealed a significantly enhanced expression of Trive160502 (pâ€‰=â€‰0.000) and Trive180426 (pâ€‰=â€‰0.031) in T. virens, Triat152366 (pâ€‰=â€‰0.027) and Triat210209 (pâ€‰=â€‰0.000) in T. atroviride, and Trire56426 (pâ€‰=â€‰0.000) in T. reesei upon contact with the host fungus (Figure 3). On the other hand, expression of Triat142946 (pâ€‰=â€‰0.000), Triat136196 (pâ€‰=â€‰0.000) in T. atroviride, Trive92622 (pâ€‰=â€‰0.000), Trive47976 (pâ€‰=â€‰0.000), Trive30459 (pâ€‰=â€‰0.034) in T. virens, and Trire70139 (pâ€‰=â€‰0.032), Trire119819 (pâ€‰=â€‰0.000) in T. reesei was significantly decreased in the presence of R. solani compared to the corresponding controls. Transcript levels of Triat290043 (pâ€‰=â€‰0.971), Triat142943 (pâ€‰=â€‰0.093), and Trire82246 (pâ€‰=â€‰0.102) were unaffected by the presence of R. solani. Again no transcript could be detected for Triat46847. Expression of Triat46847 was further assessed on both plates and in liquid minimal and full media and under different developmental stages (vegetative growth, conidiation) of the fungus. No transcript could be detected under all the conditions tested (data not shown).
Analysis of the location of the seven PAQR-encoding genes in the genome of T. atroviride revealed that three of them (Triat142946, Triat142943, Triat46847) are in close vicinity on scaffold 19 (Figure 4). This is similar in T. virens and T. reesei for the orthologues of Triat142946 and Triat142943 suggesting the possibility that the third T. atroviride gene (Triat46847), which was found not to be expressed under any of the conditions tested, may have resulted from gene duplication with subsequent inactivation.
The finding that the genes located in the genomes of both T. atroviride and T. virens between the orthologous receptor triplets Triat142946/Trive160502/Trire70139 and Triat142943/Trive92622/Trire82246 have been lost in T. reesei (Figure 4) is consistent with a reported paralogous gene expansion in T. atroviride and T. virens compared to T. reesei and other non-mycoparasitic fungi .
After the class of PTH11-like receptors, the PAQR family is the second largest GPCR class in Trichdoderma. The expansion of the PAQR family especially in T. atroviride and T. virens together with the fact that S. cerevisiae Izh2 was found to regulate fungal development in response to plant osmotin , make these receptors interesting candidates for an involvement in interspecies communication between Trichoderma and other (host) fungi and/or plants. The importance of fungal class VIII GPCRs in environmental sensing is further supported by the recent characterization of a PAQR family member of the fungus Sporothrix schenkii. SsPAQR1 was found to respond to the steroid hormone progesterone by signaling via the GÎ± subunit SSG-2 .
Trichoderma members of classes IX to XII of fungal GPCRs
A 7-transmembrane protein with a bacteriorhodopsin domain is encoded in the genome of T. atroviride. Triat210598 is orthologous to N. crassa NOP-1 and ORP-1 and A. nidulans NopA (Additional file 1). Interestingly, Triat210598 has no homologs in T. reesei and T. virens. Due to the finding that Triat210598 is located in a non-syntenic genome region it has been suggested that T. reesei and T. virens have lost this gene during evolution . This hypothesis is in agreement with recent results showing that T. reesei and T. virens are derived relative to T. atroviride, the latter resembling the more ancient state of Trichoderma.
Classes X, XI, and XII of fungal GPCRs have recently been defined in Verticillium spp. . Similar to Verticillium and other filamentous fungi such as A. nidulans, M. grisea, N. crassa, and F. graminearum, one putative PTM1-like GPCR was identified in the two mycoparasites T. atroviride and T. virens as well as the saprophyte T. reesei. Consistent with the presence of a Lung_7-TM_R domain (pfam06814) and similarity to the putative tumor necrosis factor receptor-like GPCR PTM1 of S. cerevisiae, the respective Trichoderma proteins were designated as class X members (Table 1).
One putative member related to human GPR89A was identified in the genome of each of the three Trichoderma species (Table 1). The Trichoderma proteins showed the typical structure previously described for receptors of class XI with 9 transmembrane regions and a large third cytoplasmic loop , and contain a ABA_GPCR (pfam12430; abscisic acid G protein-coupled receptor) domain.
Putative fungal receptors with similarity to family C-like GPCRs (metabotropic glutamate/pheromone receptors) have previously been defined as class XII . Similar to other filamentous ascomycetes, one putative GPCR grouping to this class was identified in each of the three Trichoderma species. Whereas the respective proteins of both T. atroviride and T. reesei exhibit the typical structure with 7 transmembrane domains and the long C-terminal tail, the T. virens homologue (Trive179509) only exhibits 6 transmembrane regions.
PTH11-Related proteins ofTrichoderma
The PTH11 receptor was first identified in M. grisea, where it is required for host surface recognition and pathogenicity . PTH11 has an extracellular amino-terminal CFEM domain followed by seven transmembrane regions and PTH11-related proteins are restricted to fungi belonging to the subphylum Pezizomycotina .
In both the mycoparasitic Trichoderma species as well as T. reesei[38, 39], the number of identified PTH11-like proteins was higher than in the saprophyte N. crassa (25 members) but lower than in the plant pathogens M. grisea (61 members) and F. graminearum (106 members) [2, 14]. Similar to the above mentioned fungi, only a subset of the identified Trichoderma proteins contained the fungal-specific cysteine-rich CFEM (pfam05730) domain (Figure 5, Additional file 2), which is characteristically present in the extracellular region of some membrane proteins with proposed roles in fungal pathogenicity. Compared to T. atroviride (38 members) and T. reesei (35 members), we found a marked expansion of PTH11-related proteins in T. virens (52 members).
Additional putative GPCRs of Trichoderma which are beyond the existing classification system of fungal GPCRs (class XIII)
Recently, a putative GPCR of Phytophtora sojae (GPR11) controlling zoospore development and virulence of P. sojae to soybean has been described . Performing a BLASTP search with GPR11 as a query against the proteomes of T. atroviride, T. virens, T. reesei, and those of N. crassa, M. grisea, and A. fumigatus revealed respective orthologues in all fungi tested. Whereas in T. atroviride three proteins were identified (Table 1), T. reesei and T. virens as well as the other ascomycetes possess two members each. All putative Trichoderma GPCRs identified this way have a DUF300 domain (domain of unknown function, pfam03619). Such a domain is also present in e.g. the class A GPCRs Cand9 and Cand10 of Arabidopsis thaliana and P. sojae GPR11. Topological analysis of the Trichoderma proteins revealed a heptahelical topology with three N-terminal transmembrane regions, a long second cytoplasmic loop followed by four transmembrane regions and a long intracellular loop at the C-terminus. As these putative GPCRs represented a separate clade in the phylogenetic analysis (Figure 1), they were assigned to a new class (class XIII, Table 1) thereby extending the classification system of fungal GPCRs to 14 classes.
A thorough examination of the genomes of the two mycoparasites T. atroviride and T. virens and the saprophyte T. reesei for putative GPCRs revealed for most classes a high conservation of their number and structure within this genus. On the other hand, remarkable differences in individual classes were found among the three Trichoderma species and among Trichoderma and other filamentous fungi. Whereas for class I to VII members, orthologous triplets with similar length and sequence are present in the genomes of the three Trichoderma species and their number is also similar to other fungi, the PAQR family has expanded especially in T. atroviride. Considering the identification of members of classes X, XI, and XII and proteins similar to the P. sojae GPR11 receptor in Trichoderma, the presented 14 classes now define the most comprehensive classification system for GPCR-like proteins of fungi. The huge diversity of GPCRs in Trichoderma spp. and especially in the mycoparasites is likely to reflect the capability of these fungi to establish various ecological niches and interactions with other organisms.
It is worth mentioning that with the exception of few members, the proteins identified as putative GPCRs in this study have only been characterized in silico. Taking into account that only three Î±, one Î² and one Î³ subunit of heterotrimeric G proteins are encoded in the Trichoderma genomes which face more than 55 GPCRs, studying the signaling output and identifying the respective intracellular interaction partners of those receptors will provide interesting insights on how these fungi adapt to their different lifestyles.
Identification of GPCR-encoding genes of Trichoderma atroviride andTrichoderma virens
Version 2 of the T. atroviride genome database  comprises 11,863 gene models on 29 scaffolds; version 2 of the T. virens genomic sequence  comprises 12,427 gene models on 93 scaffolds. For the homology-based search of GPCR-like proteins from T. atroviride and T. virens, the genomic sequences and deduced proteomes of the following fungi were used: Trichoderma reeseiAspergillus nidulans, Aspergillus fumigatus, Aspergillus oryzae, Neurospora crassa, Magnaporthe grisea, Podospora anserine, Chaetomium globosum, Fusarium graminearum, and Nectria haematococca. An e-value limit of 1e-09 was applied.
To identify putative GPCRs within the T. atroviride and T. virens proteomes that lack significant sequence similarity to known GPCR-like proteins and therefore may escape detection by homology search, a more sensitive database searching using hidden Markov models (HMM) was performed using the program HMMER (http://hmmer.janelia.org/) .
All obtained predicted proteins were analyzed with the TMHMM, ConPred II and HMMTOP algorithms [70–72] to test for the typical 7-transmembrane domain topology. For those few proteins exhibiting less than seven transmembrane domains, the respective encoding gene and flanking regions were retrieved from the genome database and examined manually. Wrongly predicted intron-exon boundaries were mainly found and manually corrected resulting in the detection of the missing transmembrane domains.
Protein alignments and phylogenetic analysis
The classification system of Lafon et al. , which classifies fungal GPCRs into nine classes according to their sequence similarity, was applied to all detected putative GPCRs of Trichoderma. In addition, members of the three additional classes identified in Verticillium spp. , and the GPR11 protein of Phytophtora sojae were used to identify and classify respective members of T. atroviride, T. virens and T. reesei. Multiple sequence alignments of the identified putative GPCR-like proteins and phylogenetic trees with a neighbor-joining approach were generated using ClustalX . A bootstrap with 1000 repetitions was included.
Cultivations and RT-qPCR analysis
T. atroviride strain P1 (ATCC 74058; teleomorph Hypocrea atroviridis), T. virens strain IMI 206040 (teleomorph Hypocrea virens), and T. reesei strain QM6a (ATCC13631; teleomorph Hypocrea jecorina) were used in this study. The fungi were cultivated at 28Â°C on either complete medium (PDA, PDB) or minimal medium (MM, containing [g/l]: MgSO4â€‰Â·â€‰7H2O 1, KH2PO4 10, (NH4)2SO4 6, tri-sodium citrate 3, FeSO4â€‰Â·â€‰7H2O 0.005, ZnSO4â€‰Â·â€‰2H2O 0.0014, CoCl2â€‰Â·â€‰6H2O 0.002, MnSO4â€‰Â·â€‰6H2O 0.0017, glucose 10) on plates and in liquid culture, respectively. Plate confrontation assays were performed by cultivating Trichoderma together with Rhizoctonia solani on PDA plates covered with a cellophane membrane at 28Â°C. After direct contact between the two fungi, mycelium of Trichoderma was harvested from the confrontation zone. For RNA isolation, 30 mg fungal mycelium was grinded in liquid nitrogen and RNA isolated using the peqGOLD TriFast Solution (PeqLab, Erlangen, Germany) according to the manufacturerÂ´s instructions.
For cDNA synthesis the Revert Aid H Minus First Strand cDNA Synthesis Kit (Fermentas, Vilnius, Lithuania) was used according to the manufacturerÂ´s instructions with a combination of an oligo(dT)18 and a random hexamer primer. The sequences for the respective primer pairs for cDNA amplification of the reference gene sar1 and the genes encoding the putative receptors of class VIII identified in the Trichoderma genomes are given in Additional file 3. Transcript quantification was performed with the following PCR program (initial denaturation for 120 s at 95Â°C, 50 cycles with 95Â°C for 20 s, 60Â°C for 20 s and 72Â°C for 20 s) on an Eppendorf (Hamburg, Germany) realplex2S Mastercycler using the IQ SYBR Green Supermix (Bio-Rad, Hercules, CA) and 25 Î¼l assays with standard MgCl2 concentration (3 mM) and with final primer concentrations of 100 nM each. All assays were carried out in 96-well plates covered with optical tape. PCR efficiency was determined from a single tube reaction set-up as described  and expression ratio was calculated according to Pfaffl . All samples were analyzed in three independent experiments with three replicates in each run. Statistical analysis was done by relative expression analysis with REST software using the Pair Wise Fixed Reallocation Randomisation Test .
SZ conceived the study, drafted the manuscript, and performed in silico analyses together with MO. SG contributed to gene identifications and performed the cultivations and RT-qPCR experiments. All authors read and approved the final manuscript.
Lafon A, Han KH, Seo JA, Yu JH, d'Enfert C: G-protein and cAMP-mediated signaling in aspergilli: a genomic perspective. Fungal Genet Biol. 2006, 43: 490-502. 10.1016/j.fgb.2006.02.001.
Li L, Wright SJ, Krystofova S, Park G, Borkovich KA: Heterotrimeric G Protein Signaling in Filamentous Fungi. Annu Rev Microbiol. 2007, 61: 423-452. 10.1146/annurev.micro.61.080706.093432.
Xue C, Hsueh YP, Heitman J: Magnificent seven: roles of G protein coupled receptors in extracellular sensing in fungi. FEMS Microbiol Rev. 2008, 32: 1010-1032. 10.1111/j.1574-6976.2008.00131.x.
Kroeze WK, Sheffler DJ, Roth BL: G-protein-coupled receptors at a glance. J Cell Sci. 2003, 116: 4867-10.1242/jcs.00902.
Dohlman H, Thorner J, Caron M, Lefkowitz R: Model systems for the study of seven-transmembrane-segment receptors. Annu Rev Bbiochem. 1991, 60: 653-688. 10.1146/annurev.bi.60.070191.003253.
Oldham WM: Structural basis of function in heterotrimeric G proteins. Quaterly Rev Biophys. 2006, 39: 117-166. 10.1017/S0033583506004306.
Gutkind JS: The pathways connecting G protein-coupled receptors to the nucleus through divergent mitogen-activated protein kinase cascades. J Biol Chem. 1839, 1998: 273.
Neer EJ: Heterotrimeric G proteins: Organizers of transmembrane signals. Cell. 1995, 80: 249-257. 10.1016/0092-8674(95)90407-7.
Oldham WM, Hamm HE: Heterotrimeric G protein activation by G-protein-coupled receptors. Nature Rev Mol Cell Biol. 2008, 9: 60-71. 10.1038/nrm2299.
Bhattacharya M, Babwah A, Ferguson S: Small GTP-binding protein-coupled receptors. Biochem Soc Trans. 2004, 32: 1040-1044. 10.1042/BST0321040.
Hall RA, Premont RT, Lefkowitz RJ: Heptahelical receptor signaling: beyond the G protein paradigm. J Cell Biol. 1999, 145: 927-10.1083/jcb.145.5.927.
Mitchell R, McCulloch D, Lutz E, Johnson M, MacKenzie C, Fennell M, Fink G, Zhou W, Sealfon SC: Rhodopsin-family receptors associate with small G proteins to activate phospholipase D. Nature. 1998, 392: 411-414. 10.1038/32937.
Xiao K, Sun J, Kim J, Rajagopal S, Zhai B, VillÃ©n J, Haas W, Kovacs JJ, Shukla AK, Hara MR, Hernandez M, Lachmann A, Zhao S, Lin Y, Cheng Y, Mizuno K, Ma'ayan A, Gygi SP, Lefkowitz RJ: Global phosphorylation analysis of Î²-arrestinâ€“mediated signaling downstream of a seven transmembrane receptor (7TMR). Proc Nat Acad Sci USA. 2010, 107: 15299-15304. 10.1073/pnas.1008461107.
Kulkarni RD, Thon MR, Pan H, Dean RA: Novel G-protein-coupled receptor-like proteins in the plant pathogenic fungus Magnaporthe grisea. Genome Biol. 2005, 6: R24-10.1186/gb-2005-6-3-r24.
Blumer KJ, Reneke JE, Courchesne WE, Thorner J: Functional domains of a peptide hormone receptor: the alpha-factor receptor (STE2 gene product) of the yeast Saccharomyces cerevisiae. Cold Spring Harb Symp Quant Biol. 1998, 53 (Pt 2): 591-603.
Chang YC, Miller GF, Kwon-Chung K: Importance of a developmentally regulated pheromone receptor of Cryptococcus neoformans for virulence. Infect Immun. 2003, 71: 4953-10.1128/IAI.71.9.4953-4960.2003.
Hagen DC, McCaffrey G, Sprague GF: Evidence the yeast STE3 gene encodes a receptor for the peptide pheromone a factor: gene sequence and implications for the structure of the presumed receptor. Proc Nat Acad Sci USA. 1986, 83: 1418-10.1073/pnas.83.5.1418.
Hsueh YP, Xue C, Heitman J: A constitutively active GPCR governs morphogenic transitions in Cryptococcus neoformans. EMBO J. 2009, 28: 1220-1233. 10.1038/emboj.2009.68.
Kim H, Borkovich KA: A pheromone receptor gene, pre 1, is essential for mating type specific directional growth and fusion of trichogynes and female fertility in Neurospora crassa. Mol Microbiol. 2004, 52: 1781-1798. 10.1111/j.1365-2958.2004.04096.x.
Krystofova S, Borkovich KA: The predicted G-protein-coupled receptor GPR-1 is required for female sexual development in the multicellular fungus Neurospora crassa. Eukaryot Cell. 2006, 5: 1503-10.1128/EC.00124-06.
Li L, Borkovich KA: GPR-4 is a predicted G-protein-coupled receptor required for carbon source-dependent asexual growth and development in Neurospora crassa. Eukaryot Cell. 2006, 5: 1287-10.1128/EC.00109-06.
Miwa T, Takagi Y, Shinozaki M, Yun CW, Schell WA, Perfect JR, Kumagai H, Tamaki H: Gpr1, a putative G-protein-coupled receptor, regulates morphogenesis and hypha formation in the pathogenic fungus Candida albicans. Eukaryot Cell. 2004, 3: 919-10.1128/EC.3.4.919-931.2004.
Seibel C, Tisch D, Kubicek CP, Schmoll M: The pheromone receptors for communication and mating in Hypocrea jecorina. Fungal Genet Biol. 2012, 49 (10): 814-824. 10.1016/j.fgb.2012.07.004.
Tanaka K, Davey J, Imai Y, Yamamoto M: Schizosaccharomyces pombe map3+ encodes the putative M-factor receptor. Mol Cell Biol. 1993, 13: 80.
Xue C, Bahn YS, Cox GM, Heitman J: G protein-coupled receptor Gpr4 senses amino acids and activates the cAMP-PKA pathway in Cryptococcus neoformans. Mol Biol Cell. 2006, 17: 667.
Yun CW, Tamaki H, Nakayama R, Yamamoto K, Kumagai H: G-protein coupled receptor from yeast Saccharomyces cerevisiae. Biochem Biophys Res Commun. 1997, 240: 287-292. 10.1006/bbrc.1997.7649.
Druzhinina IS, Seidl-Seiboth V, Herrera-Estrella A, Horwitz BA, Kenerley CM, Monte E, Mukherjee PK, Zeilinger S, Grigoriev IV, Kubicek CP: Trichoderma: the genomics of opportunistic success. Nat Rev Microbiol. 2011, 16: 749-759.
Omann M, Zeilinger S: How a mycoparasite employs g-protein signaling: using the example ofTrichoderma. Journal of Signal Transduction. 2010, 2010: 123126.
Reithner B, Brunner K, Schuhmacher R, Peissl I, Seidl V, Krska R, Zeilinger S: The G protein alpha subuniz Tga1 of Trichoderma atroviride is involved in chitinase formation and differential production of antifungal metabolites. Fungal Genet Biol. 2005, 42 (9): 749-760. 10.1016/j.fgb.2005.04.009.
Rocha-Ramirez V, Omero C, Chet I, Horwitz BA, Herrera-Estrella A: Trichoderma atroviride G protein alpha subunit gene tga1 is involved in mycoparasitic coiling and conidiation. Eukaryot Cell. 2002, 1 (4): 594-605. 10.1128/EC.1.4.594-605.2002.
Zeilinger S, Reithner B, Scala V, Peissl I, Lorito M, Mach RL: Signal transduction by Tga3, a novel G protein alpha subunit of Trichoderma atroviride. Appl Environ Microbiol. 2005, 71: 1591-10.1128/AEM.71.3.1591-1597.2005.
Mukherjee PK, Latha J, Hadar R, Horwitz BA: Role of two G protein alpha subunits, tgaA and TgaB, in the antagonism of plant pathogens by Trichoderma virens. Appl Environ Microbiol. 2004, 70 (1): 542-549. 10.1128/AEM.70.1.542-549.2004.
Schmoll M, Esquivel-Naranjo EU, Herrera-Estrella A: Trichoderma in the light of day-physiology and development. Fungal Genet Biol. 2010, 47: 909-916. 10.1016/j.fgb.2010.04.010.
Tisch D, Kubicek CP, Schmoll M: The phosducin-like protein PhLP1 impacts regulation of glycoside hydrolases and light response in Trichoderma reesei. BMC Genomics. 2011, 12: 613-10.1186/1471-2164-12-613.
Wang Y, Li A, Wang X, Zhang X, Zhao W, Dou D, Zheng X: GPR11, a putative seven-transmembrane G protein-coupled receptor, controls zoospore development and virulence of Phytophthora sojae. Eukaryot Cell. 2010, 9: 242-10.1128/EC.00265-09.
Zheng H, Zhou L, Dou T, Han X, Cai Y, Zhan X, Tang C, Huang J, Wu Q: Genome-wide prediction of G protein-coupled receptors in Verticillium spp. Fungal Biol. 2010, 114: 359-368. 10.1016/j.funbio.2010.02.008.
DeZwaan TM, Carroll AM, Valent B, Sweigard JA: Magnaporthe grisea Pth11p is a novel plasma membrane protein that mediates appressorium differentiation in response to inductive substrate cues. The Plant Cell Online. 2013, 1999: 11.
Brunner K, Omann M, Pucher ME, Delic M, Lehner SM, Domnanich P, Kratochwill K, Druzhinina I, Denk D, Zeilinger S: Trichoderma G protein-coupled receptors: functional characterisation of a cAMP receptor-like protein from Trichoderma atroviride. Curr Genet. 2008, 54: 283-299. 10.1007/s00294-008-0217-7.
Schmoll M: The information highways of a biotechnological workhorseâ€“signal transduction in Hypocrea jecorina. BMC Genomics. 2008, 9: 430-10.1186/1471-2164-9-430.
Kubicek CP, Herrera-Estrella A, Seidl-Seiboth V, Martinez DA, Druzhinina IS, Thon M, Zeilinger S, Casas-Flores S, Horwitz BA, Mukherjee PK, Mukherjee M, Kredics L, Alcaraz LD, Aerts A, Antal Z, Atanasova L, Cervantes-Badillo MG, Challacombe J, Chertkov O, McCluskey K, Coulpier F, Deshpande N, Von DÃ¶hren H, Ebbole DJ, Esquivel-Naranjo EU, Fekete E, Flipphi M, Glaser F, GÃ³mez-RodrÃguez EY, Gruber S, Han C, Henrissat B, Hermosa R, HernÃ¡ndez-OÃ±ate M, Karaffa L, Kosti I, Le Crom S, Lindquist E, Lucas S, LÃ¼beck M, LÃ¼beck PS, Margeot A, Metz B, Misra M, Nevalainen H, Omann M, Packer N, Perrone G, Uresti-Rivera EE, Salamov A, Schmoll M, Seiboth B, Shapiro H, Sukno S, Tamayo-Ramos JA, Tisch D, Wiest A, Wilkinson HH, Zhang M, Coutinho PM, Kenerley CM, Monte E, Baker SE, Grigoriev IV: Comparative genome sequence analysis underscores mycoparasitism as the ancestral life style of Trichoderma. Genome Biol. 2011, 12: R40-10.1186/gb-2011-12-4-r40.
Chaverri P, Castlebury LA, Samuels GJ, Geiser DM: Multilocus phylogenetic structure within the Trichoderma harzianum/Hypocrea lixii complex. Mol Phyl Evol. 2003, 27: 302-313. 10.1016/S1055-7903(02)00400-1.
Dodd SL, Lieckfeldt E, Samuels GJ: Hypocrea atroviridis sp. nov., the teleomorph of Trichoderma atroviride. Mycologia. 2003, 95: 27-40. 10.2307/3761959.
Lemaire K, Van de Velde S, Van Dijck P, Thevelein JM: Glucose and sucrose act as agonist and mannose as antagonist ligands of the G protein-coupled receptor Gpr1 in the yeast Saccharomyces cerevisiae. Mol Cell. 2004, 16: 293-299. 10.1016/j.molcel.2004.10.004.
Lorenz MC, Pan X, Harashima T, Cardenas ME, Xue Y, Hirsch JP, Heitman J: The G protein-coupled receptor Gpr1 is a nutrient sensor that regulates pseudohyphal differentiation in Saccharomyces cerevisiae. Genetics. 2000, 154: 609.
Gehrke A, Heinekamp T, Jacobsen ID, Brakhage AA: Heptahelical receptors GprC and GprD of Aspergillus fumigatus are essential regulators of colony growth, hyphal morphogenesis, and virulence. Appl Environ Microbiol. 2010, 76: 3989-10.1128/AEM.00052-10.
Han KH, Seo JA, Yu JH: A putative G protein coupled receptor negatively controls sexual development in Aspergillus nidulans. Mol Microbiol. 2004, 51: 1333-1345. 10.1111/j.1365-2958.2003.03940.x.
Affeldt KJ, Brodhagen M, Keller NP: Aspergillus oxylipin signaling and quorum sensing pathways depend on G protein-coupled receptors. Toxins. 2012, 4: 695-717. 10.3390/toxins4090695.
Chung KS, Won M, Lee SB, Jang YJ, Hoe KL, Kim DU, Lee JW, Kim KW, Yoo H: Isolation of a Novel Gene fromSchizosaccharomyces pombe:stm1+ Encoding a Seven-transmembrane Loop Protein That May Couple with the Heterotrimeric G 2 Protein, Gpa2. J Biol Chem. 2001, 276: 40190.
Seidl V, Song L, Lindquist E, Gruber S, Koptchinskiy A, Zeilinger S, Schmoll M, MartÃnez P, Sun J, Grigoriev I, Herrera-Estrella A, Baker SE, Kubicek CP: Transcriptomic response of the mycoparasitic fungus Trichoderma atroviride to the presence of a fungal prey. BMC Genomics. 2009, 10: 567-10.1186/1471-2164-10-567.
Omann MR, Lehner S, Escobar RodrÃguez C, Brunner K, Zeilinger S: The seven-transmembrane receptor Gpr1 governs processes relevant for the antagonistic interaction of Trichoderma atroviride with its host. Microbiology. 2012, 158: 107-118. 10.1099/mic.0.052035-0.
Chen JG, Willard FS, Huang J, Liang J, Chasse SA, Jones AM, Siderovski DP: A seven-transmembrane RGS protein that modulates plant cell proliferation. Science. 2003, 301: 1728-10.1126/science.1087790.
Tang YT, Hu T, Arterburn M, Boyle B, Bright JM, Emtage PC, Funk WD: PAQR proteins: a novel membrane receptor family defined by an ancient 7-transmembrane pass motif. J Mol Evol. 2005, 61 (3): 372-380. 10.1007/s00239-004-0375-2.
Karpichev IV, Cornivelli L, Small GM: Multiple regulatory roles of a novel Saccharomyces cerevisiae protein, encoded by YOL002c, in lipid and phosphate metabolism. J Biol Chem. 2002, 277: 19609-10.1074/jbc.M202045200.
Lyons TJ, Villa NY, Regalla LM, Kupchak BR, Vagstad A, Eide DJ: Metalloregulation of yeast membrane steroid receptor homologs. Proc Nat Acad Sc USA. 2004, 101: 5506-10.1073/pnas.0306324101.
Narasimhan ML, Coca MA, Jin J, Yamauchi T, Ito Y, Kadowaki T, Kim KK, Pardo JM, Damsz B, Hasegawa PM, Yun DJ, Bressan RA: Osmotin is a homolog of mammalian adiponectin and controls apoptosis in yeast through a homolog of mammalian adiponectin receptor. Mol Cell. 2005, 17: 171-180. 10.1016/j.molcel.2004.11.050.
Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17: 540-10.1093/oxfordjournals.molbev.a026334.
The Trichoderma atroviride genome database. http://genome.jgi-psf.org/Triat2/Triat2.home.html,
The Trichoderma virens genome database. http://genome.jgi-psf.org/TriviGv29_8_2/TriviGv29_8_2.home.html,
The Aspergillus comparative database. http://www.broadinstitute.org/annotation/genome/aspergillus_group/MultiHome.html,
The Trichoderma reesei genome database. http://genome.jgi-psf.org/Trire2/Trire2.home.html,
Gookin TE, Kim J, Assmann SM: Whole proteome identification of plant candidate G protein-coupled receptors in Arabidopsis, rice, and poplar: computational prediction and in-vivo coupling. Genome Biol. 2008, 9 (7): R120-10.1186/gb-2008-9-7-r120.
Gonzalez-Velazquez W, Gonzalez-Mendez R, Rodriguez-DelValle N: Characterization and ligand identification of a membrane progesterone receptor in fungi: existence of a novel PAQR in Sporothrix schenkii. BMC Microbiol. 2012, 12: 194-10.1186/1471-2180-12-194.
The Neurospora crassa genome database. http://www.broad.mit.edu/annotation/genome/neurospora/Home.html,
The Magnaporthe grisea genome database. http://www.broad.mit.edu/annotation/fungi/magnaporthe,
The Podospora anserina genome database. http://podospora.igmors.u-psud.fr,
The Chaetomium globosum genome database. http://www.broad.mit.edu/annotation/genome/chaetomium_globosum/Home.html,
The Fusarium graminearum genome database. http://mips.gsf.de/genre/proj/fusarium,
The Nectria haematococca genome database. http://genome.jgi-psf.org/Necha2/Necha2.home.html,
Durbin R, Eddy S, Krogh A, Mitchison G: Biological sequence analysis: probabilistic models of proteins and nucleic acids. 1998, Cambridge: Cambridge University Press
Arai M, Mitsuke H, Ikeda M, Xia JX, Kikuchi T, Satake M, Shimizu T: ConPred II: a consensus prediction method for obtaining transmembrane topology models with high reliability. Nucleic Acids Res. 2004, 32: W390-10.1093/nar/gkh380.
Krogh A, Larsson BÃˆ, Von Heijne G, Sonnhammer ELL: Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J Mol Biol. 2001, 305: 567-580. 10.1006/jmbi.2000.4315.
Tusnady GE, Simon I: The HMMTOP transmembrane topology prediction server. Bioinformatics. 2001, 17: 849-10.1093/bioinformatics/17.9.849.
Larkin M, Blackshields G, Brown NP, Chenna R, McGettigan PA, MCWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-10.1093/bioinformatics/btm404.
Tichopad A, Dilger M, Schwarz G, Pfaffl MW: Standardized determination of real time PCR efficiency from a single reaction set up. Nucleic Acids Res. 2003, 31: e122-10.1093/nar/gng122.
Pfaffl MW: A new mathematical model for relative quantification in real-time RTâ€“PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.
Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30: E36-10.1093/nar/30.9.e36.
This work was supported by the Austrian Science Fund FWF (grant V139-B20) and the Vienna Science and Technology Fund WWTF (grant LS09-036).
The authors declare that they have no competing interests.