Selected cachaça yeast strains share a genomic profile related to traits relevant to industrial fermentation processes.

The isolation and selection of yeast strains to improve the quality of the cachaça-Brazilian Spirit-have been studied in our research group. Our strategy considers Saccharomyces cerevisiae as the predominant species involved in sugarcane juice fermentation and the presence of different stressors (osmolarity, temperature, ethanol content, and competition with other microorganisms). It also considers producing balanced concentrations of volatile compounds (higher alcohols and acetate and/or ethyl esters), flocculation capacity, and ethanol production. Since the genetic bases behind these traits of interest are not fully established, the whole genome sequencing of 11 different Saccharomyces cerevisiae strains isolated and selected from different places was analyzed to identify the presence of a specific genetic variation common to cachaça yeast strains. We have identified 20,128 single-nucleotide variants shared by all genomes. Of these shared variants, 37 were new variants (being six missenses), and 4,451 were identified as missenses. We performed a detailed functional annotation (using enrichment analysis, protein-protein interaction network analysis, and database and in-depth literature searches) of these new and missense variants. Many genes carrying these variations were involved in the phenotypes of flocculation, tolerance to fermentative stresses, and production of volatile compounds and ethanol. These results demonstrate the existence of a genetic profile shared by the 11 strains under study that could be associated with the applied selective strategy. Thus, this study points out genes and variants that may be used as molecular markers for selecting strains well suited to the fermentation process, including genetic improvement by genome editing, ultimately producing high-quality beverages and adding value.IMPORTANCEThis work demonstrates the existence of new genetic markers related to different phenotypes used to select yeast strains and mutations in genes directly involved in producing flavoring compounds and ethanol, and others related to flocculation and stress resistance.

Bacterial volatile organic compounds induce adverse ultrastructural changes and DNA damage to the sugarcane pathogenic fungus Thielaviopsis ethacetica.

Due to an increasing demand for sustainable agricultural practices, the adoption of microbial volatile organic compounds (VOCs) as antagonists against phytopathogens has emerged as an eco-friendly alternative to the use of agrochemicals. Here, we identified three Pseudomonas strains that were able to inhibit, in vitro, up to 80% of mycelial growth of the phytopathogenic fungus Thielaviopsis ethacetica, the causal agent of pineapple sett rot disease in sugarcane. Using GC/MS, we found that these bacteria produced 62 different VOCs, and further functional validation revealed compounds with high antagonistic activity to T. ethacetica. Transcriptomic analysis of the fungal response to VOCs indicated that these metabolites downregulated genes related to fungal central metabolism, such as those involved in carbohydrate metabolism. Interestingly, genes related to the DNA damage response were upregulated, and micro-FTIR analysis corroborated our hypothesis that VOCs triggered DNA damage. Electron microscopy analysis showed critical morphological changes in mycelia treated with VOCs. Altogether, these results indicated that VOCs hampered fungal growth and could lead to cell death. This study represents the first demonstration of the molecular mechanisms involved in the antagonism of sugarcane phytopathogens by VOCs and reinforces that VOCs can be a sustainable alternative for use in phytopathogen biocontrol.

Ethanol production process driving changes on industrial strains.

Ethanol production has key differences between the two largest producing countries of this biofuel, Brazil and the USA, such as feedstock source, sugar concentration and ethanol titers in industrial fermentation. Therefore, it is highly probable that these specificities have led to genome adaptation of the Saccharomyces cerevisiae strains employed in each process to tolerate different environments. In order to identify particular adaptations, in this work, we have compared the genomes of industrial yeast strains widely used to produce ethanol from sugarcane, corn and sweet sorghum, and also two laboratory strains as reference. The genes were predicted and then 4524 single-copy orthologous were selected to build the phylogenetic tree. We found that the geographic location and industrial process were shown as the main evolutionary drivers: for sugarcane fermentation, positive selection was identified for metal homeostasis and stress response genes, whereas genes involved in membrane modeling have been connected with corn fermentation. In addition, the corn specialized strain Ethanol Red showed an increased number of copies of MAL31, a gene encoding a maltose transporter. In summary, our work can help to guide new strain chassis selection for engineering strategies, to produce more robust strains for biofuel production and other industrial applications.

Comparative transcriptome analysis reveals different strategies for degradation of steam-exploded sugarcane bagasse by Aspergillus niger and Trichoderma reesei.

BACKGROUND: Second generation (2G) ethanol is produced by breaking down lignocellulosic biomass into fermentable sugars. In Brazil, sugarcane bagasse has been proposed as the lignocellulosic residue for this biofuel production. The enzymatic cocktails for the degradation of biomass-derived polysaccharides are mostly produced by fungi, such as Aspergillus niger and Trichoderma reesei. However, it is not yet fully understood how these microorganisms degrade plant biomass. In order to identify transcriptomic changes during steam-exploded bagasse (SEB) breakdown, we conducted a RNA-seq comparative transcriptome profiling of both fungi growing on SEB as carbon source. RESULTS: Particular attention was focused on CAZymes, sugar transporters, transcription factors (TFs) and other proteins related to lignocellulose degradation. Although genes coding for the main enzymes involved in biomass deconstruction were expressed by both fungal strains since the beginning of the growth in SEB, significant differences were found in their expression profiles. The expression of these enzymes is mainly regulated at the transcription level, and A. niger and T. reesei also showed differences in TFs content and in their expression. Several sugar transporters that were induced in both fungal strains could be new players on biomass degradation besides their role in sugar uptake. Interestingly, our findings revealed that in both strains several genes that code for proteins of unknown function and pro-oxidant, antioxidant, and detoxification enzymes were induced during growth in SEB as carbon source, but their specific roles on lignocellulose degradation remain to be elucidated. CONCLUSIONS: This is the first report of a time-course experiment monitoring the degradation of pretreated bagasse by two important fungi using the RNA-seq technology. It was possible to identify a set of genes that might be applied in several biotechnology fields. The data suggest that these two microorganisms employ different strategies for biomass breakdown. This knowledge can be exploited for the rational design of enzymatic cocktails and 2G ethanol production improvement.

Insights into the plant polysaccharide degradation potential of the xylanolytic yeast Pseudozyma brasiliensis.

In second-generation (2G) bioethanol production, plant cell-wall polysaccharides are broken down to release fermentable sugars. The enzymes of this process are classified as carbohydrate-active enzymes (CAZymes) and contribute substantially to the cost of biofuel production. A novel basidiomycete yeast species, Pseudozyma brasiliensis, was recently discovered. It produces an endo-ß-1,4-xylanase with a higher specific activity than other xylanases. This enzyme is essential for the hydrolysis of biomass-derived xylan and has an important role in 2G bioethanol production. In spite of the P. brasiliensis biotechnological potential, there is no information about how it breaks down polysaccharides. For the first time, we characterized the secretome of P. brasiliensis grown on different carbon sources (xylose, xylan, cellobiose and glucose) and also under starvation conditions. The growth and consumption of each carbohydrate and the activity of the CAZymes of culture supernatants were analyzed. The CAZymes found in its secretomes, validated by enzymatic assays, have the potential to hydrolyze xylan, mannan, cellobiose and other polysaccharides. The data show that this yeast is a potential source of hydrolases, which can be used for biomass saccharification.

Proteomic analyses of baculovirus Anticarsia gemmatalis multiple nucleopolyhedrovirus budded and occluded virus.

Baculoviruses infect insects, producing two distinct phenotypes during the viral life cycle: the budded virus (BV) and the occlusion-derived virus (ODV) for intra- and inter-host spread, respectively. Since the 1980s, several countries have been using Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) as a biological control agent against the velvet bean caterpillar, A. gemmatalis. The genome of AgMNPV isolate 2D (AgMNPV-2D) carries at least 152 potential genes, with 24 that possibly code for structural proteins. Proteomic studies have been carried out on a few baculoviruses, with six ODV and two BV proteomes completed so far. Moreover, there are limited data on virion proteins carried by AgMNPV-2D. Therefore, structural proteins of AgMNPV-2D were analysed by MALDI- quadrupole-TOF and liquid chromatography MS/MS. A total of 44 proteins were associated with the ODV and 33 with the BV of AgMNPV-2D. Although 38 structural proteins were already known, we found six new proteins in the ODV and seven new proteins carried by the AgMNPV-2D BV. Eleven cellular proteins that were found on several other enveloped viruses were also identified, which are possibly carried with the virion. These findings may provide novel insights into baculovirus biology and their host interaction. Moreover, our data may be helpful in subsequent applied studies aiming to improve AgMNPV use as a biopesticide and a biotechnology tool for gene expression or delivery.

Pseudozyma brasiliensis sp. nov., a xylanolytic, ustilaginomycetous yeast species isolated from an insect pest of sugarcane roots.

A novel ustilaginomycetous yeast isolated from the intestinal tract of an insect pest of sugarcane roots in Ribeirão Preto, São Paulo State, Brazil, represents a novel species of the genus Pseudozyma based on molecular analyses of the D1/D2 rDNA large subunit and the internal transcribed spacer (ITS1+ITS2) regions. The name Pseudozyma brasiliensis sp. nov. is proposed for this species, with GHG001(T) ( = CBS 13268(T) = UFMG-CM-Y307(T)) as the type strain. P. brasiliensis sp. nov. is a sister species of Pseudozyma vetiver, originally isolated from leaves of vetiver grass and sugarcane in Thailand. P. brasiliensis sp. nov. is able to grow well with xylan as the sole carbon source and produces high levels of an endo-1,4-xylanase that has a higher specific activity in comparison with other eukaryotic xylanases. This enzyme has a variety of industrial applications, indicating the great biotechnological potential of P. brasiliensis.

Molecular characterization of the Aspergillus nidulans fbxA encoding an F-box protein involved in xylanase induction.

The filamentous fungus Aspergillus nidulans has been used as a fungal model system to study the regulation of xylanase production. These genes are activated at transcriptional level by the master regulator the transcriptional factor XlnR and repressed by carbon catabolite repression (CCR) mediated by the wide-domain repressor CreA. Here, we screened a collection of 42 A. nidulans F-box deletion mutants grown either in xylose or xylan as the single carbon source in the presence of the glucose analog 2-deoxy-D-glucose, aiming to identify mutants that have deregulated xylanase induction. We were able to recognize a null mutant in a gene (fbxA) that has decreased xylanase activity and reduced xlnA and xlnD mRNA accumulation. The ΔfbxA mutant interacts genetically with creAd-30, creB15, and creC27 mutants. FbxA is a novel protein containing a functional F-box domain that binds to Skp1 from the SCF-type ligase. Blastp analysis suggested that FbxA is a protein exclusive from fungi, without any apparent homologs in higher eukaryotes. Our work emphasizes the importance of the ubiquitination in the A. nidulans xylanase induction and CCR. The identification of FbxA provides another layer of complexity to xylanase induction and CCR phenomena in filamentous fungi.

Assessing the intracellular primary metabolic profile of Trichoderma reesei and Aspergillus niger grown on different carbon sources.

Trichoderma reesei and Aspergillus niger are efficient biological platforms for the production of various industrial products, including cellulases and organic acids. Nevertheless, despite the extensive research on these fungi, integrated analyses of omics-driven approaches are still missing. In this study, the intracellular metabolic profile of T. reesei RUT-C30 and A. niger N402 strains grown on glucose, lactose, carboxymethylcellulose (CMC), and steam-exploded sugarcane bagasse (SEB) as carbon sources for 48 h was analysed by proton nuclear magnetic resonance. The aim was to verify the changes in the primary metabolism triggered by these substrates and use transcriptomics data from the literature to better understand the dynamics of the observed alterations. Glucose and CMC induced higher fungal growth whereas fungi grown on lactose showed the lowest dry weight. Metabolic profile analysis revealed that mannitol, trehalose, glutamate, glutamine, and alanine were the most abundant metabolites in both fungi regardless of the carbon source. These metabolites are of particular interest for the mobilization of carbon and nitrogen, and stress tolerance inside the cell. Their concomitant presence indicates conserved mechanisms adopted by both fungi to assimilate carbon sources of different levels of recalcitrance. Moreover, the higher levels of galactose intermediates in T. reesei suggest its better adaptation in lactose, whereas glycolate and malate in CMC might indicate activation of the glyoxylate shunt. Glycerol and 4-aminobutyrate accumulated in A. niger grown on CMC and lactose, suggesting their relevant role in these carbon sources. In SEB, a lower quantity and diversity of metabolites were identified compared to the other carbon sources, and the metabolic changes and higher xylanase and pNPGase activities indicated a better utilization of bagasse by A. niger. Transcriptomic analysis supported the observed metabolic changes and pathways identified in this work. Taken together, we have advanced the knowledge about how fungal primary metabolism is affected by different carbon sources, and have drawn attention to metabolites still unexplored. These findings might ultimately be considered for developing more robust and efficient microbial factories.

Bacterial volatile organic compounds (VOCs) promote growth and induce metabolic changes in rice.

Plant growth-promoting bacteria (PGPB) represent an eco-friendly alternative to reduce the use of chemical products while increasing the productivity of economically important crops. The emission of small gaseous signaling molecules from PGPB named volatile organic compounds (VOCs) has emerged as a promising biotechnological tool to promote biomass accumulation in model plants (especially Arabidopsis thaliana) and a few crops, such as tomato, lettuce, and cucumber. Rice (Oryza sativa) is the most essential food crop for more than half of the world's population. However, the use of VOCs to improve this crop performance has not yet been investigated. Here, we evaluated the composition and effects of bacterial VOCs on the growth and metabolism of rice. First, we selected bacterial isolates (IAT P4F9 and E.1b) that increased rice dry shoot biomass by up to 83% in co-cultivation assays performed with different durations of time (7 and 12 days). Metabolic profiles of the plants co-cultivated with these isolates and controls (without bacteria and non-promoter bacteria-1003-S-C1) were investigated via 1H nuclear magnetic resonance. The analysis identified metabolites (e.g., amino acids, sugars, and others) with differential abundance between treatments that might play a role in metabolic pathways, such as protein synthesis, signaling, photosynthesis, energy metabolism, and nitrogen assimilation, involved in rice growth promotion. Interestingly, VOCs from IAT P4F9 displayed a more consistent promotion activity and were also able to increase rice dry shoot biomass in vivo. Molecular identification by sequencing the 16S rRNA gene of the isolates IAT P4F9 and E.1b showed a higher identity with Serratia and Achromobacter species, respectively. Lastly, volatilomes of these and two other non-promoter bacteria (1003-S-C1 and Escherichia coli DH5α) were evaluated through headspace solid-phase microextraction coupled with gas chromatography-mass spectrometry. Compounds belonging to different chemical classes, such as benzenoids, ketones, alcohols, sulfide, alkanes, and pyrazines, were identified. One of these VOCs, nonan-2-one, was validated in vitro as a bioactive compound capable of promoting rice growth. Although further analyses are necessary to properly elucidate the molecular mechanisms, our results suggest that these two bacterial isolates are potential candidates as sources for bioproducts, contributing to a more sustainable agriculture.

Gene Co-expression Network Reveals Potential New Genes Related to Sugarcane Bagasse Degradation in Trichoderma reesei RUT-30.

The biomass-degrading fungus Trichoderma reesei has been considered a model for cellulose degradation, and it is the primary source of the industrial enzymatic cocktails used in second-generation (2G) ethanol production. However, although various studies and advances have been conducted to understand the cellulolytic system and the transcriptional regulation of T. reesei, the whole set of genes related to lignocellulose degradation has not been completely elucidated. In this study, we inferred a weighted gene co-expression network analysis based on the transcriptome dataset of the T. reesei RUT-C30 strain aiming to identify new target genes involved in sugarcane bagasse breakdown. In total, ~70% of all the differentially expressed genes were found in 28 highly connected gene modules. Several cellulases, sugar transporters, and hypothetical proteins coding genes upregulated in bagasse were grouped into the same modules. Among them, a single module contained the most representative core of cellulolytic enzymes (cellobiohydrolase, endoglucanase, ß-glucosidase, and lytic polysaccharide monooxygenase). In addition, functional analysis using Gene Ontology (GO) revealed various classes of hydrolytic activity, cellulase activity, carbohydrate binding and cation:sugar symporter activity enriched in these modules. Several modules also showed GO enrichment for transcription factor activity, indicating the presence of transcriptional regulators along with the genes involved in cellulose breakdown and sugar transport as well as other genes encoding proteins with unknown functions. Highly connected genes (hubs) were also identified within each module, such as predicted transcription factors and genes encoding hypothetical proteins. In addition, various hubs contained at least one DNA binding site for the master activator Xyr1 according to our in silico analysis. The prediction of Xyr1 binding sites and the co-expression with genes encoding carbohydrate active enzymes and sugar transporters suggest a putative role of these hubs in bagasse cell wall deconstruction. Our results demonstrate a vast range of new promising targets that merit additional studies to improve the cellulolytic potential of T. reesei strains and to decrease the production costs of 2G ethanol.

Draft Genome Sequence of Saccharomyces cerevisiae Barra Grande (BG-1), a Brazilian Industrial Bioethanol-Producing Strain.

Here, we present the draft genome sequence of Saccharomyces cerevisiae BG-1, a Brazilian industrial strain widely used for bioethanol production from sugarcane. The 11.7-Mb genome sequence consists of 216 scaffolds and harbors 5,607 predicted protein-coding genes.

Comparative Secretome Analysis of Trichoderma reesei and Aspergillus niger during Growth on Sugarcane Biomass.

BACKGROUND: Our dependence on fossil fuel sources and concern about the environment has generated a worldwide interest in establishing new sources of fuel and energy. Thus, the use of ethanol as a fuel is advantageous because it is an inexhaustible energy source and has minimal environmental impact. Currently, Brazil is the world's second largest producer of ethanol, which is produced from sugarcane juice fermentation. However, several studies suggest that Brazil could double its production per hectare by using sugarcane bagasse and straw, known as second-generation (2G) bioethanol. Nevertheless, the use of this biomass presents a challenge because the plant cell wall structure, which is composed of complex sugars (cellulose and hemicelluloses), must be broken down into fermentable sugar, such as glucose and xylose. To achieve this goal, several types of hydrolytic enzymes are necessary, and these enzymes represent the majority of the cost associated with 2G bioethanol processing. Reducing the cost of the saccharification process can be achieved via a comprehensive understanding of the hydrolytic mechanisms and enzyme secretion of polysaccharide-hydrolyzing microorganisms. In many natural habitats, several microorganisms degrade lignocellulosic biomass through a set of enzymes that act synergistically. In this study, two fungal species, Aspergillus niger and Trichoderma reesei, were grown on sugarcane biomass with two levels of cell wall complexity, culm in natura and pretreated bagasse. The production of enzymes related to biomass degradation was monitored using secretome analyses after 6, 12 and 24 hours. Concurrently, we analyzed the sugars in the supernatant. RESULTS: Analyzing the concentration of monosaccharides in the supernatant, we observed that both species are able to disassemble the polysaccharides of sugarcane cell walls since 6 hours post-inoculation. The sugars from the polysaccharides such as arabinoxylan and ß-glucan (that compose the most external part of the cell wall in sugarcane) are likely the first to be released and assimilated by both species of fungi. At all time points tested, A. niger produced more enzymes (quantitatively and qualitatively) than T. reesei. However, the most important enzymes related to biomass degradation, including cellobiohydrolases, endoglucanases, ß-glucosidases, ß-xylosidases, endoxylanases, xyloglucanases, and α-arabinofuranosidases, were identified in both secretomes. We also noticed that the both fungi produce more enzymes when grown in culm as a single carbon source. CONCLUSION: Our work provides a detailed qualitative and semi-quantitative secretome analysis of A. niger and T. reesei grown on sugarcane biomass. Our data indicate that a combination of enzymes from both fungi is an interesting option to increase saccharification efficiency. In other words, these two fungal species might be combined for their usage in industrial processes.

Draft Genome Sequence of Pseudozyma brasiliensis sp. nov. Strain GHG001, a High Producer of Endo-1,4-Xylanase Isolated from an Insect Pest of Sugarcane.

Here, we present the nuclear and mitochondrial genome sequences of Pseudozyma brasiliensis sp. nov. strain GHG001. P. brasiliensis sp. nov. is the closest relative of Pseudozyma vetiver. P. brasiliensis sp. nov. is capable of growing on xylose or xylan as a sole carbon source and has great biotechnological potential.

A recombinant Anticarsia gemmatalis MNPV harboring chiA and v-cath genes from Choristoneura fumiferana defective NPV induce host liquefaction and increased insecticidal activity.

One of the interesting features of Anticarsia gemmatalis multiple nucleopolyhedrovirus isolate 2D (AgMNPV-2D) genome is the absence of chitinase (chiA) and cathepsin (v-cath) genes. This characteristic may be responsible for the lack of liquefaction and melanization in A. gemmatalis larvae killed by AgMNPV-2D infection. This study aimed to test the hypothesis that CHIA and V-CATH proteins from Choristonera fumiferana DEF multiple nucleopolyhedrovirus (CfDEFNPV) are able to liquefy and melanize the cuticle of A. gemmatalis larvae infected by a recombinant AgMNPV containing chiA and v-cath genes inserted in its genome. A fragment from the CfDefNPV genome containing chiA and v-cath genes was inserted into the genome of AgMNPV-2D. The recombinant virus (vAgp2100Cf.chiA/v-cath) was purified and used to infect insect cells and larvae. Transcripts of v-cath and chiA genes were detected along the infection of insect cells by qRT-PCR, from early to late phases of infection. The analysis of A. gemmatalis larvae killed by vAgp2100Cf.chiA/v-cath infection confirmed the hypothesis proposed. The vAgp2100Cf.chiA/v-cath showed higher insecticidal activity against third instar A. gemmatalis larvae when compared to AgMNPV-2D. The mean time to death was also lower for the vAgp2100Cf.chiA/v-cath when compared to AgMNPV-2D at 10 days post infection. Occlusion body production was higher in A. gemmatalis larvae infected with vAgp2100Cf.chiA/v-cath when compared to AgMNPV-2D. Enzyme assays showed higher chitinase and cysteine protease activities in insect cells and insects infected with vAgp2100Cf.chiA/v-cath when compared to AgMNPV-2D. The introduction of chiA and v-cath genes into the genome of AgMNPV improves its insecticidal activity against A. gemmatalis larvae and this recombinant virus could be used as an alternative to the wild type virus to control this important insect pest.

Analysis of the genome of Spodoptera frugiperda nucleopolyhedrovirus (SfMNPV-19) and of the high genomic heterogeneity in group II nucleopolyhedroviruses.

The genome of the most virulent among 22 Brazilian geographical isolates of Spodoptera frugiperda nucleopolyhedrovirus, isolate 19 (SfMNPV-19), was completely sequenced and shown to comprise 132,565 bp and 141 open reading frames (ORFs). A total of 11 ORFs with no homology to genes in the GenBank database were found. Of those, four had typical baculovirus promoter motifs and polyadenylation sites. Computer-simulated restriction enzyme cleavage patterns of SfMNPV-19 were compared with published physical maps of other SfMNPV isolates. Differences were observed in terms of the restriction profiles and genome size. Comparison of SfMNPV-19 with the sequence of the SfMNPV isolate 3AP2 indicated that they differed due to a 1427 bp deletion, as well as by a series of smaller deletions and point mutations. The majority of genes of SfMNPV-19 were conserved in the closely related Spodoptera exigua NPV (SeMNPV) and Agrotis segetum NPV (AgseMNPV-A), but a few regions experienced major changes and rearrangements. Synthenic maps for the genomes of group II NPVs revealed that gene collinearity was observed only within certain clusters. Analysis of the dynamics of gene gain and loss along the phylogenetic tree of the NPVs showed that group II had only five defining genes and supported the hypothesis that these viruses form ten highly divergent ancient lineages. Crucially, more than 60 % of the gene gain events followed a power-law relation to genetic distance among baculoviruses, indicative of temporal organization in the gene accretion process.

Genome of the most widely used viral biopesticide: Anticarsia gemmatalis multiple nucleopolyhedrovirus.

The genome of Anticarsia gemmatalis multiple nucleopolyhedrovirus isolate 2D (AgMNPV-2D), which is the most extensively used virus pesticide in the world, was completely sequenced and shown to have 132 239 bp (G+C content 44.5 mol%) and to be capable of encoding 152 non-overlapping open reading frames (ORFs). Three ORFs were unique to AgMNPV-2D, one of which (ag31) had similarity to eukaryotic poly(ADP-ribose) polymerases. The lack of chiA and v-cath may explain some of the success and growth of the AgMNPV biological control programme, as it may explain the high recovery of polyhedra sequestered inside dead larvae in the field, which are collected and used for further application as biological pesticides in soybean fields. The genome organization was similar to that of the Choristoneura fumiferana defective MNPV (CfDefNPV). Most of the variation between the two genomes took place near highly repetitive regions, which were also closely associated with bro-coding regions. The separation of the NPVs into groups I and II was supported by: (i) a phenogram of the complete genomes of 28 baculovirus and Heliothis zea virus 1, (ii) the most parsimonious reconstruction of gene content along the phenograms and (iii) comparisons of genomic features. Moreover, these data also reinforced the notion that group I of the NPVs can be split further into the AgMNPV lineage (AgMNPV, CfDefNPV, Epiphyas postvittana NPV, Orgyia pseudotsugata MNPV and C. fumiferana MNPV), sharing eight defining genes, and the Autographa californica MNPV (AcMNPV) lineage (AcMNPV, Rachiplusia ou NPV and Bombyx mori NPV), sharing nine defining genes.