Modularity and evolutionary constraints in a baculovirus gene regulatory network
© Oliveira et al.; licensee BioMed Central Ltd. 2013
Received: 18 February 2013
Accepted: 15 August 2013
Published: 4 September 2013
The structure of regulatory networks remains an open question in our understanding of complex biological systems. Interactions during complete viral life cycles present unique opportunities to understand how host-parasite network take shape and behave. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) is a large double-stranded DNA virus, whose genome may encode for 152 open reading frames (ORFs). Here we present the analysis of the ordered cascade of the AgMNPV gene expression.
We observed an earlier onset of the expression than previously reported for other baculoviruses, especially for genes involved in DNA replication. Most ORFs were expressed at higher levels in a more permissive host cell line. Genes with more than one copy in the genome had distinct expression profiles, which could indicate the acquisition of new functionalities. The transcription gene regulatory network (GRN) for 149 ORFs had a modular topology comprising five communities of highly interconnected nodes that separated key genes that are functionally related on different communities, possibly maximizing redundancy and GRN robustness by compartmentalization of important functions. Core conserved functions showed expression synchronicity, distinct GRN features and significantly less genetic diversity, consistent with evolutionary constraints imposed in key elements of biological systems. This reduced genetic diversity also had a positive correlation with the importance of the gene in our estimated GRN, supporting a relationship between phylogenetic data of baculovirus genes and network features inferred from expression data. We also observed that gene arrangement in overlapping transcripts was conserved among related baculoviruses, suggesting a principle of genome organization.
Albeit with a reduced number of nodes (149), the AgMNPV GRN had a topology and key characteristics similar to those observed in complex cellular organisms, which indicates that modularity may be a general feature of biological gene regulatory networks.
Cellular pathways and gene regulatory networks (GRNs) are complex systems that emerge under natural selection and bare distinct evolutionary proprieties such as modularity. Modularity is an effective mechanism for keeping perturbations confined while preserving the complete system . It is widely observed in various organisms and is possibly a fundamental biological design principle [2–4]. Emergent proprieties in networks also impose constraints on individual genes, mostly due to epistasis, leading to gene conservation. Important enzymes and other essential proteins – as reported by flux balance analysis [5, 6] – tend to vary less than those under lower functional load , indicating that the flow of matter through metabolic networks expounds an evolutionary constraint imposed on components of any given pathway. These findings where confirmed for yeast, in which highly connected enzymes evolve more slowly than less connected ones . Viruses intertwine their metabolic functions with those of the host cell, which allows infected cell to be understood as superorganisms, comprising the minimal essential condition for viral replication and biomagnification ; yet how viruses organize their gene network remains an open question.
Baculoviruses provide a rich environment where to investigate these relations, since these viruses have many completely sequenced genomes, each one encoding around one hundred plus open reading frames (ORFs) with transcription strategies and defined overall temporal regulation. The Baculoviridae constitute a family of invertebrate viruses that infect mainly insects of the order Lepidoptera, with large circular, covalently closed, double-stranded DNA genome [10, 11]. The analysis of fifty-seven complete baculovirus genomes has shown 37 genes probably shared among all of them [12–17], most of which are involved in essential processes, such as replication, transcription and oral infectivity. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) has been used in Brazil and other countries to control the velvet bean caterpillar Anticarsia gemmatalis (Lepidoptera: Noctuidae), an important pest of soybean crops [18, 19]. The prototype AgMNPV isolate 2D (AgMNPV-2D) [20, 21] has a genome of 132,239 bp, which may encode 152 ORFs .
Gene expression and replication cycle of baculovirus appear to be regulated by viral factors and the cellular milieu, where distinct gene classes are thought to be trans-activated (directly or indirectly) by closely related transcription complexes, causing expression to be controlled by distinct promoters, activated in a temporal concerted fashion during the infection cycle . Three main transcriptional phases can be distinguished: early (immediate early and delayed early), late and very late. The early phase precedes the onset of viral DNA replication and includes transcription of genes involved in host modulation, viral DNA replication and regulation of delayed early and late gene expression. The late and very late phases follow DNA replication and include the expression of genes required for virus assembly and occlusion (structural genes) .
An interesting feature of baculovirus transcription is that some ORFs can be transcribed in tandem, spanning up to seven units [24–26]. Specific genomic regions encode a wide variety of transcripts with different lengths, the so-called ‘overlapping transcripts’ [27–29]. Several tandem transcripts were mapped to different genomic regions of Autographa californica nucleopolyhedrovirus (AcMNPV) [27, 30] and Bombyx mori nucleopolyhedrovirus (BmNPV) genomes . The expression profiles of six AgMNPV genes were investigated: egt; p10; v-trex; helicase; iap-3 and p74. Nevertheless, there was no comprehensive study done so far on its complete gene content transcription and regulation during the replication cycle. Therefore, it is of interest to address the role of gene expression regulation in shaping baculovirus genome organization and evolution. Here we analyzed transcriptional organization of AgMNPV compared to randomized datasets in two distinct cell lines with different infection kinetics . Our analysis focused on: (i) the viral transcriptome profile, (ii) the structural properties of its gene regulatory network (GRN), (iii) genomic arrangement of transcripts of distinct viral cycle phases and, (iv) association of transcripts to distinct temporal promoter types.
Viral genes are expressed sooner and at higher amounts in its permissive cell line
The inferred gene regulatory network has non-random proprieties
Since we had three independent replicates of infection experiments lasting 48 hours each, we used a single 144 hrs time series to infer a putative AgMNPV-2D gene regulatory network (GRN) from 149 ORFs unfolded during infection in UFL-AG-286 and likewise, a single 144 hrs time series for IPLB-SF-9. Under the assumption of a recurrent temporal regulation program for viral gene expression, this approach should maximize true gene-to-gene cross-correlation signal present in the data, while mimicking three successive infection cycles in a row. Ranked lists of directed node interactions were generated with GENIE3 from expression data of AgMNPV-2D infecting UFL-AG-286, infecting IPLB-SF-9 and the respective randomized data (UFL-AG-286rand and IPLB-SF-9rand). One thousand links were necessary on average to assemble a single connected component with 149 nodes for all GRN (including those generated from random datasets). Nevertheless, the complete network for AgMNPV-2D on UFL-AG-286 needed 1120 edges to include all nodes into a single component, which was in agreement with similar but random GRNs, which had an average of 1030 edges. Preliminary exercises with random GRNs revealed that the data generated from our complete lists of directed node interactions generated with GENIE3 lacked essential complex features of known biological networks such as modularity. Since we needed more realistic null models to better discern features of real networks from the data, we decided to keep only the relevant links that would be sufficient to generate a single component network while generating randomized networks by shuffling the real data time points. This procedure obliterated the time structure of the real data but kept its value range and sampling density.
The network architecture displays increasing modularity
It is of utter importance to emphasize that node linkages in our GRNs depict the relationships among regulatory elements inferred from gene expression profile (transcriptional data). Hence, it has to be made clear that the networks we inferred convey the implied associations among transcription regulation elements per se, and not necessarily gene product function relationships (such as in an interactome). Fittingly, genetic regulatory functions are best understood when viewed in terms of intergenic linkages of diverse modalities, and these non-linear functions are not visible at the level of any individual gene .
The gene communities (modules) have redundant functions
Gene function impacts on expression level variation and on GRN topology
Given the apparent assortativity (that is, the preferential attachment of nodes by shared properties) of both gene function and promoter types, we then compared expression profiles among genes with correlated functions, such as: (i) replication: dnapol, lef-1, lef-2 and helicase; (ii) transcription: lef-4, lef-5, lef-8, lef-9, p47 and vlf-1; (iii) structural proteins: pif-1, pif-2, gp41, odv-e56, odv-ec27, p74, p6.9, vp91, vp39, vp1054 and p33; (iv) core genes with unknown function; (v) some immediate-early genes that do not belong to the 37 conserved genes (cg30, ie-0, ie-1, ie-2, me53 and pe38); and carried out the statistical analysis (ANOVA) among groups followed by Tukey’s multiple comparison test, assuming a significance level at p < 0.05. We noticed that early genes and those associated with DNA replication had less difference in their temporal and quantitative profiles between the two cell lines (Figure 1). This result suggested that core conserved functions do appear to have more invariance in the temporal structure of their expression program. Crucially, this temporal consistency was also associated with a trend on the association of a key GRN topologic feature and gene variability, since core functions (replication and capsid) had a higher median betweenness centrality (BC = 0.0019, Additional file 4: Figure S2) with a large variance and a significantly lower genetic diversity (θ = 0.42 ± 0.035), while satellite (all the others) functions had lower median betweenness centrality (BC = 0.0013) with large variance and a significantly higher genetic diversity (θ = 0.47 ± 0.044, Additional file 4: Figure S2).
The spatial organization of regions of overlapping transcription (ROTs) is important and conserved
To better infer the relation between time expression and gene position, we built a syntenic map for the complete genome of thirteen Group I alphabaculoviruses (Additional file 5: Figure S3 and Additional file 6: Table S3) with twenty-two highly conserved ‘regions of overlapping transcription’ (ROTs) defined as multiply transcribed genomic loci encoding tandem RNA transcripts of different lengths. These ROTs had almost the same overall order and relative orientation in all genomes, even when genomic reorganization events, such as inversions and gene deletions took place. When compared to other genomic regions, for all 13 genomes investigated, we found in ROTs a higher frequency of ORFs with very short or absent intergenic spacers (overlapping genes) (χ2 = 65.041, p < 0.00001; d.f. = 1). If ROTs entail the synthesis of overlapping transcripts in alphabaculoviruses, it is possible that approximately 40% of the AgMNPV-2D genes (61 of 152) could be transcribed in tandem.
Expression of structural and regulatory genes vary more on less permissive cells
The use of two different infection kinetics allowed us to observe contrasting aspects of the transcriptional program of AgMNPV-2D in vitro, providing comprehensive, time-related information on the transcription dynamic and quantitative information on the expression of 149 of its 152 predicted ORFs. The ‘sooner and higher’ expression pattern in UFL-AG-286 compared to IPLB-SF-9 could be explained by the different states of adaptation of the virus to distinct cellular environments, since the infection in IPLB-SF-9 (a less permissive cell line) could be understood as a ‘perturbation’ of the AgMNPV-2D GRN. Nevertheless, we have not noticed significant difference in most genes related with DNA replication, such as dnapol (ORF 65), lef-1 (ORF 19) and lef-2 (ORF 3) that were expressed at the same time and at similar levels in both cell lines. While replication-associated genes had similar expression, most structural genes had a significant delay in IPLB-SF-9 and notable differences in the magnitude of expression (Additional file 1: Table S1 and Additional file 6: Table S3). Hence, the expression profiles of the structural genes found in IPLB-SF-9 could help explain the delay in viral morphogenesis routinely observed in our laboratory that was also reported for IPLB-SF-21 cells .
We consistently observed earlier expression for most of the genes than expected based on the observed for other baculoviruses. For example, the p6.9 gene (ORF 96) was reported as late, expressed at 12 h p.i. in the AcMNPV . P6.9 is a small protein (6.9 kDa), basic, arginine-rich and it is involved in the condensation of viral DNA. It is not surprising to find this gene expressed earlier in the AgMNPV-2D viral cycle, since DNA replication begins at 2 h p.i.  and P6.9 is known to be necessary to package and stabilize the newly synthesized viral genome . In addition, Iwanaga and colleagues  also found p6.9 mRNA and other transcripts earlier in infection, which corroborates our results and earlier observations on other DNA viruses, such as FV-3 and CIV [47–49]. Crucially, the difference we observed could be due to the use of real-time PCR, which is much more efficient and requires fewer targets than hybridization techniques (such as Northern blot), detecting low amounts of RNA initially present in each sampling time. This precocity could simply reflect differences between baculoviruses, given the earlier onset of DNA replication in AgMNPV-2D (2 h p.i.) compared to others, such as AcMNPV (5 h p.i.) , SfMNPV (8 h p.i.)  and LdMNPV (18 h p.i.) .
Furthermore, we also noticed that genes with more than one copy in the genome (such as iaps, pifs and ptps) had distinct expression profiles. For example: (i) iap-3 (ORF 34) was expressed earlier in UFL-AG-286 cells; (ii) iap-2 (ORF 70) simultaneously, whereas; (iii) iap-1 (ORF 40) showed expression only later, in both cell lines, being more expressed than iap-2 and iap-3. Gene duplication and redundancy are important mechanisms that allow adaptive evolution of genomes, and represent 8-20% of the gene content in eukaryotes . We would argue that this allows one of the copies of a gene to evolve with lesser functional constraint , which eventually could help it acquiring a new functionality (such as, iap and bro genes).
Conserved genes have higher centrality in a modular AgMNPV-2D GRN
By comparing with our randomized models the inferred GRN in UFL-AG-286 cells, where five communities of nodes were identified, it is notable that it had a significantly different modular structure. This was made clear by the trend in both the average connectivity and the average coefficient distribution values for nodes as function of the number of neighbors they had (Additional file 3: Figure S1). By looking at the AgMNPV-2D GRN as a single component we observed that nodes with highest importance, both in terms of BC and SC, included relevant host-interaction functions such as gp64 (ORF 124) and egt (ORF 18). The gp64 codes for an envelope glycoprotein paramount for effective membrane fusion during the secondary infection that is highly conserved in the group I nucleopolyhedroviruses (NPVs) . Interestingly, the expression of the egt gene in Lymantria dispar multiple nucleopolyhedrovirus (LdMNPV) was shown to be the genetic basis of climbing behavior and tree top disease in silkworm infected with baculovirus, which facilitates preferential predation of sick larvae while maximizing viral spread in the environment .
A relevant property of our inferential procedure was that there was a positive correlation of reduced genetic diversity (θ) and importance of the gene in the GRN, as expected for key components of metabolic networks [7, 8]. Core conserved functions also had more temporal consistency of expression and distinct topologic features in the GRN. These coherent results support our modeling effort, since they constitute independent evidence for a relationship among phylogenetic data of baculovirus genes and complex GRN attributes inferred from expression data alone.
The observed modular heterogeneity that separates key genes that are functionally related on different communities could help accommodate redundancy and shield important functions, therefore maximizing GRN robustness by compartmentalization . For example, genes possibly associated with anti-apoptotic function (iap genes)  were well distributed among three communities (I, II and IV), while dnapol and helicase nested in different communities (III and V). Likewise, assigning several different functions to a given module would also provide robustness by increasing functional redundancy . Nevertheless, some functional associations were evident. Community I was highly populated by structural proteins associated with capsid and virion formation and promoters with early & late motifs (Figure 4). Interestingly, in the AgMNPV-2D 48 of 56 nodes in community I (13 of which are core genes) are conserved in almost all group I alphabaculoviruses sequenced so far. This feature was also observed in community II, in which one third of the nodes consisted of highly conserved genes, many of which were genes associated with DNA replication. The clustering of genes in community II with unidentified regulatory motifs could entail promoter motifs not yet described , which is made more compelling giving the high connectivity among its members. Spatially or chemically isolated modules composed of several cellular components and carrying discrete functions are considered fundamental building blocks of cellular organization, but their presence in highly integrated biochemical networks lacks quantitative support. In a key comparative study, the metabolic networks of 43 distinct organisms were shown to have modular organization and in E. coli, the hierarchical modularity had a good correlation with metabolic functionality . Spatially or chemically isolated modules composed of several cellular components and carrying discrete functions are considered fundamental building blocks of cellular organization, but their incidence and generality needs to be better investigated across a wider spectrum of organismic complexity, such as in viruses.
Unfortunately, at this time we are not able to wire cell-coded nodes in the AgMNPV-2D GRN. Nevertheless, we would argue that they should have not only a crucial role in the viral cycle, but also should be connected to the real viral GRN . However, viral GRNs encompass both the cellular and viral encoded functions necessary to complete the viral multiplication cycle . Accordingly, preliminary data on RNA suppression subtractive hybridization (SSH) experiments (Oliveira et al., in preparation) indicated that at least 50 cellular distinct genes were down regulated during infection at 24 h.p.i. (Additional file 11: Table S7). By using both Gene Ontology (http://www.geneontology.org/) and EGene (http://www.coccidia.icb.usp.br/egene/) programs while using the eggNOG2.0 database (http://eggnog.embl.de/), we found sequences coding mainly for proteins related to metabolism or protein modification (33.9%), ion transport and energy (12.0%), nucleic acids metabolism (11.3%) and several other functions at a lower frequency. But importantly, we found cellular genes at lower frequency, which are of direct relevance to viral host interaction such as: (i) the significantly hypo-expressed anterior fat body protein (AFP), (ii) defensin and, (iii) heat shock 90 (hsp90) (Additional file 11: Table S7). At this juncture, it is hard to assess the extent at which the hypo-expression of cellular genes does not simply reflect cellular degradation due to viral disruption of cellular homeostasis, or is due to interactions with viral gene products. Notwithstanding, these preliminary results could imply that at least 50 additional cellular genes should be wired to the viral GRN for a proper description of the viral life-cycle.
The conserved genome architecture of Group I baculovirus has implications for gene expression
The presence of highly conserved ROTs in closely related baculoviruses, observed in our sinteny analyses, suggests a complex pattern of gene expression regulation, possibly also present in other baculoviruses [28, 31]. The intergenic spacers among ORFs in ROTs are either very short or absent, and overlapping genes are overwhelmingly located in these genomic regions (χ2 = 65.041, p < 0.00001; d.f. = 1). This would provide a competitive advantage, given that keeping genomes short would speed up genome replication at the cost of maximizing interdependence among superimposed genes . Interestingly, we also found evidence that the virus tries to keep similar attributes dispersed around the genome, mainly avoiding their tandem organization (Additional file 9: Table S5, Additional file 10: Table S6 and Additional file 11: Table S7) with a few notable exceptions. A simple explanation for this observation would be that by avoiding clusters of similar attributes, the genome maximizes robustness of the viral GRN. Moreover, whilst ROTs are somehow conserved among group I NPV, our findings suggest that physical proximity is not a reliable predictor of temporal patterns of gene expression, since in the AgMNPV-2D there was no correlation between the position of genes in the genome and their temporal expression (Additional file 7: Figure S4). These results are in line with the possibility that tandem transcripts from ROTs in baculovirus are not necessarily translated at once [24–26]. Moreover, long mRNAs of several genes have been observed late in infection [27, 28, 31, 63].
The region comprising from ORF 129 (p24) to 133 (alk-exo) in AcMNPV encodes for diverse tandem transcripts , which share the same T-rich 3’ end termination site . These ORFs are transcribed as monocistrons at the onset of infection (up until around 2 h p.i.) and later are transcribed in tandem . The synthesis of long transcripts has been related to regulatory functions, such as promoter occlusion  and RNA interference (RNAi) [66, 67]. Promoter occlusion is caused by RNA polymerase transcriptional complex activity during the transcription of upstream genes, blocking the monocistronic expression of downstream genes . The expression of antisense sequences in tandem RNAs has been suggested as a mechanism of down-regulation [68, 69]. Moreover, it has been assumed that overlapping transcripts should have a functional relevance for regulation of gene expression [27–31]. Possibly, monocistronic expression is sufficiently high as to hamper mRNA detection from tandem transcripts, which would help explaining the lack of correlation shown in Additional file 7: Figure S4. Fittingly, it has been suggested that regulatory elements functioning as ‘internal ribosome entry sites’ (IRES) are functional in the IPLB-SF-9 cell line . Furthermore, the capsid protein VP1054 is probably synthesized by IRES-mediated translation , since its monocistronic transcript was not found in AcMNPV and BmNPV gene expression programs .
The transcriptome of the AgMNPV-2D indicated that conserved viral functions showed expression synchronicity, higher betweeness centrality and reduced genetic diversity, consistent with evolutionary constraints present in complex biological systems [7, 8]. Most ORFs were expressed at ‘sooner and higher’ levels in a more permissive cell line. In 13 group I alphabaculovirus genomes, we found 22 highly conserved regions of overlapping transcription, whose extended 3’ UTR could play a role in gene expression regulation . The inferred GRN had a modular topology, comprising five communities of highly interconnected nodes, consisting of genes with different functions, promoter motifs, and physical location in the genome. This modular heterogeneity was suggestive of architectural redundancy that could promote robustness [57, 59]. The fact that a similar architecture was observed on a simpler GRN, such the one we report for the AgMNPV-2D, tends to support the idea that hierarchical modularity may be an extensive and generic self-similarity feature of system-level biological organization .
Primer design and validation
Specific primer pairs for each ORF annotated in the AgMNPV-2D genome  were designed using Oligo Analysis Software v. 6.8 (Molecular Biology Insights) (Additional file 12: Table S8). Each pair was tested by PCR amplification and sequencing. Consensus sequences were verified by BLAST searching to check cross-hybridization against each other and the AgMNPV-2D genome. After this initial validation step, each pair of primers was further tested in real-time PCR reactions to ensure its specificity as determined by a unimodal-melting curve. After the specific amplification was confirmed, primer pairs were also validated by five points dilution (10-2, 10-4, 10-6, 10-8 and 10-10 or 10-3, 10-4, 10-5, 10-6 and 10-7) amplification by real-time PCR. Reaction efficiency of 0.9 and the correlation coefficient (r2) values of 0.99 were assumed as minimal quality standards for data acquisition and melting curves were also analyzed to confirm specificity. ORF 135 (bro-h) was not included in the study because there were unspecific products for three alternative primer pairs designed inside and upstream of its coding region.
Cell lines and infection kinetic: virus infection, RNA extraction and reverse transcription
The UFL-AG-286  and IPLB-SF-9  cell lines were maintained in Grace medium (Gibco) supplemented with 10% fetal bovine serum (Gibco). These permissible cell lines were used since AgMNPV has different infection kinetics in them . UFL-AG-286 was derived from embryonic tissue of the lepidopteron Anticarsia gemmatalis, which is natural host of AgMNPV [19, 21, 73–76], while IPLB-SF-9 is a clone of IPLB-SF-21 cell line, established from pupal ovaries of Spodoptera frugiperda, another Noctuidae species from a different genus . The virus strain AgMNPV-2D  was multiplied only in UFL-AG-286 cells, and both cell lines were used in the kinetic infections. Cultures of 105 UFL-AG-286 e IPLB-SF-9 cells were infected at a multiplicity of infection (MOI) of 10 and incubated at 28°C for 1h to allow infection synchronization. Unattached virions were then removed and fresh culture medium was added to the infected cells. Mock-infected samples were treated in the same manner as virus-infected cells, but with fresh culture medium instead of viral suspension. All experiments were done in triplicate. Total RNA was extracted from infected and mock cells at 0, 1, 3, 5, 7, 9, 11, 13, 24 and 48 hp.i. using RNeasy Plus Kit (Qiagen) according to the manufacturer’s recommendations. RNA samples were subjected to DNase I treatments (DNA-free, Ambion) according to the manufacturer’s protocols. 100 ng of each RNA sample was reverse transcribed using both SuperScript III Reverse (Invitrogen) and High-Capacity cDNA Archive kit (Applied Biosystems) using oligo (dT) and random primers (respectively) to allow maximum conversion of mRNA into cDNA. The cDNAs were mixed and used in the real-time PCR reactions.
We used real-time PCR to quantify mRNAs concentration during infection in cell cultures. This was done because of its accurate qualitative and quantitative measurements of the amount of transcripts with high sensibility and reproducibility [78, 79], outperforming microarray in these regards . Reactions were performed along with the points of the standard curve and the cDNA samples from kinetic infections, to keep the same conditions for standards and experimental samples. Reactions mixes contained 1.0 μl of cDNA from each time point (or point dilution), 7.5 μl of iQ SYBR Green Supermix kit (Bio-Rad) and 0.5 μM of each primer, in a final volume of 15 μl. The cycling conditions were: 2 min at 96°C, followed by 40 cycles of 30 s at 96°C, 30 s at 50°C /52°C /54°C /57°C (depending on the Tm of primer pair), and 40 s at 72°C. Melting curve analysis was performed increasing the temperature of the last cycle (72°C) until reaching 96°C, 1°C per cycle, 5 s at each cycle. Amplification, detection, and data analysis were performed using the Rotor-Gene 3000 systems (Corbett Life Science).
Promoter analysis and synteny in regions of overlapping transcription (ROTs)
Promoters were searched in regions comprising 200 bp upstream of the putative start codon of each ORF by screening for known promoter elements. We chose 200 bp since this was determined by several independent studies to harbor all known late promoters (within the first 80 bp upstream from the first ATG) and most known functional early promoters of baculovirus genes. Early promoter (E) indicates a TATA box sequence (TATA or TATAWWW, W= A or T) followed by a CANT motif downstream, transcribed by the host RNA polymerase II . To screen the late promoter motifs transcribed by the viral encoded RNA polymerase, the conserved pattern used was DTAAG [23, 82]. Moreover, to investigate the genomic organization of clusters of genes with similar transcriptional profiles that we found in the AgMNPV, we did complete genome overlays including 13 group I alphabaculovirus with the Artemis Comparative Tool , while focusing on the conservation of synteny of 22 known regions of overlapping transcription (ROTs) that where first described in AcMNPV  and BmNPV .
A conceptual problem inherent to transcriptional data analyses is the assumption that gene expression dynamics is linear, deterministic and constant, while in fact it is mostly non-linear, infrequent and random [84, 85]. Expression of genes most often happens as pulses of indeterminate duration , and depends on the interplay of different biological components (transcription factors, promoters motifs, RNA polymerases, splicing factors, etc.) at a given temporal (i.e., developmental phase) and spatial (i.e., tissular) context . Moreover, levels of mRNA expression and its encoded protein translation are not necessarily correlated , while remarkable differences are observed in mRNA half-life . Notwithstanding, clustering methods were used to analyze baculovirus gene expression transcriptional data [46, 60, 89]; but because transcript numbers experience iterative exponential fluctuations in time during successive infection cycles, gene expression is better understood as a non-linear system [84, 85], taking place on a log-scaled attractor [90–92]. We argue that these properties make simple least-square distance estimates among individual gene expression time-series uninformative and biased, being inadequate for process description [26, 93]. Therefore, we used a prediction algorithm for the AgMNPV-2D gene regulatory network that is based on multiple regression and tree-based ensemble methods using our gene expression data (Additional files 1 and 2: Tables S1 and S2) . We chose the method implemented in GENIE3 instead of clustering because it: (i) was shown to perform well on reconstructing complex GRNs; (ii) does not make any assumption about the nature of gene regulation; (iii) infers directed GRNs; (iv) can deal with combinatorial and non-linear interactions among genes and; (v) provides a confidence ranking of regulatory links. Ranked lists of regulatory links were calculated with GENIE3 loaded in the free software environment for statistical computing and graphics R v.2.10.0 , then uploaded the ranked list output in Cytoscape 2.8.0 for network visualization and analysis . To evaluate the importance of node connection in the modeled GRN, we used graph-theoretic methods implemented in the Hub Objects Analyzer (Hubba) plug-in for Cytoscape, which were shown for example to recover known yeast interactome data with greater than 70% accuracy . With Hubba, we used Maximum Neighborhood Component (MNC) and Density of Maximum Neighborhood Component (DMNC) algorithms. This approach was helpful to obtain a minimal set of most important edges necessary to fully connect all viral nodes in the GRN.
Measuring complex network attributes
Modular networks have subsets of nodes that are densely connected within and weekly connected between subsets [61, 98]. Therefore, betweenness centrality (BC), which estimates the number of shortest paths from all vertices to all others that pass through a node, is a more informative measure of topologic importance of a node than connectivity (k), which indicates the number of connections of a node irrespective of its neighbor’s connectivity properties. Likewise, we also estimated other relevant topologic attributes that inform on network architecture such as: (i) the average clustering coefficient distribution (ACC) that gives the average of the clustering coefficients for all nodes with all quantities of neighbors, helping to find modularity in networks; (ii) the average neighborhood connectivity (ANC) that is the average connectivity of all neighbors of a node, which also informs on the sub-structuring of a network; (iii) stress centrality (SC) [99, 100] that is the number of shortest paths passing through a node and; (iv) modularity (Q), which measures the number of edges falling within groups minus the expected number in an equivalent network with edges placed at random, quantifying the level of network compartmentalization in communities of nodes. While comparing genomic and GRN organizations we also investigated if three attributes (GRN communities, promoter motifs and functional categories) were spatially organized in the AgMNPV-2D. Hence, we estimated their probability of arbitrary co-occurrence by a Monte Carlo (MC) procedure, in which the attributes were randomized 100 million times without replacement, using a PERL script available from the authors upon request. During our randomization exercises, the co-occurrences were incremented when an attribute randomly sampled from the list was identical to the previously sampled. Thus, the chance to find co-occurrence events of an attribute depends on its frequency in the list. During MC simulations this approach should introduce compositional biases, since it does not reflect the real gene accretion process during baculovirus genome evolution  but critically, it does generate virtual attribute arrangements with the same total number of components as the real genome.
Estimating the evolutionary change on each AgMNPV-2D ORF
We also sought to generate node attributes such as molecular evolution (i.e., genetic diversity), to correlate with the estimated network parameters. To verify this property, we inferred the genetic diversity (θ) in a comparative set of 121 coding sequences shared among 4 group I alphabaculoviruses: AgMNPV (Anticarsia gemmatalis MNPV, DQ813662), CfDEFNPV (Choristoneura fumiferana defective NPV, AY327402), EppoMNPV (Epiphyas postvittana MNPV, AY043265) and CfMNPV (Choristoneura fumiferana MNPV, AF512031). These baculoviruses were chosen due to their close phylogenetic relationship, and the homologues ORFs were selected according to Oliveira et al.. Nucleotide sequences were extracted from GenBank database using the annotation software Artemis . Multiple sequence alignments for each gene were generated using ClustalW  and manually adjusted using Bioedit . Finally, we used a phylogeny-based Markov Chain Monte Carlo (MCMC) Bayesian method implemented in the Coalesce 1.5 beta program  to estimate the genetic diversity (θ = Ne.μ) for each viral gene alignment and compared θ to betweenness centrality (BC) estimates obtained for ‘core’ (replication and capsid) and ‘satellite’ (all the others) functions . With Coalesce, an initial Watterson estimate of θ was used and the MCMC consisted of 10 short chains and 2 long chains of 4000 steps each, during which the viral trees for each gene were optimized and θ estimated.
JVCO, CTB and AI hold FAPESP scholarships (04/12456-0, 09/16740-8 and 12/04818-5), AFB and CCMF hold CAPES-MSc and PhD scholarships and PMAZ holds a CNPq-PQ scholarship. This work was supported financially by FAPESP (Fundação de Amparo a Pesquisa do Estado de São Paulo, process: 2007/55282-0).
- Kitano H: Biological robustness. Nat Rev Genet. 2004, 5: 826-837.PubMedView ArticleGoogle Scholar
- Schlosser G, Wagner GP: Modularity in development and evolution. 2004, London: University of Chicago PressGoogle Scholar
- West-Eberhard MJ: Developmental Plasticity and Evolution. 2003, New York: Oxford University PressGoogle Scholar
- Hartwell LH, Hopfield JJ, Leibler S, Murray AW: From molecular to modular cell biology. Nature. 1999, 402: C47-C52. 10.1038/35011540.PubMedView ArticleGoogle Scholar
- Schilling CH, Edwards JS, Palsson BO: Toward metabolic phenomics: analysis of genomic data using flux balances. Biotechnol Prog. 1999, 15: 288-295. 10.1021/bp9900357.PubMedView ArticleGoogle Scholar
- Varma A, Palsson BO: Metabolic capabilities of Escherichia coli: I. synthesis of biosynthetic precursors and cofactors. J Theor Biol. 1993, 165: 477-502. 10.1006/jtbi.1993.1202.PubMedView ArticleGoogle Scholar
- Wagner A: Does selection mold molecular networks?. Sci STKE. 2003, 2003: PE41-PubMedGoogle Scholar
- Vitkup D, Kharchenko P, Wagner A: Influence of metabolic network structure and function on enzyme evolution. Genome Biol. 2006, 7: R39-10.1186/gb-2006-7-5-r39.PubMedPubMed CentralView ArticleGoogle Scholar
- Krakauer DC, Zanotto PMA: Viral Individuality and Limitations of the Life Concept. Protocells: Bridging Nonliving and Living Matter. Edited by: Rasmussen S, Bedau M, Packard N, Krakauer D, Stadler P, Deamer D. 2009, Cambridge: MIT Press, 712-Google Scholar
- O’Reilly DR, Miller LK, Luckow VA: Baculovirus Expression Vectors. 1993, New York: A Laboratory ManualGoogle Scholar
- Granados RR, Federici BA: Biology Of Baculoviruses. 1987, Boca Raton, Florida, USA: CRC PressGoogle Scholar
- Garavaglia MJ, Miele SA, Iserte JA, Belaich MN, Ghiringhelli PD: The ac53, ac78, ac101, and ac103 genes are newly discovered core genes in the family Baculoviridae. J Virol. 2012, 86: 12069-12079. 10.1128/JVI.01873-12.PubMedPubMed CentralView ArticleGoogle Scholar
- Garcia-Maruniak A, Maruniak JE, Zanotto PM, Doumbouya AE, Liu JC, Merritt TM, Lanoie JS: Sequence analysis of the genome of the Neodiprion sertifer nucleopolyhedrovirus. J Virol. 2004, 78: 7036-7051. 10.1128/JVI.78.13.7036-7051.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Herniou EA, Olszewski JA, Cory JS, O’Reilly DR: The genome sequence and evolution of baculoviruses. Annu Rev Entomol. 2003, 48: 211-234. 10.1146/annurev.ento.48.091801.112756.PubMedView ArticleGoogle Scholar
- Lauzon HA, Lucarotti CJ, Krell PJ, Feng Q, Retnakaran A, Arif BM: Sequence and organization of the Neodiprion lecontei nucleopolyhedrovirus genome. J Virol. 2004, 78: 7023-7035. 10.1128/JVI.78.13.7023-7035.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Miele SA, Garavaglia MJ, Belaich MN, Ghiringhelli PD: Baculovirus: molecular insights on their diversity and conservation. Int J Evol Biol. 2011, 2011: 379424-PubMedPubMed CentralView ArticleGoogle Scholar
- Slack J, Arif BM: The baculoviruses occlusion-derived virus: virion structure and function. Adv Virus Res. 2007, 69: 99-165.PubMedView ArticleGoogle Scholar
- Szewczyk B, Hoyos-Carvajal L, Paluszek M, Skrzecz I, de Lobo Souza M: Baculoviruses– re-emerging biopesticides. Biotechnol Adv. 2006, 24: 143-160. 10.1016/j.biotechadv.2005.09.001.PubMedView ArticleGoogle Scholar
- Moscardi F: Assessment of the application of baculoviruses for control of Lepidoptera. Annu Rev Entomol. 1999, 44: 257-289. 10.1146/annurev.ento.44.1.257.PubMedView ArticleGoogle Scholar
- Maruniak JE, Garcia-Maruniak A, Souza ML, Zanotto PM, Moscardi F: Physical maps and virulence of Anticarsia gemmatalis nucleopolyhedrovirus genomic variants. Arch Virol. 1999, 144: 1991-2006. 10.1007/s007050050720.PubMedView ArticleGoogle Scholar
- Maruniak JE: Molecular biology of Anticarsia gemmatalis baculovirus. Mem Inst Oswaldo Cruz. 1989, 84: 107-111.View ArticleGoogle Scholar
- Oliveira JV, Wolff JL, Garcia-Maruniak A, Ribeiro BM, de Castro ME, de Souza ML, Moscardi F, Maruniak JE, Zanotto PM: Genome of the most widely used viral biopesticide: anticarsia gemmatalis multiple nucleopolyhedrovirus. J Gen Virol. 2006, 87: 3233-3250. 10.1099/vir.0.82161-0.PubMedView ArticleGoogle Scholar
- Lu A, Miller LK: Regulation of baculovirus late and very late gene expression. The Baculoviruses. Edited by: Miller LK. 1997, New York: Plenum, 193-216.View ArticleGoogle Scholar
- Happ B, Li J, Doerfler W: Proteins encoded in the 81.2- to 85.0-map-unit fragment of Autographa californica nuclear polyhedrosis virus DNA can be translated in vitro and in Spodoptera frugiperda cells. J Virol. 1991, 65: 89-97.PubMedPubMed CentralGoogle Scholar
- Oellig C, Happ B, Muller T, Doerfler W: Overlapping sets of viral RNAs reflect the array of polypeptides in the EcoRI J and N fragments (map positions 81.2 to 85.0) of the Autographa californica nuclear polyhedrosis virus genome. J Virol. 1987, 61: 3048-3057.PubMedPubMed CentralGoogle Scholar
- Smith I: Misleading messengers? Interpreting baculovirus transcriptional array profiles. J Virol. 2007, 81: 7819-7820. 10.1128/JVI.00615-07. author reply 7820–7811PubMedPubMed CentralView ArticleGoogle Scholar
- Friesen PD, Miller LK: Temporal regulation of baculovirus RNA: overlapping early and late transcripts. J Virol. 1985, 54: 392-400.PubMedPubMed CentralGoogle Scholar
- Kool M, Vlak JM: The structural and functional organization of the Autographa californica nuclear polyhedrosis virus genome. Arch Virol. 1993, 130: 1-16. 10.1007/BF01318992.PubMedView ArticleGoogle Scholar
- Lubbert H, Doerfler W: Transcription of overlapping sets of RNAs from the genome of Autographa californica nuclear polyhedrosis virus: a novel method for mapping RNAs. J Virol. 1984, 52: 255-265.PubMedPubMed CentralGoogle Scholar
- Mainprize TH, Lee K, Miller LK: Variation in the temporal expression of overlapping baculovirus transcripts. Virus Res. 1986, 6: 85-99. 10.1016/0168-1702(86)90059-6.PubMedView ArticleGoogle Scholar
- Katsuma S, Kang W, Shin-i T, Ohishi K, Kadota K, Kohara Y, Shimada T: Mass identification of transcriptional units expressed from the Bombyx mori nucleopolyhedrovirus genome. J Gen Virol. 2011, 92: 200-203. 10.1099/vir.0.025908-0.PubMedView ArticleGoogle Scholar
- Rodrigues JC, De Souza ML, O’Reilly D, Velloso LM, Pinedo FJ, Razuck FB, Ribeiro B, Ribeiro BM: Characterization of the ecdysteroid UDP-glucosyltransferase (egt) gene of Anticarsia gemmatalis nucleopolyhedrovirus. Virus Genes. 2001, 22: 103-112. 10.1023/A:1008142621359.PubMedView ArticleGoogle Scholar
- Razuck FB, Ribeiro B, Vargas JH, Wolff JL, Ribeiro BM: Characterization of the p10 gene region of Anticarsia gemmatalis nucleopolyhedrovirus. Virus Genes. 2002, 24: 243-247. 10.1023/A:1015328516018.PubMedView ArticleGoogle Scholar
- Slack JM, Shapiro M: Anticarsia gemmatalis multicapsid nucleopolyhedrovirus v-trex gene encodes a functional 3’ to 5’ exonuclease. J Gen Virol. 2004, 85: 2863-2871. 10.1099/vir.0.80109-0.PubMedView ArticleGoogle Scholar
- de Lima L, Pinedo FJ, Ribeiro BM, Zanotto PM, Wolff JL: Identification, expression and phylogenetic analysis of the Anticarsia gemmatalis multicapsid nucleopolyhedrovirus (AgMNPV) Helicase gene. Virus Genes. 2004, 29: 345-352. 10.1007/s11262-004-7438-8.PubMedView ArticleGoogle Scholar
- Carpes MP, de Castro ME, Soares EF, Villela AG, Pinedo FJ, Ribeiro BM: The inhibitor of apoptosis gene (iap-3) of Anticarsia gemmatalis multicapsid nucleopolyhedrovirus (AgMNPV) encodes a functional IAP. Arch Virol. 2005, 150: 1549-1562. 10.1007/s00705-005-0529-6.PubMedView ArticleGoogle Scholar
- Belaich MN, Rodriguez VA, Bilen MF, Pilloff MG, Romanowski V, Sciocco-Cap A, Ghiringhelli PD: Sequencing and characterisation of p74 gene in two isolates of Anticarsia gemmatalis MNPV. Virus Genes. 2006, 32: 59-70. 10.1007/s11262-005-5846-z.PubMedView ArticleGoogle Scholar
- Castro ME, Souza ML, Araujo S, Bilimoria SL: Replication of Anticarsia gemmatalis nuclear polyhedrosis virus in four lepidopteran cell lines. J Invertebr Pathol. 1997, 69: 40-45. 10.1006/jipa.1996.4624.PubMedView ArticleGoogle Scholar
- Sieburth P, Maruniak JE: Growth characteristics of a continuous cell line from the velvetbean caterpillar, Anticarsia gemmatalis hübner (lepidoptera: Noctuidae). In Vitro Cel Dev Biol Plant. 1988, 24: 195-198. 10.1007/BF02623546.View ArticleGoogle Scholar
- Vaughn JL, Goodwin RH, Tompkins GJ, McCawley P: The establishment of two cell lines from the insect Spodoptera frugiperda (Lepidoptera; Noctuidae). In vitro. 1977, 13: 213-217. 10.1007/BF02615077.PubMedView ArticleGoogle Scholar
- Davidson E, Levin M: Gene regulatory networks. Proc Natl Acad Sci USA. 2005, 102: 4935-10.1073/pnas.0502024102.PubMedPubMed CentralView ArticleGoogle Scholar
- Cohen D, Marek M, Davies B, Vlak JM, van Oers M: Encyclopedia of Autographa californica Nucleopolyhedrovirus Genes. Virol Sin. 2009, 24: 359-414. 10.1007/s12250-009-3059-7.View ArticleGoogle Scholar
- Baculovirus Molecular Biology: Baculovirus Molecular Biology. http://www.ncbi.nlm.nih.gov/books/NBK49500/,
- Wilson ME, Mainprize TH, Friesen PD, Miller LK: Location, transcription, and sequence of a baculovirus gene encoding a small arginine-rich polypeptide. J Virol. 1987, 61: 661-666.PubMedPubMed CentralGoogle Scholar
- Kelly DC, Brown DA DA, Ayres MD, Allen CJ, Walker IO: Properties of the Major Nucleocapsid Protein of Heliothis zea Singly Enveloped Nuclear Polyhedrosis Virus. J Gen Virol. 1983, 64: 399-408. 10.1099/0022-1317-64-2-399.View ArticleGoogle Scholar
- Iwanaga M, Takaya K, Katsuma S, Ote M, Tanaka S, Kamita SG, Kang W, Shimada T, Kobayashi M: Expression profiling of baculovirus genes in permissive and nonpermissive cell lines. Biochem Biophys Res Commun. 2004, 323: 599-614. 10.1016/j.bbrc.2004.08.114.PubMedView ArticleGoogle Scholar
- D’Costa SM, Yao H, Bilimoria SL: Transcription and temporal cascade in Chilo iridescent virus infected cells. Arch Virol. 2001, 146: 2165-2178. 10.1007/s007050170027.PubMedView ArticleGoogle Scholar
- Chinchar VG, Han J, Mao J, Brooks I, Srivastava K: Instability of frog virus 3 mRNA in productively infected cells. Virology. 1994, 203: 187-192. 10.1006/viro.1994.1473.PubMedView ArticleGoogle Scholar
- Chinchar VG, Yu W: Metabolism of host and viral mRNAs in frog virus 3-infected cells. Virology. 1992, 186: 435-443. 10.1016/0042-6822(92)90008-D.PubMedView ArticleGoogle Scholar
- Tjia ST, Carstens EB, Doerfler W: Infection of Spodoptera frugiperda cells with Autographa californica nuclear polyhedrosis virus II. The viral DNA and the kinetics of its replication. Virology. 1979, 99: 399-409. 10.1016/0042-6822(79)90018-7.PubMedView ArticleGoogle Scholar
- Liu HS, Bilimoria SL: Infected cell specific protein and viral DNA synthesis in productive and abortive infections of Spodoptera frugiperda nuclear polyhedrosis virus. Arch Virol. 1990, 115: 101-113. 10.1007/BF01310626.PubMedView ArticleGoogle Scholar
- Riegel CI, Slavicek JM: Characterization of the replication cycle of the Lymantria dispar nuclear polyhedrosis virus. Virus Res. 1997, 51: 9-17. 10.1016/S0168-1702(97)00075-0.PubMedView ArticleGoogle Scholar
- Moore RC, Purugganan MD: The early stages of duplicate gene evolution. Proc Natl Acad Sci USA. 2003, 100: 15682-15687. 10.1073/pnas.2535513100.PubMedPubMed CentralView ArticleGoogle Scholar
- Ohta T: Time for spreading of compensatory mutations under gene duplication. Genetics. 1989, 123: 579-584.PubMedPubMed CentralGoogle Scholar
- Blissard GW, Wenz JR: Baculovirus gp64 envelope glycoprotein is sufficient to mediate pH-dependent membrane fusion. J Virol. 1992, 66: 6829-6835.PubMedPubMed CentralGoogle Scholar
- Hoover K, Grove M, Gardner M, Hughes DP, McNeil J, Slavicek J: A gene for an extended phenotype. Science. 2011, 333: 1401-10.1126/science.1209199.PubMedView ArticleGoogle Scholar
- Foster KR: The sociobiology of molecular systems. Nat Rev Genet. 2011, 12: 193-203. 10.1038/nrg2903.PubMedView ArticleGoogle Scholar
- Fan TJ, Han LH, Cong RS, Liang J: Caspase family proteases and apoptosis. Acta Biochim Biophys Sin (Shanghai). 2005, 37: 719-727. 10.1111/j.1745-7270.2005.00108.x.View ArticleGoogle Scholar
- Viana MP, Tanck E, Beletti ME, Costa Lda F: Modularity and robustness of bone networks. Mol Biosyst. 2009, 5: 255-261. 10.1039/b814188f.PubMedView ArticleGoogle Scholar
- Jiang SS, Chang IS, Huang LW, Chen PC, Wen CC, Liu SC, Chien LC, Lin CY, Hsiung CA, Juang JL: Temporal transcription program of recombinant Autographa californica multiple nucleopolyhedrosis virus. J Virol. 2006, 80: 8989-8999. 10.1128/JVI.01158-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabasi AL: Hierarchical organization of modularity in metabolic networks. Science. 2002, 297: 1551-1555. 10.1126/science.1073374.PubMedView ArticleGoogle Scholar
- Hurst LD, Pal C, Lercher MJ: The evolutionary dynamics of eukaryotic gene order. Nat Rev Genet. 2004, 5: 299-310.PubMedView ArticleGoogle Scholar
- Lubbert H, Doerfler W: Mapping of early and late transcripts encoded by the Autographa californica nuclear polyhedrosis virus genome: is viral RNA spliced?. J Virol. 1984, 50: 497-506.PubMedPubMed CentralGoogle Scholar
- Jin J, Guarino LA: 3’-end formation of baculovirus late RNAs. J Virol. 2000, 74: 8930-8937. 10.1128/JVI.74.19.8930-8937.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Adhya S, Gottesman M: Promoter occlusion: transcription through a promoter may inhibit its activity. Cell. 1982, 29: 939-944. 10.1016/0092-8674(82)90456-1.PubMedView ArticleGoogle Scholar
- Dinger ME, Gascoigne DK, Mattick JS: The evolution of RNAs with multiple functions. Biochimie. 2011, 93: 2013-2018. 10.1016/j.biochi.2011.07.018.PubMedView ArticleGoogle Scholar
- Hughes TA: Regulation of gene expression by alternative untranslated regions. Trends Genet. 2006, 22: 119-122. 10.1016/j.tig.2006.01.001.PubMedView ArticleGoogle Scholar
- Friesen PD, Miller LK: Divergent transcription of early 35- and 94-kilodalton protein genes encoded by the HindIII K genome fragment of the baculovirus Autographa californica nuclear polyhedrosis virus. J Virol. 1987, 61: 2264-2272.PubMedPubMed CentralGoogle Scholar
- Ooi BG, Miller LK: Transcription of the baculovirus polyhedrin gene reduces the levels of an antisense transcript initiated downstream. J Virol. 1990, 64: 3126-3129.PubMedPubMed CentralGoogle Scholar
- Kang ST, Leu JH, Wang HC, Chen LL, Kou GH, Lo CF: Polycistronic mRNAs and internal ribosome entry site elements (IRES) are widely used by white spot syndrome virus (WSSV) structural protein genes. Virology. 2009, 387: 353-363. 10.1016/j.virol.2009.02.012.PubMedView ArticleGoogle Scholar
- Olszewski J, Miller LK: Identification and characterization of a baculovirus structural protein, VP1054, required for nucleocapsid formation. J Virol. 1997, 71: 5040-5050.PubMedPubMed CentralGoogle Scholar
- Summers MD, Smith GE: A manual of methods for baculovirus vectors and insect cell culture procedures. Texas Agric Exper Sta Bull. 1987, Issue 1555: 1-57.Google Scholar
- Abot AR, Moscardi F, Fuxa JR, Sosa-Gomez DR, Richter AR: Susceptibility of populations of Anticarsia gemmatalis from Brazil and the United States to a nuclear polyhedrosis virus. JEntomolSci. 1995, 30: 62-69.Google Scholar
- Abot AR, Moscardi F, Fuxa JR, Sosa-Gomez DR, Richter AR: Development of resistance by Anticarsia gemmatalis from Brazil and the United States to a nuclear polyhedrosis virus under laboratory selection pressure. Biol Control. 1996, 7: 126-130. 10.1006/bcon.1996.0075.View ArticleGoogle Scholar
- Allen GE, Knell JD: A nuclear polyhedrosis virus of Anticarsia gemmatalis: I, ultrastructure, replication, and pathogenicity. Fl Entomol. 1977, 60 (3): 233-240. 10.2307/3493914.View ArticleGoogle Scholar
- Fuxa JR, Richter AR: Distance and rate of spread of Anticarsia gemmatalis (Lepidoptera: Noctuidae) nuclear polyhedrosis virus released into soybean. Environ Entomol. 1994, 23: 1308-1316.View ArticleGoogle Scholar
- Johnson DW, Maruniak JE: Physical map of Anticarsia-gemmatalis nuclear polyhedrosis-virus (AgMNPV-2) DNA. J Gen Virol. 1989, 70: 1877-1883. 10.1099/0022-1317-70-7-1877.View ArticleGoogle Scholar
- Klein D: Quantification using real-time PCR technology: applications and limitations. Trends Mol Med. 2002, 8: 257-260. 10.1016/S1471-4914(02)02355-9.PubMedView ArticleGoogle Scholar
- Mackay IM, Arden KE, Nitsche A: Real-time PCR in virology. Nucleic Acids Res. 2002, 30: 1292-1305. 10.1093/nar/30.6.1292.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang D-H, Barari M, Arif BM, Krell PJ: Development of an oligonucleotide-based DNA microarray for transcriptional analysis of Choristoneura fumiferana nucleopolyhedrovirus (CfMNPV) genes. J Virol Methods. 2007, 143: 175-185. 10.1016/j.jviromet.2007.03.007.PubMedView ArticleGoogle Scholar
- Friesen PD: Regulation of baculovirus early gene expression. The Baculoviruses. Edited by: Miller LK. 1997, New York: Plenum, 141-170.View ArticleGoogle Scholar
- Guarino LA, Xu B, Jin J, Dong W: A virus-encoded RNA polymerase purified from baculovirus-infected cells. J Virol. 1998, 72: 7985-7991.PubMedPubMed CentralGoogle Scholar
- Carver T, Berriman M, Tivey A, Patel C, Bohme U, Barrell BG, Parkhill J, Rajandream MA: Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. Bioinformatics. 2008, 24: 2672-2676. 10.1093/bioinformatics/btn529.PubMedPubMed CentralView ArticleGoogle Scholar
- Kurakin A: Self-organization vs Watchmaker: stochastic gene expression and cell differentiation. Dev Genes Evol. 2005, 215: 46-52. 10.1007/s00427-004-0448-7.PubMedView ArticleGoogle Scholar
- Ross IL, Browne CM, Hume DA: Transcription of individual genes in eukaryotic cells occurs randomly and infrequently. Immunol Cell Biol. 1994, 72: 177-185. 10.1038/icb.1994.26.PubMedView ArticleGoogle Scholar
- Larson DR: What do expression dynamics tell us about the mechanism of transcription?. Curr Opin Genet Dev. 2011, 21: 591-599. 10.1016/j.gde.2011.07.010.PubMedPubMed CentralView ArticleGoogle Scholar
- Maier T, Guell M, Serrano L: Correlation of mRNA and protein in complex biological samples. FEBS Lett. 2009, 583: 3966-3973. 10.1016/j.febslet.2009.10.036.PubMedView ArticleGoogle Scholar
- Schwanhausser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, Chen W, Selbach M: Global quantification of mammalian gene expression control. Nature. 2011, 473: 337-342. 10.1038/nature10098.PubMedView ArticleGoogle Scholar
- Yamagishi J, Isobe R, Takebuchi T, Bando H: DNA microarrays of baculovirus genomes: differential expression of viral genes in two susceptible insect cell lines. Arch Virol. 2003, 148: 587-597. 10.1007/s00705-002-0922-3.PubMedView ArticleGoogle Scholar
- Huang S: Cell lineage determination in state space: a systems view brings flexibility to dogmatic canonical rules. PLoS Biol. 2010, 8: e1000380-10.1371/journal.pbio.1000380.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang S, Eichler G, Bar-Yam Y, Ingber DE: Cell fates as high-dimensional attractor states of a complex gene regulatory network. Phys Rev Lett. 2005, 94: 128701-PubMedView ArticleGoogle Scholar
- Kauffman SA: The origins of order: self-organization and selection in evolution. 1993, New York & Oxford: Oxford Univsersity PressGoogle Scholar
- D’Haeseleer P, Liang S, Somogyi R: Genetic network inference: from co-expression clustering to reverse engineering. Bioinformatics. 2000, 16: 707-726. 10.1093/bioinformatics/16.8.707.PubMedView ArticleGoogle Scholar
- Huynh-Thu VA, Irrthum A, Wehenkel L, Geurts P: Inferring regulatory networks from expression data using tree-based methods. PLoS ONE. 2010, 5: e12776-10.1371/journal.pone.0012776.PubMedPubMed CentralView ArticleGoogle Scholar
- The R Project for Statistical Computing. available at http://www.r-project.org/
- Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T: Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2010, 27: 431-432.PubMedPubMed CentralView ArticleGoogle Scholar
- Lin YY, Qi Y, Lu JY, Pan X, Yuan DS, Zhao Y, Bader JS, Boeke JD: A comprehensive synthetic genetic interaction network governing yeast histone acetylation and deacetylation. Genes Dev. 2008, 22: 2062-2074. 10.1101/gad.1679508.PubMedPubMed CentralView ArticleGoogle Scholar
- Holme P: Metabolic robustness and network modularity: a model study. PLoS ONE. 2011, 6 (2): e16605-10.1371/journal.pone.0016605.PubMedPubMed CentralView ArticleGoogle Scholar
- Brandes U: A faster algorithm for betweenness centrality. Journal Math Sociol. 2001, 25: 163-177. 10.1080/0022250X.2001.9990249.View ArticleGoogle Scholar
- Shimbel A: Structural parameters of communication networks. Bull Math Biophys. 1953, 15: 501-507. 10.1007/BF02476438.View ArticleGoogle Scholar
- Zanotto PMA, Krakauer DC: Complete genome viral phylogenies suggests the concerted evolution of regulatory cores and accessory satellites. PLoS ONE. 2008, 3: e3500-10.1371/journal.pone.0003500.View ArticleGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation. Bioinformatics. 2000, 16: 944-945. 10.1093/bioinformatics/16.10.944.PubMedView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed CentralView ArticleGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 95-98.Google Scholar
- Fetching the COALESCE program. http://evolution.genetics.washington.edu/lamarc/coalesce.html,
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.