Advancing Eucalyptus genomics: identification and sequencing of lignin biosynthesis genes from deep-coverage BAC libraries
© Paiva et al; licensee BioMed Central Ltd. 2011
Received: 1 June 2010
Accepted: 4 March 2011
Published: 4 March 2011
Eucalyptus species are among the most planted hardwoods in the world because of their rapid growth, adaptability and valuable wood properties. The development and integration of genomic resources into breeding practice will be increasingly important in the decades to come. Bacterial artificial chromosome (BAC) libraries are key genomic tools that enable positional cloning of important traits, synteny evaluation, and the development of genome framework physical maps for genetic linkage and genome sequencing.
We describe the construction and characterization of two deep-coverage BAC libraries EG_Ba and EG_Bb obtained from nuclear DNA fragments of E. grandis (clone BRASUZ1) digested with Hind III and BstY I, respectively. Genome coverages of 17 and 15 haploid genome equivalents were estimated for EG_Ba and EG_Bb, respectively. Both libraries contained large inserts, with average sizes ranging from 135 Kb (Eg_Bb) to 157 Kb (Eg_Ba), very low extra-nuclear genome contamination providing a probability of finding a single copy gene ≥ 99.99%. Libraries were screened for the presence of several genes of interest via hybridizations to high-density BAC filters followed by PCR validation. Five selected BAC clones were sequenced and assembled using the Roche GS FLX technology providing the whole sequence of the E. grandis chloroplast genome, and complete genomic sequences of important lignin biosynthesis genes.
The two E. grandis BAC libraries described in this study represent an important milestone for the advancement of Eucalyptus genomics and forest tree research. These BAC resources have a highly redundant genome coverage (> 15×), contain large average inserts and have a very low percentage of clones with organellar DNA or empty vectors. These publicly available BAC libraries are thus suitable for a broad range of applications in genetic and genomic research in Eucalyptus and possibly in related species of Myrtaceae, including genome sequencing, gene isolation, functional and comparative genomics. Because they have been constructed using the same tree (E. grandis BRASUZ1) whose full genome is being sequenced, they should prove instrumental for assembly and gap filling of the upcoming Eucalyptus reference genome sequence.
Renowned for their fast growth, valuable wood properties and wide adaptability, Eucalyptus species are among the most planted hardwoods in the world. The genus Eucalyptus includes over 700 species , including the most planted species E. grandis and E. urophylla (section Transversaria), E. globulus (section Maidenaria) and E. camaldulensis (section Exsertaria), all belonging to the subgenus Symphyomyrtus. Native to Australia, these fast growing trees were rapidly introduced into India, France, Chile, Brazil, South Africa, and Portugal in the first quarter of the 1800 s  and were promptly adopted for plantation forestry. Nowadays, Eucalyptus and their hybrids are among the world's leading sources of woody biomass and are the main hardwoods used for pulpwood and timber. In particular, E. grandis is grown in tropical and subtropical regions and E. globulus in limited temperate regions of the world while their hybrids are among the most widely used for industrial plantations because of their rapid growth rate, their adaptability to diverse ecological conditions and their good quality wood fiber.
Because of their high commercial value, Eucalyptus species are major targets for genetic improvement. Nevertheless from a genetics perspective, Eucalyptus are still in the early stages of domestication when compared to crop species and this fact has important implications when applying genomics approaches to understand the structural and functional biology of the genome [3–5].
In the last 15 years, several studies have led to a better understanding of the Eucalyptus genome and the development of an important set of genetic/genomic tools, which can be used to enhance future breeding efforts. Eucalyptus sp. are diploid plants with a haploid chromosome number of 11 [6–8]. Grattapaglia and Bradshaw  estimated the haploid genome size of Eucalyptus species to range ranging from 370 to 700 million base pairs (Mbp), with the Symphyomyrtus species having on average a haploid genome size of 650 Mbp. The high level of genetic diversity, the ability to generate large progeny sets, the relatively small genome size and low proportion of repetitive DNA facilitated genetic linkage mapping in Eucalyptus. Most Eucalyptus linkage maps [3–5], have been constructed from highly heterozygous parents and segregating half-sib of full-sib families, with progenies up to 200 individuals often from different species. The development of co-dominant markers such as microsatellites allowed comparative mapping studies among Eucalyptus species. Extensive map synteny and co linearity were detected between species of the subgenus Symphomyrtus. QTL mapping in Eucalyptus has been applied to the identification of genetic loci associated with variation in biomass productivity, stem form, wood properties (wood density and composition, fiber traits, bark composition), vegetative propagation, biotic/abiotic stress responses, development, foliar chemistry, inbreeding depression, and transcript level [4, 5]. Expressed sequence tag (EST) catalogues are available [11–16] and private consortia generated thousands of other ESTs. The complete nucleotide sequence of the chloroplast genome from E. globulus is also available , and is similar to other angiosperms, with an inverted repeat (IR) separated by a large single copy (LCS) region, and a small single copy (SSC) region.
Currently in progress, the first draft of a whole genome shotgun assembly of the Eucalyptus genome (E. grandis, clone BRASUZ1; http://www.eucagen.org) is expected to be publicly available by mid 2010. This effort will represent a major achievement and can be used as the first Eucalyptus reference sequence for future genomic undertakings .
Bacterial artificial chromosome (BAC) libraries have served as an essential genomic tool to enable positional cloning of important traits, synteny evaluation, and the development of genome framework physical maps for genetic linkage and genome sequencing . A number of BAC libraries have been reported for woody plants, including forest trees such as Liriodendron tulipifera, Pinus pinaster[20, 21], Populus trichocarpa[22, 23] and P. tremuloides[23, 24]. The first E. grandis BAC library was built in Brazil in the context of the Genolyptus project, supplying full genomic clones of key genes as well as end-sequences used to derive microsatellite and STS markers to help anchoring selected BAC clones to the existing linkage maps . Unfortunately, due to intellectual property restrictions, this library in not publicly available and no detailed report has been published yet about this specific work.
In this study, we report the construction and characterization of two deep-coverage BAC libraries (15-17×) from E. grandis (genotype "BRASUZ1") whose genome is currently being sequenced, becoming the public reference genome for the genus Eucalyptushttp://www.eucagen.org. The libraries described here i.e. EG_Ba (constructed from Hind III restricted nuclear DNA fragments) and EG_Bb library (constructed from Bst YI-digested DNA), are the first large-insert DNA libraries publicly available for genus Eucalyptushttp://www.genome.arizona.edu/orders/. The Roche GS FLX technology was successfully used to sequence five BAC clones allowing us to report here the whole sequence of the E. grandis chloroplast genome as well as the genomic sequences of important lignin biosynthetic genes.
Results and discussion
BAC library characterization
Two genomic BAC libraries (EG_Ba and EG_Bb) were constructed using partially digested (Hind III or BstY I) and size-selected nuclear DNA isolated from E. grandis (genotype BRASUZ1) and the pAGIBAC1 cloning vector, as described in Materials and Methods. Two libraries were constructed, using different restriction enzymes to avoid a biased distribution of clones along the Eucalyptus genome [26–29]. Each library, EG_Ba (Hind III) and EG_Bb (BstY I), contains 73,728 robotically picked clones arrayed into 192 384-well microtiter plates.
BAC library screening
To characterize these BAC libraries and facilitate clone identification, we prepared high-density macroarrays on nylon filters from a subset of the libraries representing 8.5× and 7.5× of genome coverage for EG_Ba and EG_Bb libraries, respectively. Membranes were hybridized with a series of pooled probes representing the chloroplast and mitochondria genome as well as with probes derived from lignin biosynthesis related genes.
Contamination with extranuclear genomes was estimated at 0.55% of the total number of BAC clones in both libraries. BAC clones containing E. grandis chloroplast sequences represented about 0.48% of all BAC clones in both libraries, lower than the estimates for several other plant species libraries [32–37]. The mitochondrial genome was represented by 0.07% of our BAC clones; slightly higher as compared to the 0.03%, 0.012%, and 0.04% found in coffee , tomato , and banana  BAC libraries, respectively.
BAC libraries screening. cDNA probes used were either involved in the lignin biosynthetic pathway (PAL, C4 H, HCT, C3 H, CCoAOMT, CCR, F5H/CAld5H, CAD) or regulating this pathway (MYB1 and 2, RAC1).
Positive clones Array hybridization (a)
Ratio (b/a) %
Phenylalanine ammonia lyase (PAL)
Cinnamate 4-hydroxylase (C4H)
4 coumarate CoA ligase (4CL)
Hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyl transferase (HTC)
p-coumarate 3-hydroxylase (C3H)
Caffeoyl-CoA O- methyltransferase (CCoaOMT)
Cinnamoyl CoA reductase (CCR)
Ferulate 5 hydroxylase/coniferaldehyde 5 hydroxylase (F5H/CAld5H)
Caffeic acid/5- hydroxyoniferaldehyde O-methyltransferase (COMT)
Cinnamyl Alcohol dehydrogenase
Rho-related small GTP-binding protein
R2R3 MYB transcription factor
R2R3 MYB transcription factor
To remove false positives and target a single genetic locus, we performed an additional confirmation by PCR, using specific primer pairs designed from available E. gunnii available cDNA sequences. On average, the estimated proportions of 62% and 57% of the hybridization positive clones were confirmed by PCR screening for the EG_Ba and EG_Bb libraries, respectively (Table 1).
The results of the hybridization experiments compared to those obtained by PCR validation suggest that these genes may be present in duplicate or belong to multigene families present in the Eucalyptus genome, in agreement with the EST analysis of Rengel et al.  that found different unigene members for some of these genes.
Sequencing of selected BAC clones
Sequencing (454FLX) and assembly of selected BAC clones.
% in the assembly
Large Contigs (> 500Kb)
The shotgun sequencing of these BAC clones produced on average 29,491 high quality reads per clone sequenced with a mean read length of 261nt. These clones were sequenced to different levels of sequence coverage ranging from 44.6× to 62.1×.
On average 98% of these reads were used to assemble the full sequences of each clone into a minimum of five and a maximum of 36 contigs, for EG_Ba_35H24 and EG_Ba_18G23, respectively. The number of large contigs (> 500 bp) varied among clones from 4 (for clone EG_Ba_35H24) to 25 (for clone EG_Bb_94G18). These contigs were then reassembled into one (clone EG_Ba_2B15) to six (clone EG_Bb_94G18) scaffolds, allowing the reconstruction of the full sequence of all BAC clones except clone EG_Bb_94G18 (6 scaffolds). In this latter case, the presence of repetitive sequences prevented our ability to order and orient four of the six. Such a problem was already reported in barley by Wicker et al. . Despite the relatively restricted number of clones sequenced, our results suggest the feasibility of using 454 sequencing for rapid and cost-effective sequencing and assembly of Eucalyptus BAC clones. The increased length of the 454 reads currently achievable with the Titanium chemistry (expected size ~400-550 bp) should result in regions of high-quality finished genomic sequences. BAC clone sequences were deposited into GenBank (accession ID HM366540 to HM366544).
Characterization of genomic nuclear BAC clone sequences
RepeatMasker http://www.repeatmasker.org/ was used to estimate GC content and search for interspersed repeats and low complexity DNA sequences in the four BAC sequences (588,509 bp, 13 scaffolds). The four BAC clone sequences revealed a low number of transposable elements, and low percentage of low complexity sequences (Additional file 1) when compared to 10 to 35% in other plant genomes analyzed (The Arabidopsis Genome Initiative 2000; International Rice Genome Sequencing Project 2005). However, these estimates cannot be generalized to the entire E. grandis genome due to the small number of BAC clones analyzed as well as the ascertainment bias resulting from the selection of gene-containing BACs. Furthermore, it is known that the distribution of these repetitive elements vary along the genomes. Indeed, clone EG_Ba_2B15 presented the lowest repeat content (1.2% of 152,083 bp) while the highest levels were found for clones EG_Bb_94G18 (5.5% of 129,018 bp) and clone EG_Ba_11K15 (5.4% of 137,697 bp). Six putative retroelements were found within the scaffold sequences of clone EG_Ba_18G23 (1.2% total length, 3 LINES and 3 LTRs) whereas none were identified in clone EG_Ba_2B15. Four LTR class retroelements were found within the scaffold sequences of clone EG_Bb_94G18 (2.5% total length, 3 Ty1/Copia and one Gypsy/DIRS1) and EG_Ba_11K15 (3,6% total length, 4 Gypsy/DIRS1). Low complexity sequences covered 0.8%, to 1.9% of the BAC clones sequence scaffolds analysed
Using Sputnik (Abajian, 1994, http://www.cbib.u-bordeaux2.fr/pise/sputnik.html) a total of 88 microsatellites also called simple sequence repeats (SSR) were found within the BAC scaffold sequences and these can be developed as new genetic markers. SSR markers have been extensively used in linkage analysis and comparative mapping of Eucalyptus species[10, 46–48], genetic fingerprinting , population genetics , and for clonal fidelity assessment . One possible application of these new BAC-derived genetic markers could be the anchoring of a physical map to the available Eucalyptus genetic maps. Within the BAC clone sequence scaffolds we found 57 perfect Class I SSRs, that are more likely to be polymorphic  as SSR mutations tend to be positively correlated with SSR length . This class of SSRs was found to occur on average every 10.3 kb within the sequences of the four nuclear BAC clones sequenced, from a minimum of 7.2 kb [EG_Bb_94G18] to 15.3 kb [EG_Ba_11K15]. These differences could reflect the non-random distribution of SSR in plant genomes [52, 54, 55]. The frequency of ClassI SSRs observed in this study was similar to that observed by Mun et al.  for selected Medicago truncatula BAC clone sequences, but very low when compared to that observed in another tree species Populus, where SSRs occurred on average every 2.5 kb . However, one cannot discard the possibility that the existence of gaps generated by the presence of repetitive sequences, where 454 sequencing has trouble going through, and/or due to a low sequence coverage in the region, might potentially underestimate the SSR frequency. Furthermore, the frequency of SSRs within the BAC clones is also very low when compared to that observed by Ceresini et al. ) and Rabello et al.  that reported a frequency of one Class I SSR every 2.5-2.7 kb in Eucalyptus cDNA libraries.
Comparison of the structure and sequence of the CCR and CAD2 genes from E. grandis and E. gunnii
Cinnamoyl CoA reductase (CCR; EC 126.96.36.199), that catalyzes the conversion of cinnamoyl CoA esters to their corresponding cinnamaldehydes and Cinnamyl alcohol dehydrogenase (CAD; EC 1.1.195) that catalyzes the conversion of these aldehydes to the corresponding alcohols are considered key enzymes in lignin biosynthesis. For instance, it was previously found that CCR down-regulated plants had lower lignin levels than controls [58, 59] and the extractability of the lignin polymer was improved in CAD down-regulated plants [60, 61]. Sequences of positive BAC clones for CAD and CCR probes were analyzed for gene identification via homology searches. Searches with BLASTN performed against a non-redundant databases (e-value cut off of 1e-5) allowed us to easily identify the E. grandis homologous genes to the E. gunnii CCR gene (within scaffold #1 of EG_Ba_2B15) and the CAD2 gene (within scaffold#2 of clone EG_Ba_11K15).
Non-coding regions were also very well conserved between E. gunnii and E. grandis CAD and CCR sequences, respectively. The CAD promoter regions exhibited sequence identity of 96% for the first 2.5 kb. Whereas the sequence conservation was lower between the CCR promoter sequences, showing 88% sequence identity within the first 500 pb upstream the transcription start. The alignment was interrupted by an insertion of 444 bp in the E. grandis sequence and the identity level in the remaining 5' sequence dropped to 85%.
Concerning introns, identities between E. grandis and E. gunnii increased from 90% to 97%, the less conserved being intron 4, showing 90% and 93% identity for CAD2 and CCR respectively. Intron 1 in the E. grandis CAD2 gene exhibited a deletion of 71 nt after position 122 as compared to the corresponding intron in E. gunnii.
The indels, reported here, between the CCR promoters and between introns 1 of CAD2 could be used to develop markers to discriminate these two species, as they seem highly conserved within each species (Additional file 2).
The E. grandis chloroplast genome characterization
The two Eucalyptus BAC libraries described in this study represent an important milestone for the advancement of Eucalyptus genomics and forest tree research. These BAC resources have a highly redundant genome coverage (> 15×), contain large average inserts (157 and 135 kb) and have a very low percentage of clones with organellar DNA or empty vectors. This indicates that these publicly available BAC libraries are suitable for a broad range of applications in genetic and genomic research in Eucalyptus and possibly in related species of Myrtaceae, including genome sequencing, gene isolation, functional and comparative genomics. The analysis of ~0.6 Mb of BAC clone sequences generated by Roche GS FLX sequencing technology provided an overview of the composition of the Eucalyptus nuclear genome and the feasibility of using this high-throughput technology for low-cost and efficient sequencing and assembly of the targeted Eucalyptus sequences. SSRs identified within the BAC clone sequences could be used to develop new genetic markers for multiple genotyping purposes. In addition, we report the full chloroplast genome sequence of E. grandis (160,137 bp) allowing comparison of this genome with E. globulus and other plant species. Comparative analysis of the CAD2 and CCR genes between E. grandis and E. gunnii showed a high conservation of the structure of genes as well as a high identity both in the coding and non coding sequences.
Plant material and DNA preparation
High Molecular Weight (HMW) DNA was prepared from young leaves of clonally propagated plants of tree BRASUZ1 grown in partial shade at the Federal University of Viçosa, Brazil. BRASUZ1 is a S1 individual (one generation of selfing) of an elite E. grandis tree originally derived from seed lots collected in Coffs Harbor (Australia). BRASUZ1 was developed by Suzano Papel e Celulose Co. in 1987 and confirmed as S1 by microsatellite genotyping (D. Grattapaglia, unpublished). It is a fast growing tree that flowers normally and does not exhibit any sign of inbreeding depression. This tree is being used by the Eucalyptus Genome project currently under production at JGI (Joint Genome Institute) for its lower genomic heterozygosity as compared to a regular outcrossed tree to facilitate assembly. Tender, expanded leaves were collected during a time period of two months and kept frozen at -80°C. For each extraction, approximately 50 grams of frozen leaf tissue was ground to powder in liquid nitrogen with a mortar and pestle used to prepare megabase-size DNA embedded in agarose plugs as described by Zhang et al.  Agarose plugs containing high molecular weight (HMW) nuclei DNA were subsequently sent submerged in ethanol 95% to the Arizona Genomic Institute (AGI) for library construction.
BAC library construction and BAC clone characterization
BAC libraries were constructed  using modifications recently described for Oryza. DNA digestions were performed with 0.5 Unit Hind III in 100 μL reaction volume (EG_Ba library) and 0.8 Unit Bst YI also in 100 μL reaction volumeFor (EG_Bb libary) to obtain the appropriate partial digestion conditions. the EG_Ba BAC library, pulsed-field gel electrophoresis (PFGE) size-selected restriction fragments were ligated to the pAGIBAC1 vector (a modified pIndigoBAC536Blue with an additional Swa I site), while for the EG_Bb library, size-selected BstY I digested fragments were ligated to BamHI digested vector. Ligation products were transformed into DH10B T1 phage resistant Escherichia coli cells (Invitrogen, Carlsbad, CA) and plated on LB agar that contained chloramphenicol (12.5 μg.mL-1), X-gal (20 mg.mL-1) and IPTG (0.1 M). Clones were robotically picked into 384-well plates containing LB freezing media. Plates were incubated for 18 h, replicated and then frozen at -80°C for long term storage.
To estimate insert sizes, 5 μL aliquots of extracted BAC DNA were digested with 5 U of Not I enzyme for 3 hrs at 37°C. The digestion products were resolved by PFGE (CHEF-DRIII system, Bio-Rad) in a 1% agarose gel in TBE buffer 0.5×. Electrophoresis was carried out for 16 hours at 14°C with an initial switch time of 5 sec, a final switch time of 15 sec, in a voltage gradient of 6 V.cm-1. Insert sizes were compared to those of the MidRange I PFG Marker (New England Biolabs).
BAC library screening
a) High Density Filter production and hybridization
A subset (plates 1-96) of libraries EG_Ba and EG_Bb was used for screening. High density colony filters for both libraries were prepared using a Genetix Q-bot (Genetix, New Milton, Hampshire, United Kingdom). Each 22.5×22.5 cm filter (Hybond-N+; Amersham Bioscience, Piscataway, NJ, USA) contained 36,864 independent clones arrayed in a double spotted 6×6 pattern. Spotted nylon membranes were transferred to solid LB agar plates supplemented with 12.5 μg/mL chloramphenicol, and bacterial clones allowed to grow overnight at 37°C, followed by transferred to 4°C just until colonies became confluent. Membranes were kept for 4 min onto a Whatmann 3 MM paper saturated with denaturation buffer (0.5 M NaOH, 1.5N NaCl), treated for 10 min at 100°C, and neutralized for 10 min on Whatmann 3 MM paper saturated with neutralization buffer (1.5 M Tris-HCL pH7.4, 1.5N NaCl). Immediately membranes were incubated at 37°C for 45 min with 250 mg/L proteinase K in 100 mM Tris-HCl pH 8, 50 mM EDTA, 0.5N NaCl. Finally, membranes were dried for 45 min at 80°C and UV-crosslinked (120,000 μJ.cm-2 for 50 sec).
Labelling of organellar or nuclear probes was performed with [α-33P] dCTP (Perkin-Elmer, Waltham, Massachusetts, USA) on DNA (150 ng amplified DNA fragment) using the Ready-To-GO DNA Labeling kit. The dCTP (GE Healthcare, Waukesha, WI) and the non-incorporated nucleotides were removed using Illustra-ProbeQuant G50 (GE Healthcare, Waukesha, WI) according to manufacture's instructions.
Before performing hybridization, membranes were incubated for 30 min at 50°C in 6× SSC. Filters were pre-hybridized for 2 hours in 50 mL of hybridization buffer (6XSSC, 5X Denhardt, 0.5% SDS, 100 μg.mL-1 denatured salmon sperm DNA) at 68°C. Hybridizations were carried out in high stringency conditions at 68°C overnight using 50 mL of fresh hybridization buffer supplemented with a minimum of 107cpm purified and denatured probe per mL of buffer. Membranes were washed for 15 min at 50°C in SSC2x and SDS 0.1% buffer, followed by a second wash at 50°C for 30 min in SSC 0.5X and SDS 0.1% buffer. Finally, they were wrapped in a plastic film, exposed to the General Purpose PhosphorImager screen (Amersham Bioscience, Piscataway, NJ, USA) for a period of three days, and finally scanned using a Storm System (Amersham Bioscience (Amersham Bioscience, Piscataway, NJ, USA), set to a resolution of 50 μm.
b) Organellar DNA content estimation
To estimate the percentage of chloroplast and mitochondrial DNA content in each library, one high-density filter from each library was screened separately with a pool of five E. globulus probes for chloroplast genes psbA, psbB, psbD, rbcL, ndhB and with a pool of probes for mitochondrial genes ccb256, ccb452, cox3. Chloroplast probes were obtained by amplification of E. globulus DNA using specific chloroplast primers (Additional file 3). Positive clones for chloroplast and mitochondria probes were identified and 10% of these clones were used in the PCR validation step. BAC clone (EG_Ba_35H24] was selected to be full sequenced, and was subsequently sequenced by Roche GS FLX technology (454 Life Sciences, Branford, CT, USA) by Cogenics company(Grenoble, France).
c) Screening for selected genes
Each library was screened with different probes for genes of lignification biosynthesis pathway and for three regulatory genes (EguRAC1; EguMyb1, EguMyb2). Characteristics of each gene probe are described in Table 1. Clone inserts were amplified using universal primers M13 (5'CAGGAAACAGCTATGACC3') and M13 reverse (5'TGTAAAACGACGGCCAGT3'), and the identity of the amplicons checked by sequencing before use in the hybridization. Hybridization positive clones were then validated by PCR using specific gene primers (Additional file 4).
BAC clones sequencing
Once selected, BAC clones were isolated, and the presence of the gene of interest was re-confirmed by PCR on isolated colony. BAC insert size was determined as described above. BAC DNA extractions, sequencing library generation, and sequencing were performed by Cogenics (Grenoble, France). BAC insert were sequenced by pyrosequencing using a Roche GS FLX Life Sciences instrument (Branford, CT, USA) [70, 71]. The Newbler software (454 Life Sciences, Branford, CT, USA) was used to perform de novo assembly of the Roche GS FLX. 454 These assemblies may still contain indel errors even at high coverage.
E. grandis chloroplast genome sequence annotation
The E. grandis chloroplast genome was annotated using the E. globulus chloroplast genome (NCBI accession id NC_008115) as the reference genome and following the annotation pipeline detailed in Addicional file 5.
Since these two genomes correspond to circular genomes, it was not guaranteed that the linear sequence obtained from sequencing was cut at the same genome position. An important step before using Multiple Sequence Alignment (MSA) algorithms on genome annotation of circular genomes includes the circularization and rotation of the genomes that are to be aligned. MSA algorithms were developed to deal with linear genomic sequences and in this sense they are very sensitive to the location where the genomic sequence begins. To improve the MSA alignment of both E. grandis and E. globulus, these two genomes were first processed using the CSA algorithm . This algorithm identified the largest chain of non-repeated longest subsequences common to these two circular genomes. The genomes were then rotated and made linear for MSA purposes.
After the pre-processing step, the Dual Organellar GenoMe Annotator (DOGMA) software package  was used to perform BLAST searches against a custom database of plant chloroplast genomes.
To conclude the annotation procedure, the Artemis Comparison Tool (ACT)  was used. With this tool it was possible to visualize and solve any the remaining annotation inconsistencies, by comparing the annotation obtained from DOGMA with the E. globulus chloroplast genome sequence. Each gene structure was manually checked to define putative exons.
This work was partially supported by Fundação para a Ciência e Tecnologia (Portugal) (GenEglobwq project - PTDC/AGR-GPL/66564/2006). Authors also acknowledge the fundings from the ERA-PG EUCANET project n° ANR-06-ERAPG-10-03, from the Brazilian Ministry of Science and Technology (GENOLYPTUS project - FINEP grant 1755-01 and CNPq grant 520489/02-0) and a CNPq research fellowship to DG. The authors gratefully acknowledge Nathalie Ladouce for help on cDNA clones selection and amplification for probe preparation.
- Brooker MIH: A new classification of the genus Eucalyptus L'Her (Myrtaceae). Australian Systematic Botany. 2000, 13 (1): 79-148. 10.1071/SB98008.View ArticleGoogle Scholar
- Doughty RW: The Eucalyptus. A natural and commercial history of the gum tree. 2000, Baltimore USA, London UK: Johns Hopkins Universitary PressGoogle Scholar
- Poke FS, Vaillancourt RE, Potts BM, Reid JB: Genomic research in Eucalyptus. Genetica. 2005, 125 (1): 79-101. 10.1007/s10709-005-5082-4.PubMedView ArticleGoogle Scholar
- Myburg A, Potts B, Marques C, Kirst M, Gion J, Grattapaglia D, Grima-Pettenati J:Eucalyptus. Genome mapping and molecular breeding in plants. Edited by: C K. 2007, New York, NY, USA: Springer, 7: 115-160. full_text.Google Scholar
- Grattapaglia D, Kirst M: Eucalyptus applied genomics: from gene sequences to breeding tools. New Phytologist. 2008, 179 (4): 911-929. 10.1111/j.1469-8137.2008.02503.x.PubMedView ArticleGoogle Scholar
- Eldridge K, Davidson J, Harwood C, van Wyk G: Eucalypt domestication and breeding. 1993, Oxford, UK: Clarendon PressGoogle Scholar
- Potts BM, Wiltshire RJE: Eucalypt genetics and genecology. Eucalypt ecology: individuals to ecosystems. Cambridge, UK. Edited by: Williams J, Woinarski J. 1997, 56-91.Google Scholar
- Grattapaglia D, Bradshaw HD: Nuclear-DNA content of commercially important Eucalyptus species and hybrids. Canadian Journal of Forest Research-Revue Canadienne De Recherche Forestiere. 1994, 24 (5): 1074-1078. 10.1139/x94-142.View ArticleGoogle Scholar
- Grattapaglia D, Sederoff R: Genetic-linkage maps of Eucalyptus grandis and Eucalyptus urophylla using pseudo testcross mapping strategy and RAPD markers. Genetics. 1994, 137 (4): 1121-1137.PubMedPubMed CentralGoogle Scholar
- Marques CM, Brondani RPV, Grattapaglia D, Sederoff R: Conservation and synteny of SSR loci and QTLs for vegetative propagation in four Eucalyptus species. Theoretical and Applied Genetics. 2002, 105 (2-3): 474-478. 10.1007/s00122-002-0899-z.PubMedView ArticleGoogle Scholar
- Foucart C, Paux E, Ladouce N, San-Clemente H, Grima-Pettenati J, Sivadon P: Transcript profiling of a xylem vs phloem cDNA subtractive library identifies new genes expressed during xylogenesis in Eucalyptus. New Phytologist. 2006, 170 (4): 739-752. 10.1111/j.1469-8137.2006.01705.x.PubMedView ArticleGoogle Scholar
- Paux E, Tamasloukht M, Ladouce N, Sivadon P, Grima-Pettenati J: Identification of genes preferentially expressed during wood formation in Eucalyptus. Plant Molecular Biology. 2004, 55 (2): 263-280. 10.1007/s11103-004-0621-4.PubMedView ArticleGoogle Scholar
- Paux E, Carocha V, Marques C, de Sousa AM, Borralho N, Sivadon P, Grima-Pettenati J: Transcript profiling of Eucalyptus xylem genes during tension wood formation. New Phytologist. 2005, 167 (1): 89-100. 10.1111/j.1469-8137.2005.01396.x.PubMedView ArticleGoogle Scholar
- Rengel D, Clemente HS, Servant F, Ladouce N, Paux E, Wincker P, Couloux A, Sivadon P, Grima-Pettenati J: A new genomic resource dedicated to wood formation in Eucalyptus. BMC Plant Biology. 2009, 9-Google Scholar
- Rasmussen-Poblete S, Valdes J, Gamboa MC, Valenzuela PDT, Krauskopf E: Generation and analysis of an Eucalyptus globulus cDNA library constructed from seedlings subjected to low temperature conditions. Electronic Journal of Biotechnology. 2008, 11 (2): 10.2225/vol11-issue2-fulltext-14.
- Keller G, Marchal T, SanClemente H, Navarro M, Ladouce N, Wincker P, Couloux A, Teulieres C, Marque C: Development and functional annotation of an 11,303-EST collection from Eucalyptus for studies of cold tolerance. Tree Genetics & Genomes. 2009, 5 (2): 317-327.View ArticleGoogle Scholar
- Steane DA: Complete nucleotide sequence of the chloroplast genome from the Tasmanian blue gum, Eucalyptus globulus (Myrtaceae). DNA Research. 2005, 12 (3): 215-220. 10.1093/dnares/dsi006.PubMedView ArticleGoogle Scholar
- Monaco AP, Larin Z: YACS, BACS, PACS and MACS - artificial chromosomes. Trends in Biotechnology. 1994, 12 (7): 280-286. 10.1016/0167-7799(94)90140-6.PubMedView ArticleGoogle Scholar
- Liang HY, Fang EG, Tomkins JP, Luo MZ, Kudrna D, Kim HR, Arumuganathan K, Zhao SY, Leebens-Mack J, Schlarbaum SE, et al: Development of a BAC library for yellow-poplar (Liriodendron tulipifera) and the identification of genes associated with flower development and lignin biosynthesis. Tree Genetics & Genomes. 2007, 3 (3): 215-225.View ArticleGoogle Scholar
- Bautista R, Villalobos DP, Diaz-Moreno S, Canton FR, Canovas FM, Claros MG: Toward a Pinus pinaster bacterial artificial chromosome library. Annals of Forest Science. 2007, 64 (8): 855-864. 10.1051/forest:2007060.View ArticleGoogle Scholar
- Bautista R, Villalobos DP, Diaz-Moreno S, Canton ER, Canovas EM, Claros MG: New strategy for Pinus pinaster genomic library construction in bacterial artificial chromosomes. Investigacion Agraria-Sistemas Y Recursos Forestales. 2008, 17 (3): 238-249.View ArticleGoogle Scholar
- Stirling B, Newcombe G, Vrebalov J, Bosdet I, Bradshaw HD: Suppressed recombination around the MXC3 locus, a major gene for resistance to poplar leaf rust. Theoretical and Applied Genetics. 2001, 103 (8): 1129-1137. 10.1007/s001220100721.View ArticleGoogle Scholar
- Tuskan GA, Gunter LE, Yang ZMK, Yin TM, Sewell MM, DiFazio SP: Characterization of microsatellites revealed by genomic sequencing of Populus trichocarpa. Canadian Journal of Forest Research-Revue Canadienne De Recherche Forestiere. 2004, 34 (1): 85-93. 10.1139/x03-283.View ArticleGoogle Scholar
- Fladung M, Kaufmann H, Markussen T, Hoenicka H: Construction of a Populus tremuloides Michx. BAC library. Silvae Genetica. 2008, 57 (2): 65-69.Google Scholar
- Grattapaglia D, Alfenas A, Coelho A, Bearzoti E, Pappas G, Pasquali G, Pereira G, Colodette J, Gomide J, Bueno J, et al: Building resources for molecular breeding of Eucalyptus. International IUFRO Conference - Eucalyptus in a changing world: 2004: 2004; Aveiro, Portugal. 2004, 20-32.Google Scholar
- Zhang HB, Choi SD, Woo SS, Li ZK, Wing RA: Construction and characterization of two rice bacterial artificial chromosome libraries from the parents of a permanent recombinant inbred mapping population. Molecular Breeding. 1996, 2 (1): 11-24. 10.1007/BF00171348.View ArticleGoogle Scholar
- Chen Q, Sun S, Ye Q, McCuine S, Huff E, Zhang HB: Construction of two BAC libraries from the wild Mexican diploid potato, Solanum pinnatisectum, and the identification of clones near the late blight and Colorado potato beetle resistance loci. Theoretical and Applied Genetics. 2004, 108 (6): 1002-1009. 10.1007/s00122-003-1513-8.PubMedView ArticleGoogle Scholar
- Tao Q, Wang A, Zhang HB: One large-insert plant-transformation-competent BIBAC library and three BAC libraries of Japonica rice for genome research in rice and other grasses. Theoretical and Applied Genetics. 2002, 105 (6-7): 1058-1066. 10.1007/s00122-002-1057-3.PubMedView ArticleGoogle Scholar
- Wu CC, Nimmakayala P, Santos FA, Springman R, Scheuring C, Meksem K, Lightfoot DA, Zhang HB: Construction and characterization of a soybean bacterial artificial chromosome library and use of multiple complementary libraries for genome physical mapping. Theoretical and Applied Genetics. 2004, 109 (5): 1041-1050. 10.1007/s00122-004-1712-y.PubMedView ArticleGoogle Scholar
- Ammiraju JSS, Luo MZ, Goicoechea JL, Wang WM, Kudrna D, Mueller C, Talag J, Kim H, Sisneros NB, Blackmon B, et al: The Oryza bacterial artificial chromosome library resource: Construction and analysis of 12 deep-coverage large-insert BAC libraries that represent the 10 genome types of the genus Oryza. Genome Research. 2006, 16 (1): 140-147. 10.1101/gr.3766306.PubMedPubMed CentralView ArticleGoogle Scholar
- Clarke L, Carbon J: Colony bank containing synthetic plasmids representative of of entire Escherichia coli genome. Cell. 1976, 9 (1): 91-99. 10.1016/0092-8674(76)90055-6.PubMedView ArticleGoogle Scholar
- Luo MZ, Wang YH, Frisch D, Joobeur T, Wing RA, Dean RA: Melon bacterial artificial chromosome (BAC) library construction using improved methods and identification of clones linked to the locus conferring resistance to melon Fusarium wilt (Fom-2). Genome. 2001, 44 (2): 154-162. 10.1139/gen-44-2-154.PubMedView ArticleGoogle Scholar
- Leroy T, Marraccini P, Dufour M, Montagnon C, Lashermes P, Sabau X, Ferreira LP, Jourdan I, Pot D, Andrade AC, et al: Construction and characterization of a Coffea canephora BAC library to study the organization of sucrose biosynthesis genes. Theoretical and Applied Genetics. 2005, 111 (6): 1032-1041. 10.1007/s00122-005-0018-z.PubMedView ArticleGoogle Scholar
- Stevens MR, Coleman CE, Parkinson SE, Maughan PJ, Zhang HB, Balzotti MR, Kooyman DL, Arumuganathan K, Bonifacio A, Fairbanks DJ, et al: Construction of a quinoa (Chenopodium quinoa Willd.) BAC library and its use in identifying genes encoding seed storage proteins. Theoretical and Applied Genetics. 2006, 112 (8): 1593-1600. 10.1007/s00122-006-0266-6.PubMedView ArticleGoogle Scholar
- Cavagnaro PF, Chung SM, Szklarczyk M, Grzebelus D, Senalik D, Atkins AE, Simon PW: Characterization of a deep-coverage carrot (Daucus carota L.) BAC library and initial analysis of BAC-end sequences. Molecular Genetics and Genomics. 2009, 281 (3): 273-288. 10.1007/s00438-008-0411-9.PubMedView ArticleGoogle Scholar
- Terol J, Naranjo MA, Ollitrault P, Talon M: Development of genomic resources for Citrus clementina: Characterization of three deep-coverage BAC libraries and analysis of 46,000 BAC end sequences. BMC Genomics. 2008, 9-Google Scholar
- Nilmalgoda SD, Cloutier S, Walichnowski AZ: Construction and characterization of a bacterial artificial chromosome (BAC) library of hexaploid wheat (Triticum aestivum L.) and validation of genome coverage using locus-specific primers. Genome. 2003, 46 (5): 870-878. 10.1139/g03-067.PubMedView ArticleGoogle Scholar
- Hamilton CM, Frary A, Xu YM, Tanksley SD, Zhang HB: Construction of tomato genomic DNA libraries in a binary-BAC (BIBAC) vector. Plant Journal. 1999, 18 (2): 223-229. 10.1046/j.1365-313X.1999.00433.x.View ArticleGoogle Scholar
- Safar J, Noa-Carrazana JC, Vrana J, Bartos J, Alkhimova O, Sabau X, Simkova H, Lheureux F, Caruana ML, Dolezel J, et al: Creation of a BAC resource to study the structure and evolution of the banana (Musa balbisiana) genome. Genome. 2004, 47 (6): 1182-1191. 10.1139/g04-062.PubMedView ArticleGoogle Scholar
- Lacombe E, Hawkins S, VanDoorsselaere J, Piquemal J, Goffner D, Poeydomenge O, Boudet AM, GrimaPettenati J: Cinnamoyl CoA reductase, the first committed enzyme of the lignin branch biosynthetic pathway: Cloning, expression and phylogenetic relationships. Plant Journal. 1997, 11 (3): 429-441. 10.1046/j.1365-313X.1997.11030429.x.PubMedView ArticleGoogle Scholar
- Poeydomenge O, Marolda M, Boudet AM, Grimapettenati J: Nucleotide sequence of cDNA-enconding mitochondrial malate-dehydrogenase from Eucalyptus. Plant Physiology. 1995, 107 (4): 1455-1456. 10.1104/pp.107.4.1455.PubMedPubMed CentralView ArticleGoogle Scholar
- Gion JM, Rech P, Grima-Pettenati J, Verhaegen D, Plomion C: Mapping candidate genes in Eucalyptus with emphasis on lignification genes. Molecular Breeding. 2000, 6 (5): 441-449. 10.1023/A:1026552515218.View ArticleGoogle Scholar
- Goicoechea M, Lacombe E, Legay S, Mihaljevic S, Rech P, Jauneau A, Lapierre C, Pollet B, Verhaegen D, Chaubet-Gigot N, et al: Eg MYB2, a new transcriptional activator from Eucalyptus xylem, regulates secondary cell wall formation and lignin biosynthesis. Plant Journal. 2005, 43 (4): 553-567. 10.1111/j.1365-313X.2005.02480.x.PubMedView ArticleGoogle Scholar
- Legay S, Lacombe E, Goicoechea M, Briere C, Seguin A, Mackay J, Grima-Pettenati J: Molecular characterization of Eg MYB1, a putative transcriptional repressor of the lignin biosynthetic pathway. Plant Science. 2007, 173 (5): 542-549. 10.1016/j.plantsci.2007.08.007.View ArticleGoogle Scholar
- Wicker T, Schlagenhauf E, Graner A, Close TJ, Keller B, Stein N: 454 sequencing put to the test using the complex genome of barley. BMC Genomics. 2006, 7-Google Scholar
- Brondani RPV, Brondani C, Tarchini R, Grattapaglia D: Development, characterization and mapping of microsatellite markers in Eucalyptus grandis and E. urophylla. Theoretical and Applied Genetics. 1998, 97 (5-6): 816-827. 10.1007/s001220050961.View ArticleGoogle Scholar
- Brondani RPV, Brondani C, Grattapaglia D: Towards a genus-wide reference linkage map for Eucalyptus based exclusively on highly informative microsatellite markers. Molecular Genetics and Genomics. 2002, 267 (3): 338-347. 10.1007/s00438-002-0665-6.PubMedView ArticleGoogle Scholar
- Brondani RPV, Williams ER, Brondani C, Grattapaglia D: A microsatellite-based consensus linkage map for species of Eucalyptus and a novel set of 230 microsatellite markers for the genus. BMC Plant Biology. 2006, 6-Google Scholar
- Kirst M, Cordeiro CM, Rezende G, Grattapaglia D: Power of microsatellite markers for fingerprinting and parentage analysis in Eucalyptus grandis breeding populations. Journal of Heredity. 2005, 96 (2): 161-166. 10.1093/jhered/esi023.PubMedView ArticleGoogle Scholar
- Steane DA, Jones RC, Vaillancourt RE: A set of chloroplast microsatellite primers for Eucalyptus (Myrtaceae). Molecular Ecology Notes. 2005, 5 (3): 538-541. 10.1111/j.1471-8286.2005.00981.x.View ArticleGoogle Scholar
- Tripathi SB, Mathish NV, Gurumurthi K: Use of genetic markers in the management of micropropagated Eucalyptus germplasm. New Forests. 2006, 31 (3): 361-372. 10.1007/s11056-005-8677-9.View ArticleGoogle Scholar
- Temnykh S, DeClerck G, Lukashova A, Lipovich L, Cartinhour S, McCouch S: Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): Frequency, length variation, transposon associations, and genetic marker potential. Genome Research. 2001, 11 (8): 1441-1452. 10.1101/gr.184001.PubMedPubMed CentralView ArticleGoogle Scholar
- Ellegren H: Microsatellites: Simple sequences with complex evolution. Nature Reviews Genetics. 2004, 5 (6): 435-445. 10.1038/nrg1348.PubMedView ArticleGoogle Scholar
- Mun JH, Kim DJ, Choi HK, Gish J, Debelle F, Mudge J, Denny R, Endre G, Saurat O, Dudez AM, et al: Distribution of microsatellites in the genome of Medicago truncatula: A resource of genetic markers that integrate genetic and physical maps. Genetics. 2006, 172 (4): 2541-2555. 10.1534/genetics.105.054791.PubMedPubMed CentralView ArticleGoogle Scholar
- Yin TM, Zhang XY, Gunter LE, Li SX, Wullschleger SD, Huang MR, Tuskan GA: Microsatellite primer resource for Populus developed from the mapped sequence scaffolds of the Nisqually-1 genome. New Phytologist. 2009, 181 (2): 498-503. 10.1111/j.1469-8137.2008.02663.x.PubMedView ArticleGoogle Scholar
- Ceresini PC, Silva C, Missio RF, Souza EC, Fischer CN, Guillherme IR, Gregorio I, da Silva EHT, Cicarelli RMB, da Silva MTA, et al: Satellyptus: Analysis and database of microsatellites from ESTs of Eucalyptus. Genetics and Molecular Biology. 2005, 28 (3): 589-600. 10.1590/S1415-47572005000400014.View ArticleGoogle Scholar
- Rabello E, de Souza AN, Saito D, Tsai SM: In silico characterization of microsatellites in Eucalyptus spp.: Abundance, length variation and transposon associations. Genetics and Molecular Biology. 2005, 28 (3): 582-588. 10.1590/S1415-47572005000400013.View ArticleGoogle Scholar
- Piquemal J, Chamayou S, Nadaud I, Beckert M, Barriere Y, Mila I, Lapierre C, Rigau J, Puigdomenech P, Jauneau A, et al: Down-regulation of caffeic acid O-methyltransf erase in maize revisited using a transgenic approach. Plant Physiology. 2002, 130 (4): 1675-1685. 10.1104/pp.012237.PubMedPubMed CentralView ArticleGoogle Scholar
- Piquemal J, Lapierre C, Myton K, O'Connell A, Schuch W, Grima-Pettenati J, Boudet AM: Down-regulation of cinnamoyl-CoA reductase induces significant changes of lignin profiles in transgenic tobacco plants. Plant Journal. 1998, 13 (1): 71-83. 10.1046/j.1365-313X.1998.00014.x.View ArticleGoogle Scholar
- Halpin C, Knight ME, Grimapettenati J, Goffner D, Boudet A, Schuch W: Purification and characterization of cinnamyl alcohol dehydrogenase from tobacco stems. Plant Physiology. 1992, 98 (1): 12-16. 10.1104/pp.98.1.12.PubMedPubMed CentralView ArticleGoogle Scholar
- Halpin C, Knight ME, Foxon GA, Campbell MM, Boudet AM, Boon JJ, Chabbert B, Tollier MT, Schuch W: Manipulation of lignin quality by down regulation of cunnamyl alcohol dehydrogenase. Plant Journal. 1994, 6 (3): 339-350. 10.1046/j.1365-313X.1994.06030339.x.View ArticleGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: The European molecular biology open software suite. Trends in Genetics. 2000, 16 (6): 276-277. 10.1016/S0168-9525(00)02024-2.PubMedView ArticleGoogle Scholar
- Jansen RK, Kaittanis C, Saski C, Lee SB, Tomkins J, Alverson AJ, Daniell H: Phylogenetic analyses of Vitis (Vitaceae) based on complete chloroplast genome sequences: effects of taxon sampling and phylogenetic methods on resolving relationships among rosids. BMC Evolutionary Biology. 2006, 6-Google Scholar
- Byrne M, Moran GF, Tibbits WN: Restriction map and material inherithance of chloroplast DNA in Eucalyptus nitens. Journal of Heredity. 1993, 84 (3): 218-220.Google Scholar
- Wakasugi T, Tsudzuki J, Ito S, Nakashima K, Tsudzuki T, Sugiura M: Loss of all NDH genes as determined by sequencing the entire chloroplast genome of the blackpine Pinus thunbergii. Proceedings of the National Academy of Sciences of the United States of America. 1994, 91 (21): 9794-9798. 10.1073/pnas.91.21.9794.PubMedPubMed CentralView ArticleGoogle Scholar
- Tuskan GA, DiFazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, et al: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313 (5793): 1596-1604. 10.1126/science.1128691.PubMedView ArticleGoogle Scholar
- Zhang HB, Zhao XP, Ding XL, Paterson AH, Wing RA: Preparation of megabase-size DNA from plant nuclei. Plant Journal. 1995, 7 (1): 175-184. 10.1046/j.1365-313X.1995.07010175.x.View ArticleGoogle Scholar
- Luo M, RA W: An improved method for plant BAC library construction. Plant Functional Genomics. Edited by: E G. 2003, Totowa, NJ.: Human Press Inc, 3-20. full_text.View ArticleGoogle Scholar
- Duminil J, Pemonge MH, Petit RJ: A set of 35 consensus primer pairs amplifying genes and introns of plant mitochondrial DNA. Molecular Ecology Notes. 2002, 2 (4): 428-430. 10.1046/j.1471-8286.2002.00263.x.View ArticleGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen ZT, et al: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437 (7057): 376-380.PubMedPubMed CentralGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen ZT, et al: Genome sequencing in microfabricated high-density picolitre reactors (vol 437, pg 376, 2005). Nature. 2006, 441 (7089): 120-120. 10.1038/nature04726.View ArticleGoogle Scholar
- Fernandes F, Pereira L, Freitas AT: CSA: An efficient algorithm to improve circular DNA multiple alignment. BMC Bioinformatics. 2009, 10: 230-10.1186/1471-2105-10-230.PubMedPubMed CentralView ArticleGoogle Scholar
- Wyman SK, Jansen RK, Boore JL: Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004, 20 (17): 3252-3255. 10.1093/bioinformatics/bth352.PubMedView ArticleGoogle Scholar
- Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J: ACT: the Artemis comparison tool. Bioinformatics. 2005, 21 (16): 3422-3423. 10.1093/bioinformatics/bti553.PubMedView ArticleGoogle Scholar