- Short Report
- Open Access
Adult Drosophila melanogaster evolved for antibacterial defense invest in infection-induced expression of both humoral and cellular immunity genes
© McGraw et al; licensee BioMed Central Ltd. 2011
- Received: 16 June 2011
- Accepted: 23 August 2011
- Published: 23 August 2011
While the transcription of innate immunity genes in response to bacterial infection has been well-characterised in the Drosophila model, we recently demonstrated the capacity for such transcription to evolve in flies selected for improved antibacterial defense. Here we use this experimental system to examine how insects invest in constitutive versus infection-induced transcription of immunity genes. These two strategies carry with them different consequences with respect to energetic and pleiotropic costs and may be more or less effective in improving defense depending on whether the genes contribute to humoral or cellular aspects of immunity.
Contrary to expectation we show that selection preferentially increased the infection-induced expression of both cellular and humoral immunity genes. Given their functional roles, infection induced increases in expression were expected for the humoral genes, while increases in constitutive expression were expected for the cellular genes. We also report a restricted ability to improve transcription of immunity genes that is on the order of 2-3 fold regardless of total transcription level of the gene.
The evolved increases in infection-induced expression of the cellular genes may result from specific cross talk with humoral pathways or from generalised strategies for enhancing immunity gene transcription. A failure to see improvements in constitutive expression of the cellular genes suggests either that increases might come at too great a cost or that patterns of expression in adults are decoupled from the larval phase where increases would be most effective. The similarity in fold change increase across all immunity genes may suggest a shared mechanism for the evolution of increased transcription in small, discrete units such as duplication of cis-regulatory elements.
- Constitutive Expression
- Immunity Gene
- Cellular Gene
- Antibacterial Defense
Recently, using selection experiments, we examined the ability of Drosophila melanogaster to evolve in response to systemic infection by the opportunistic bacterial pathogen, Pseudomonas aeruginosa . The evolved flies demonstrated an improved defense as measured by an increase in fly survival in response to infection from 15 to 70%. The genetic basis of the evolved response was examined using microarrays and was characterised by increased expression for many immunity related genes. The evolved defense was costly, as evidenced by decreases in egg viability and longevity and the rapid loss of the increased survival when selection was removed. While the original study demonstrated the involvement of both cellular and humoral aspects of insect immunity in the evolved defense, the design did not allow for the partitioning of the total transcriptional change into its constitutive versus infection-induced components. These two avenues for transcriptional change should have different energetic or pleiotropic costs and  may be more or less effective in improving immunity depending on the functional role of the gene. The original study also only focused on genetic changes in males and so did not determine if females were arriving at their improved defense by the same genetic mechanism.
The humoral immune response in insects is characterized by cascades of genes whose expression is controlled by a range of recognition proteins that bind pathogen cell wall components and end with the production of a range of secreted host proteins that target bacteria for destruction . Exposure to infection can increase the transcription of these genes substantially, ie from 10-100 fold for some of the genes encoding antimicrobial peptides . We hypothesized that selection for improved defense would lead to further increases in the transcription of these genes upon exposure to infection. This would represent a targeted strategy, producing these molecules only when needed at a point in the infection process where they normally act against invading bacteria. The transient and inducible nature of this aspect of immunity might suggest, however, that flies are already be operating this response at or near maximal transcriptional capacity, potentially limiting the size of gains we might see. Increasing the constitutive expression of these humoral genes, though more costly, could also be beneficial. A greater number of circulating antimicrobial peptides prior to exposure could confer an enhanced ability to contain or curtail infection and recent work has suggested that fly genotypes with greater constitutive expression of these genes are better able to control infection .
The genes that define the cellular response of immunity underpin processes like encapsulation, lytic ability and phagocytosis. The cellular response has been shown to be highly effective in early infection, including the recent suggestion that the inducible humoral response may exist only for "clean up" while the cellular response serves as the frontline of antibacterial defense . The expression of these genes early in development can determine subsequent immunological traits in adults like the number and quality of circulating hemocytes [3, 7]. Many of these genes may also have pleiotropic roles to play with respect to early development [8, 9]. In our study, flies with an evolved defense also developed more rapidly, which may have resulted from such pleiotropy. Several key studies in recent years have highlighted the cost of investing specifically in cellular immunity with respect to various trade-offs in life history traits [6, 10–12]. In contrast to the secreted products of the humoral response, we would expect that only increases in the constitutive expression of these genes would result in improved antibacterial defense. Any infection induced increases in expression should theoretically be "too late" or too slow convert to infection fighting ability . Given the significant improvement in survival of selected flies, the cost of this improvement, the concomitant increase in developmental rate and the likely point of action of these proteins early in development , we expected to see evolution of increased constitutive expression for the cellular immunity genes.
Here we compare the expression patterns of immunity genes between evolved and unevolved flies in both the constitutive and infection-induced states. In so doing, we are able to partition the relative investment into constitutive versus infection-induced transcription of genes representing both humoral and cellular aspects of immunity due to selection and understand how investment in these different strategies may form the basis of fly immunity.
Fly and bacterial culture and the selection regime are as previously reported  but are described briefly below. Female flies generated from the previous selection experiment are then screened here for immunity gene transcriptional profiles as an addition to the original study.
Fly and Bacterial Culture
Brisbane (BNE) base stock was founded from 26 females D. melanogaster caught around the University of Queensland St Lucia campus in August 2006. The flies were initially treated with 0.5% penicillin and streptomycin in the diet for one generation  and then passaged without antibiotic on standard cornmeal diet supplemented with excess yeast for more than 10 generations before the start of the selection experiment. P. aeruginosa PA01 was cultured in LB medium supplemented with 100 mg ml-1 ampicillin at 37°C on a rotory shaker . For infection, the concentration of an overnight bacterial culture was adjusted to an optical density (OD) of 0.5 ± 0.05 measured spectrometrically at 600 nm. The culture was then diluted 100 fold using sterile LB. This OD was determined at the start of the selection experiments to achieve a fly population kill rate of 80-90%. To determine the number of bacterial cells used to infect each fly, individual flies were ground with a pestle in an Eppendorf tube with 200 μl of 10 mM MgSO4 and spread on LB plates supplemented with 100 mg ml-1 ampicillin immediately after infection . The average infective dose per fly was 3109 ± 635 bacterial cells.
Selection Regime and Measurement of Defense
The base stock population of flies was split into 3 selected and 3 unselected lines. Selected lines were infected each generation with P. aeruginosa PA01 and the survivors allowed to populate the subsequent generation. Selection was applied for 10 generations. For each round of selection, 8 sub-replicate populations consisting of 20 flies per sex were infected with P. aeruginosa. Mated flies aged to 4-7 days old were anaesthetized with CO2 and infected as previously described by dipping a sterile needle in the bacterial culture and piercing the intrathoracic region of the fly . Fly mortality was then monitored for each population over 48 hours. Survivors from each of the 8 sub-replicates were pooled into a single population to seed the subsequent generation. The unselected lines were not infected during the breeding regime, but were exposed to the same bottleneck in population size as their paired selected lines by randomly selecting a set of individuals to found the next generation. Survival in response to infection was monitored each generation in selected lines and at two time points, G6 and G10 in the unselected lines. While survival rates in the selected lines rose rapidly from 15% at G0 to 70% by G10, survival rates of unselected lines remained low. Selection was then removed at G10 for 5 generations and survival of the selected lines reverted to pre selection levels .
Immunity Gene Expression at G10 after selection
There were three replicate lines in the original selection experiment (S1-S3) that all exhibited similar levels of antibacterial defense after selection. For this study we have focused on just one of the lines (S1) and its paired unselected line (C1). At G10, 4 to 7 day old whole female flies from these two lines were either collected directly for constitutive measures of expression or first infected with P. aegurinosa to measure infection-induced expression. For the infected treatment, flies were collected 8 hours post infection with P. aeruginosa. Five pools of 20 flies each were collected to represent each treatment for all combinations of selected/unselected × infected/uninfected. All flies were snap frozen in liquid nitrogen and extracted for Total RNA using Trizol (Invitrogen Corp., Carlsbad, CA) as per the manufacturer's instructions. Total RNA was treated with 2 μl of DNase I (Roche) for 30 minutes at 37°C in a 20 μl reaction to eliminate genomic DNA.
Real-time PCR (RT-PCR) was used to measure the expression of a set of 8 genes identified from our previous microarray study as they were shown to be involved in improved defense to P. aeruginosa infection in male flies . Four of the genes represent the humoral (DiptB, dros4, PGRP-SB1, PGRP-SD) and 4 represent the cellular (Bc, eater, Sr-CI, TepII) immune response. Primer sequences of 6 of the genes were as reported in the initial study. In addition we designed primers for PGRP-SB1 (FBtr0075348); Forward 5'-3' TTAGCTCTATCCGCCAATGC, Reverse 5'-3' CCCTTGTGATCCGACTGAAT that generated a 191 bp product and DiptB (FBtr0086621); Forward 5'-3' CTGGCATATGCTCCCAATTT, Reverse 5'-3' ATAGGGTCCACCAAGGTGCT that generated a 198 bp product. The gene fus, (FBgn0023441), which is involved in an epidermal growth factor receptor signaling pathway , was employed a reference gene (Forward 5'-3' AAAGTGGTGGAAGCAACAGG, Reverse 5'-3' CGCACACAAACTCGAAAAGA, 158 bp product) as it is was shown not to change in response to infection or selection in the original microarray analysis . There is also no evidence that female expression should differ from that of males per the meta analysis of 13 studies (Male/Female ratio = 1.07, P = 0.21) in the Sebida database .
Approximately 0.5 μg of total RNA was first reverse transcribed using random primers and SuperScript III reverse transcriptase (Invitrogen) according to manufacturer's protocols. Then RT-PCR was performed on Rotor-gene 6000 (Corbett Life Science, Sydney, NSW) using Platinum®SYBR®Green (Invitrogen Inc, Carlsbad, CA) according to manufacturer's instructions. For each sample, a mastermix of 2 μl RNase-free water, 5 μl of SYBR Supermix and 0.5 μl of each primer (10 μM) was added to 2 μl of cDNA. Three technical replicates were run for each biological replicate. The cycling protocol was as follows; 1 cycle UDG incubation at 50°C for 2 minutes, 1 cycle Taq activation at 95°C for 2 minutes, 40 cycles of denaturation at 95°C for 5 s, annealing at 60°C for 5 s, extension at 72°C for 15 s, fluorescence acquisition at 78°C, and 1 cycle of melt curve analysis from 68-95°C in 1°C steps. Mean Cycle Threshold (CT) and mean amplification efficiency (E) per biological replicate was calculated across technical replicates using Rotor-Gene 6000 Series Software ver.1.7.75 (Info-ZIP Pty Ltd.) for analysis.
Response to Selection
Gene specific general linear models for the effects of selection and infection on expression of humoral genes
Gene specific general linear models for the effects of selection and infection on expression of cellular genes
Comparison to Previous Study
The genes studied here were identified from a previous microarray study in males where their transcription was shown to slightly increase (~2 fold on average) in response to selection in the presence of infection . The expression of several of these genes (PGRP-SD, Sr-CI, TepII, eater, Bc) was also confirmed in males using qRT-PCR. In general selected females demonstrate increases in transcription of the same immunity genes (Figure 1) with the exception of TepII and Bc. The magnitude of the female transcriptional response to selection also appears greater, but direct comparisons of expression values between males and females are likely to be confounded by their variation in size and cell count (in particular with respect to the embryo) and any potential sex bias in the reference gene.
In contrast to expectation, we demonstrated that selection either produced similar investments into constitutive and infection-induced transcription (dros4) or invested preferentially in the inducible response (PGRP-SB1, DiptB, Sr-CI, eater, Bc). For PGRP-SB1, DiptB and dros4 heightened inducible expression should confer improved humoral immunity . A failure to see selection induced increases in constitutive expression of these genes may be explained either by pleiotropic or energetic costs or a decoupling of expression between larval and adult stages. Recent work examining the constitutive expression of the immunity genes, Dipt and dro, has shown that they vary with respect to how they are controlled, demonstrating decoupled and coupled patterns of expression, respectively between larvae and adults . If the genes studied here have decoupled control, any changes in the constitutive expression in larvae, may not be present in the adult where we surveyed expression. In contrast, dros4 may be a candidate for coupled control where constitutive changes in larvae are maintained into adulthood.
For the cellular genes, the benefits of increased transcription upon exposure to infection are not clear. Previous genome wide screens of transcription have not revealed induced expression of Sr-CI or eater upon immune challenge with bacteria [4, 5, 23, 24]. Both genes encode transmembrane proteins in insect plasmatocytes that act as scavenger receptors and are required for recognition of pathogens and effective phagocytosis [25, 26]. As they are not secreted, it is difficult to understand how their enhanced transcription upon infection would produce immediate benefits with respect to microbial defense. The increased transcription seen here in response to selection may instead serve to directly trigger transcription of other more effective genes. Although the molecular cross talk between the various immunity pathways is not well characterised , it has been demonstrated for example that expression of eater is required for complete activation of Defensin . Alternatively, the heightened expression of these genes may not be adaptive. The selected lines may have arrived at genetic solutions where transcription of immunity genes is heightened across the board in a non-specific manner via transcription factors. As above the failure to see increases in constitutive expression in the cellular genes suggests they may be strong candidates for a decoupled model of larval to adult genetic control. In any case, these findings suggest that there is not a global increase in the constitutive expression of cellular immunity genes into adulthood that may explain the rapid development via pleiotropy or the costly tradeoffs seen in host fitness like longevity and egg viability. These costs may instead result from larval phase transcription investments.
The transcription of only one gene, TepII, clearly did not respond to selection. TepII's capacity to respond to bacterial challenge [4, 5, 13, 24] and in particular infection with a particular P. aeruginosa strain  has been documented previously. There is also evidence that the transcriptional regulation of TepII is at least in partly under the control of the Toll & Imd pathways . More akin to the humoral genes, the encoded protein, is secreted and acts as an opsonin during the process of targeting foreign bacterial invaders for destruction . Production of TepII is required for phagocytosis of bacteria  and is involved with wound healing . A failure to improve transcription levels of TepII upon selection may indicate a lack of genetic diversity in the population, although this is unlikely given the response of all other genes. There may also be negative pleiotropic effects associated with heightened transcription. Evidence of TepII's involvement in cellular functions other than immunity is lacking, although it expressed in oocytes . Interestingly, a number of studies have demonstrated strong signatures of positive selection in the protein coding regions of TepII [33–35] where it is thought that variable sequence diversity and the production of multiple isoforms  may underpin the ability to bind diverse pathogens. The evolution of such diversity within the gene rather may provide a more effective avenue to improved defense for the insect rather than increased transcription. This is in contrast with a recent study on eater, suggesting that selection for intronic variation rather than protein-coding diversity may have formed the basis of a selective sweep in Drosophila .
The fold increases in transcription brought about by selection ranged only from 2 to 3 for all of the various genes. Given that the genes vary widely in their total level of transcription, it is unlikely that the limiters to the degree of transcriptional improvement can be ascribed to cost. Instead it suggests that a shared evolutionary mechanism for discrete increases in transcriptional output may be operating for multiple genes or at least independently at least once for each immune pathway represented here. One possible model includes an evolved increase in the multiplicity of cis-regulatory regions [38, 39] as that has been documented previously for yeast , Drosophila  and C. elegans .
The authors wish to thank Simon Blomberg for advice on data analysis. This research was supported by a UQ development grant to EAM.
- Ye YH, Chenoweth SF, McGraw EA: Effective but costly, evolved mechanisms of defense against a virulent opportunistic pathogen in Drosophila melanogaster. PLoS Pathog. 2009, 5 (4): e1000385-10.1371/journal.ppat.1000385.PubMedPubMed CentralView ArticleGoogle Scholar
- Zuk M, Stoehr AM: Immune defense and host life history. Am Nat. 2002, 160: S9-S22. 10.1086/342131.PubMedView ArticleGoogle Scholar
- Lemaitre B, Hoffmann J: The host defense of Drosophila melanogaster. Annu Rev Immunol. 2007, 25: 697-743. 10.1146/annurev.immunol.25.022106.141615.PubMedView ArticleGoogle Scholar
- De Gregorio E, Spellman PT, Rubin GM, Lemaitre B: Genome-wide analysis of the Drosophila immune response by using oligonucleotide microarrays. Proc Natl Acad Sci USA. 2001, 98 (22): 12590-12595. 10.1073/pnas.221458698.PubMedPubMed CentralView ArticleGoogle Scholar
- Sackton T: Comparative profiling of the transcriptional response to infection in two species of Drosophila by sort-read cDNA sequencing. BMC Genomics. 2009, 10: 259-10.1186/1471-2164-10-259.PubMedPubMed CentralView ArticleGoogle Scholar
- Fedorka KM, Zuk M, Mousseau TA: Immune suppression and the cost of reproduction in the ground cricket, Allonemobius socius. Evolution. 2004, 58: 2478-2485.PubMedView ArticleGoogle Scholar
- Braun A, Hoffmann JA, Meister M: Analysis of the Drosophila host defense in domino mutant larvae, which are devoid of hemocytes. Proc Natl Acad Sci USA. 1998, 95 (24): 14337-14342. 10.1073/pnas.95.24.14337.PubMedPubMed CentralView ArticleGoogle Scholar
- Wasserman SA: Toll signalling: the enigma variations. Curr Op Genet Dev. 2000, 10: 497-502. 10.1016/S0959-437X(00)00118-0.PubMedView ArticleGoogle Scholar
- Wood W, Jacinto A: Drosophila melanogaster embryonic haemocytes: masters of multitasking. Nat Rev Mol Cell Biol. 2007, 8 (7): 542-551. 10.1038/nrm2202.PubMedView ArticleGoogle Scholar
- Ahtiainen JJ, Alatalo RV, Kortet R, Rantala MJ: A trade-off between sexual signalling and immune function in a natural population of the drumming wolf spider Hygrolycosa rubrofasciata. J Evol Biol. 2005, 18: 985-991. 10.1111/j.1420-9101.2005.00907.x.PubMedView ArticleGoogle Scholar
- Simmons LW, Roberts B: Bacterial immunity traded for sperm viability in male crickets. Science. 2005, 309 (5743): 2031-10.1126/science.1114500.PubMedView ArticleGoogle Scholar
- Kraaijeveld AR, Godfray HC: Trade-off between parasitoid resistance and larval competitive ability in Drosophila melanogaster. Nature. 1997, 389 (6648): 278-280. 10.1038/38483.PubMedView ArticleGoogle Scholar
- Meister M, Lagueux M: Drosophila blood cells. Cell Microbiol. 2003, 5 (9): 573-580. 10.1046/j.1462-5822.2003.00302.x.PubMedView ArticleGoogle Scholar
- Anderson AR, Hoffmann AA, McKechnie SW, Umina PA, Weeks AR: The latitudinal cline in the In(3R)Payne inversion polymorphism has shifted in the last 20 years in Australian Drosophila melanogaster populations. Mol Ecol. 2005, 14 (3): 851-858. 10.1111/j.1365-294X.2005.02445.x.PubMedView ArticleGoogle Scholar
- Huston WM, Jennings MP, McEwan AG: The multicopper oxidase of Pseudomonas aeruginosa is a ferroxidase with a central role in iron acquisition. Mol Microbiol. 2002, 45 (6): 1741-1750. 10.1046/j.1365-2958.2002.03132.x.PubMedView ArticleGoogle Scholar
- D'Argenio DA, Gallagher LA, Berg CA, Manoil C: Drosophila as a model host for Pseudomonas aeruginosa infection. J Bacteriol. 2001, 183 (4): 1466-1471. 10.1128/JB.183.4.1466-1471.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Vodovar N, Acosta C, Lemaitre B, Boccard F: Drosophila: a polyvalent model to decipher host-pathogen interactions. Trends Microbiol. 2004, 12 (5): 235-242. 10.1016/j.tim.2004.03.007.PubMedView ArticleGoogle Scholar
- Wakabayashi-Ito N, Belvin MP, Bluestein DA, Anderson KV: fusilli, an essential gene with a maternal role in Drosophila embryonic dorsal-ventral patterning. Dev Biol. 2001, 229: 44-54. 10.1006/dbio.2000.9954.PubMedView ArticleGoogle Scholar
- Gnad F, Parsch J: Sebida: a database for the functional and evolutionary analysis of genes with sex-biased expression. Bioinformatics. 2006, 22: 2577-2579. 10.1093/bioinformatics/btl422.PubMedView ArticleGoogle Scholar
- Simon P: Q-Gene: processing quantitative real-time RT-PCR data. Bioinformatics. 2003, 19 (11): 1439-1440. 10.1093/bioinformatics/btg157.PubMedView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J Roy Statist Soc Ser B. 1995, 57 (1): 289-300.Google Scholar
- Fellous S, Lazzaro BP: Potential for evolutionary coupling and decoupling of larval and adult immune gene expression. Mol Ecol. 2011, 20: 1558-1567. 10.1111/j.1365-294X.2011.05006.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Apidianakis Y, Mindrinos MN, Xiao W, Lau GW, Baldini RL, Davis RW, Rahme LG: Profiling early infection responses:Pseudomonas aeruginosa eludes host defenses by suppressing antimicrobial peptide gene expression. Proc Natl Acad Sci USA. 2005, 102 (7): 2573-2578. 10.1073/pnas.0409588102.PubMedPubMed CentralView ArticleGoogle Scholar
- Irving P, Troxler L, Heuer TS, Belvin M, Kopczynski C, Reichhart JM, Hoffmann JA, Hetru C: A genome-wide analysis of immune responses in Drosophila. Proc Natl Acad Sci USA. 2001, 98 (26): 15119-15124. 10.1073/pnas.261573998.PubMedPubMed CentralView ArticleGoogle Scholar
- Kocks C, Cho JH, Nehme N, Ulvila J, Pearson AM, Meister M, Strom C, Conto SL, Hetru C, Stuart LM, et al: Eater, a transmembrane protein mediating phagocytosis of bacterial pathogens in Drosophila. Cell. 2005, 123 (2): 335-346. 10.1016/j.cell.2005.08.034.PubMedView ArticleGoogle Scholar
- Ramet M, Pearson A, Manfruelli P, Li X, Koziel H, Gobel V, Chung E, Krieger M, Ezekowitz RA: Drosophila scavenger receptor CI is a pattern recognition receptor for bacteria. Immunity. 2001, 15 (6): 1027-1038. 10.1016/S1074-7613(01)00249-7.PubMedView ArticleGoogle Scholar
- Erturk-Hasdemir D, Silverman N: Eater: a big bite into phagocytosis. Cell. 2005, 123: 190-192. 10.1016/j.cell.2005.10.005.PubMedView ArticleGoogle Scholar
- Brennan CA, Delaney JR, Schneider DS, Anderson KV: Psidin is required in Drosophila blood cells for both phagocytic degradation and immune activation in the fat body. Curr Biol. 2007, 17: 67-72. 10.1016/j.cub.2006.11.026.PubMedView ArticleGoogle Scholar
- Blandin S, Levashina EA: Thioester-containing proteins and insect immunity. Mol Immunol. 2004, 40: 903-908. 10.1016/j.molimm.2003.10.010.PubMedView ArticleGoogle Scholar
- Stroschein-Stevenson SL, Foley E, O'Farrell PH, Johnson AD: Identification of Drosophila gene products required for phagocytosis of Candida albicans. PLoS Biol. 2006, 4 (1): e4-10.1371/journal.pbio.0040004.PubMedPubMed CentralView ArticleGoogle Scholar
- Vierstraete E, Verleyen P, Baggerman G, D'Hertog W, Van den Bergh G, Arckens L, De Loof A, Schoofs L: A proteomic approach for the analysis of instantly released wound and immune proteins in Drosophila melanogaster hemolymph. Proc Natl Acad Sci USA. 2004, 101: 470-475. 10.1073/pnas.0304567101.PubMedPubMed CentralView ArticleGoogle Scholar
- Nakahara K, Kim K, Sciulli C, Dowd SR, Minden JS, Carthew RW: Targets of microRNA regulation in the Drosophila oocyte proteome. Proc Natl Acad Sci USA. 2005, 102: 12023-12028. 10.1073/pnas.0500053102.PubMedPubMed CentralView ArticleGoogle Scholar
- Jiggins FM, Kim KW: A screen for immunity genes evolving under positive selection in Drosophila. J Evol Biol. 2007, 20: 965-970. 10.1111/j.1420-9101.2007.01305.x.PubMedView ArticleGoogle Scholar
- Sackton TB, Lazzaro BP, Schlenke TA, Evans JD, Hultmark D, Clark AG: Dynamic evolution of the innate immune system in Drosophila. Nat Genet. 2007, 39: 1461-1468. 10.1038/ng.2007.60.PubMedView ArticleGoogle Scholar
- Heger A, Ponting CP: Evolutionary rate analyses of orthologs and paralogs from 12 Drosophila genomes. Genome Res. 2007, 17: 1837-1849. 10.1101/gr.6249707.PubMedPubMed CentralView ArticleGoogle Scholar
- Lagueux M, Perrodou E, Levashina EA, Capovilla M, Hoffmann JA: Constitutive expression of a complement-like protein in toll and JAK gain-of-function mutants of Drosophila. Proc Natl Acad Sci USA. 2000, 97 (21): 11427-11432.PubMedPubMed CentralView ArticleGoogle Scholar
- Juneja P, Lazzaro BP: Haplotype structure and expression divergence at the Drosophila cellular immune gene eater. Mol Biol Evol. 2010, 27: 2284-2299. 10.1093/molbev/msq114.PubMedPubMed CentralView ArticleGoogle Scholar
- Gerland U, Hwa T: On the selection and evolution of regulatory DNA motifs. J Mol Evol. 2002, 55: 386-400. 10.1007/s00239-002-2335-z.PubMedView ArticleGoogle Scholar
- Wray GA, Hahn MW, Amores A, Balhoff JP, Pizer M, Rockman MV, Romano LA: The evolution of transcriptional regulation in eukaryotes. Mol Biol Evol. 2003, 20: 1377-1419. 10.1093/molbev/msg140.PubMedView ArticleGoogle Scholar
- Paixao T, Azevedo RB: Redundancy and the evolution of cis-regulatory element multiplicity. PLoS Comput Biol. 2010, 6: e1000848-10.1371/journal.pcbi.1000848.PubMedPubMed CentralView ArticleGoogle Scholar
- Ludwig MZ, Patel NH, Kreitman M: Functional analysis of eve stripe 2 enhancer evolution in Drosophila: rules governing conservation and change. Devel. 1998, 125: 949-958.Google Scholar
- O'Meara MM, Bigelow H, Flibotte S, Etchberger JF, Moerman DG: Cis-regulatory mutations in the Caenorhabditis elegans homeobox gene locus cog-1 affect neuronal development. Genetics. 2009, 181: 1679-1686. 10.1534/genetics.108.097832.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.