Comparative Genomics Tools Research Articles

Kerfuffle: a web tool for multi-species gene colocalization analysis

BackgroundThe evolutionary pressures that underlie the large-scale functional organization of the genome are not well understood in eukaryotes. Recent evidence suggests that functionally similar genes may colocalize (cluster) in the eukaryotic genome, suggesting the role of chromatin-level gene regulation in shaping the physical distribution of coordinated genes. However, few of the bioinformatic tools currently available allow for a systematic study of gene colocalization across several, evolutionarily distant species. Furthermore, most tools require the user to input manually curated lists of gene position information, DNA sequence or gene homology relations between species. With the growing number of sequenced genomes, there is a need to provide new comparative genomics tools that can address the analysis of multi-species gene colocalization.ResultsKerfuffle is a web tool designed to help discover, visualize, and quantify the physical organization of genomes by identifying significant gene colocalization and conservation across the assembled genomes of available species (currently up to 47, from humans to worms). Kerfuffle only requires the user to specify a list of human genes and the names of other species of interest. Without further input from the user, the software queries the e!Ensembl BioMart server to obtain positional information and discovers homology relations in all genes and species specified. Using this information, Kerfuffle performs a multi-species clustering analysis, presents downloadable lists of clustered genes, performs Monte Carlo statistical significance calculations, estimates how conserved gene clusters are across species, plots histograms and interactive graphs, allows users to save their queries, and generates a downloadable visualization of the clusters using the Circos software. These analyses may be used to further explore the functional roles of gene clusters by interrogating the enriched molecular pathways associated with each cluster.ConclusionsKerfuffle is a new, easy-to-use and publicly available tool to aid our understanding of functional genomics and comparative genomics. This software allows for flexibility and quick investigations of a user-defined set of genes, and the results may be saved online for further analysis. Kerfuffle is freely available at http://atwallab.org/kerfuffle, is implemented in JavaScript (using jQuery and jsCharts libraries) and PHP 5.2, runs on an Apache server, and stores data in flat files and an SQLite database.

Microbial Comparative Genomics: An Overview of Tools and Insights Into The Genus Corynebacterium

Next generation sequencing (NGS) made it possible to provide whole genome sequences of pathogenic and commercially significant organisms in limited time, and with minimal cost. Computational comparative genomics is necessary, given that we sequence thousands of organisms every day, but our follow-up knowledge is still very limited. Nevertheless, genomic information from a single genome is insufficient to provide insights into the life style and extended view of the gene pool of a species. Multiple genomes could enrich our understanding of the relatedness of, and variations in organisms. Consequently, comparative genomic analysis remains powerful tools for identifying the orthologous genes in species, presence and absence of specific genes, evolutionary signals, and candidate regions associated with pathogenicity. Furthermore, pangenomic strategies, together with subtractive genomics, help in highlighting the inter- and intra-species relationships, conserved core and, pan-genome for characterizing virulence factors, drug targets and vaccine candidates. In this article, we present an overview of microbial comparative genomics pre-requisites: sequencing technologies, alignment tools, annotation pipelines, databases and resources, visualization and comparative genomic tools, and strategies. Finally, we present comparative genomic and functional analysis based insights and recent findings in genus Corynebacterium.

MIPS PlantsDB: a database framework for comparative plant genome research

The rapidly increasing amount of plant genome (sequence) data enables powerful comparative analyses and integrative approaches and also requires structured and comprehensive information resources. Databases are needed for both model and crop plant organisms and both intuitive search/browse views and comparative genomics tools should communicate the data to researchers and help them interpret it. MIPS PlantsDB (http://mips.helmholtz-muenchen.de/plant/genomes.jsp) was initially described in NAR in 2007 [Spannagl,M., Noubibou,O., Haase,D., Yang,L., Gundlach,H., Hindemitt, T., Klee,K., Haberer,G., Schoof,H. and Mayer,K.F. (2007) MIPSPlantsDB–plant database resource for integrative and comparative plant genome research. Nucleic Acids Res., 35, D834–D840] and was set up from the start to provide data and information resources for individual plant species as well as a framework for integrative and comparative plant genome research. PlantsDB comprises database instances for tomato, Medicago, Arabidopsis, Brachypodium, Sorghum, maize, rice, barley and wheat. Building up on that, state-of-the-art comparative genomics tools such as CrowsNest are integrated to visualize and investigate syntenic relationships between monocot genomes. Results from novel genome analysis strategies targeting the complex and repetitive genomes of triticeae species (wheat and barley) are provided and cross-linked with model species. The MIPS Repeat Element Database (mips-REdat) and Catalog (mips-REcat) as well as tight connections to other databases, e.g. via web services, are further important components of PlantsDB.

A re-assigned American mink (Neovison vison) map optimal for genome-wide studies

Our previously published second generation genetic map for the American mink (Neovison vison) has been used and redesigned in its best for genome-wide studies with maximum of efficiency. A number of 114 selected markers, including 33 newly developed microsatellite markers from the CHORI-231 mink Bacterial Artificial Chromosome (BAC) library, have been genotyped in a two generation population composed of 1200 individuals. The outcome reassigns the position of some markers on the chromosomes and it produces a more reliable map with a convenient distance between markers. A total of 104 markers mapped to 14 linkage groups corresponding to the mink autosomes. Six markers are unlinked and four markers are allocated to the X chromosome by homology but no linkage was detected. The sex-average linkage map spans 1192 centiMorgans (cM) with an average intermarker distance of 11.4cM and 1648cM when the ends of the linkage groups and the autosomal unlinked markers are added. Sex-specific genetic linkage maps were also generated. The male sex-specific map had a total length of 1014.6cM between the linked markers and an average inter-marker interval of 9.7cM. The female map has a corresponding length of 1378.6cM and an average inter-marker interval of 13.3cM. The study is complemented with additional anchorage for most of the chromosomes of the map by BAC in situ hybridization with clones containing microsatellites strategically selected from the various parts of the genome. This map provides an improved tool for genetic mapping and comparative genomics in mink, also useful for the future assembly of the mink genome sequence when this will be taken forward.

Gene Duplication and Positive Selection Explains Unusual Physiological Roles of the Relaxin Gene in the European Rabbit

The relaxin gene family is a group of genes involved in different physiological roles, most of them related to reproduction. In vertebrates the genes in this family are located in three separate chromosomal locations, and have been called relaxin family locus (RFL) A, B, and C. Among mammals the RFLA and RFLC are the most conserved as no gene copy-number variation has been observed thus far. The RFLB locus is also conserved on most mammals other than primates, where there are several gene gains and losses. Interestingly, the relaxin gene found on the RFLB locus in the European rabbit has acquired a novel role. In addition to the classical reproductive roles, this gene is expressed in tracheobronchial epithelial cells and its expression has been linked to squamous differentiation. We reconstructed the evolutionary history of the European rabbit RFLB locus using the tools of comparative genomics and molecular evolution. We found that the European rabbit possess a RFLB locus which is unique among mammals in that there are five tandemly arranged relaxin gene copies, which contrast with the single relaxin copy gene found in most mammals. In addition we also found that the ancestral pre-duplication gene was subject to the action of positive selection, and several amino acid sites were identified under the action of natural selection including the sites B12 and B13 which are part of the receptor recognition and binding site.

Exploring the Solution Space of Sorting by Translocations

Genome sorting is a way to trace back the genomes of two species to their common ancestor. The study of genome rearrangements allows to better understand the process of evolution and is an important tool in comparative genomics. Given two genomes A and B, the goal is to find a shortest sequence of exchanges of non-empty chromosome ends that transform A into B. The length of such a shortest sequence is called the translocation distance between A and B, and the problem of finding the optimal sequence of translocations is called the sorting by translocations problem. For most problem instances, however, several minimum length sequences or optimal sequences of translocation exist, and in the absence of any additional information, no one is of greater value than the others. The problem of finding all optimal solutions for sorting by translocations is thus a natural generalization of sorting by translocations. In this paper, an algorithm is presented to find all solutions for sorting by translocations.

Tracing Lifestyle Adaptation in Prokaryotic Genomes

Lifestyle adaptation of microbes due to changes in their ecological niches or acquisition of new environments is a major driving force for genetic changes in their respective genomes. Moving into more specialized niches often results in the acquisition of new gene sets via horizontal gene transfer to utilize previously unavailable metabolites, while genetic ballast is shed by gene loss and/or gene inactivation. In some cases, larger genome rearrangements can be observed, such as the incorporation of whole genetic islands, providing a range of new phenotypic capabilities. Until recently these changes could not be comprehensively followed and identified due to the lack of complete microbial genome sequences. The advent of high-throughput DNA sequencing has dramatically changed the scientific landscape and today microbial genomes have become increasingly abundant. Currently, more than 2,900 genomes are published and more than 11,000 genome projects are listed in the Genomes Online Database‡. Although this wealth of information provides many new opportunities to assess microbial functionality, it also creates a new array of challenges when a comparison between multiple microbial genomes is required. Here, functional genome distribution (FGD) is introduced, analyzing the diversity between microbes based on their predicted ORFeome. FGD is therefore a comparative genomics approach, emphasizing the assessments of gene complements. To further facilitate the comparison between two or more genomes, degrees of amino-acid similarities between ORFeomes can be visualized in the Artemis comparison tool, graphically depicting small and large scale genome rearrangements, insertion and deletion events, and levels of similarity between individual open reading frames. FGD provides a new tool for comparative microbial genomics and the interpretation of differences in the genetic makeup of bacteria.

Anti‐inflammatory properties of intestinal Bifidobacterium strains isolated from healthy infants

Certain Bifidobacterium strains have been shown to inhibit inflammatory responses in intestinal epithelial cells. However, the precise mechanisms of these effects, including the chemical nature of the active compounds, remain to be elucidated. Here partial characterization of the anti-inflammatory properties of Bifidobacterium strains isolated from feces of healthy infants is reported. It was found that conditioned media (CM) of all strains studied are capable of attenuating tumor necrosis factor-α (TNF-α) and lipopolysaccharide- (LPS) induced inflammatory responses in the HT-29 cell line. In contrast, neither killed bifidobacterial cells, nor cell-free extracts showed such activities. Further investigations resulted in attribution of this activity to heat-stable, non-lipophilic compound(s) resistant to protease and nuclease treatments and of molecular weight less than 3 kDa. The anti-inflammatory effects were dose- and time-dependent and associated with inhibition of IκB phosphorylation and nuclear factor-κ light chain enhancer of activated B cells (NF-κB)-dependent promoter activation. The combined treatments of cells with CMs and either LPS or TNF-α, but not with CMs alone, resulted in upregulation of transforming growth factor-β1, IκBζ, and p21(CIP) mRNAs. Our data suggest certain species-specificities of the anti-inflammatory properties of bifidobacteria. This observation should prompt additional validation studies using larger set of strains and employing the tools of comparative genomics.

Meiotic maps of sockeye salmon derived from massively parallel DNA sequencing

BackgroundMeiotic maps are a key tool for comparative genomics and association mapping studies. Next-generation sequencing and genotyping by sequencing are speeding the processes of SNP discovery and the development of new genetic tools, including meiotic maps for numerous species. Currently there are limited genetic resources for sockeye salmon, Oncorhynchus nerka. We develop the first dense meiotic map for sockeye salmon using a combination of novel SNPs found in restriction site associated DNA (RAD tags) and SNPs available from existing expressed sequence tag (EST) based assays.ResultsWe discovered and genotyped putative SNPs in 3,430 RAD tags. We removed paralogous sequence variants leaving 1,672 SNPs; these were combined with 53 EST-based SNP genotypes for linkage mapping. The map contained 29 male and female linkage groups, consistent with the haploid chromosome number expected for sockeye salmon. The female map contains 1,057 loci spanning 4,896 cM, and the male map contains 1,118 loci spanning 4,220 cM. Regions of conservation with rainbow trout and synteny between the RAD based rainbow trout map and the sockeye salmon map were established.ConclusionsUsing RAD sequencing and EST-based SNP assays we successfully generated the first high density linkage map for sockeye salmon.

Using Sybil for interactive comparative genomics of microbes on the web

Motivation: Analysis of multiple genomes requires sophisticated tools that provide search, visualization, interactivity and data export. Comparative genomics datasets tend to be large and complex, making development of these tools difficult. In addition to scalability, comparative genomics tools must also provide user-friendly interfaces such that the research scientist can explore complex data with minimal technical expertise.Results: We describe a new version of the Sybil software package and its application to the important human pathogen Streptococcus pneumoniae. This new software provides a feature-rich set of comparative genomics tools for inspection of multiple genome structures, mining of orthologous gene families and identification of potential vaccine candidates.Availability: The S.pneumoniae resource is online at http://strepneumo-sybil.igs.umaryland.edu. The software, database and website are available for download as a portable virtual machine and from http://sourceforge.net/projects/sybil.Contact: driley@som.umaryland.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Phytozome: a comparative platform for green plant genomics

The number of sequenced plant genomes and associated genomic resources is growing rapidly with the advent of both an increased focus on plant genomics from funding agencies, and the application of inexpensive next generation sequencing. To interact with this increasing body of data, we have developed Phytozome (http://www.phytozome.net), a comparative hub for plant genome and gene family data and analysis. Phytozome provides a view of the evolutionary history of every plant gene at the level of sequence, gene structure, gene family and genome organization, while at the same time providing access to the sequences and functional annotations of a growing number (currently 25) of complete plant genomes, including all the land plants and selected algae sequenced at the Joint Genome Institute, as well as selected species sequenced elsewhere. Through a comprehensive plant genome database and web portal, these data and analyses are available to the broader plant science research community, providing powerful comparative genomics tools that help to link model systems with other plants of economic and ecological importance.

Characterization of the genome of bald cypress

BackgroundBald cypress (Taxodium distichum var. distichum) is a coniferous tree of tremendous ecological and economic importance. It is a member of the family Cupressaceae which also includes cypresses, redwoods, sequoias, thujas, and junipers. While the bald cypress genome is more than three times the size of the human genome, its 1C DNA content is amongst the smallest of any conifer. To learn more about the genome of bald cypress and gain insight into the evolution of Cupressaceae genomes, we performed a Cot analysis and used Cot filtration to study Taxodium DNA. Additionally, we constructed a 6.7 genome-equivalent BAC library that we screened with known Taxodium genes and select repeats.ResultsThe bald cypress genome is composed of 90% repetitive DNA with most sequences being found in low to mid copy numbers. The most abundant repeats are found in fewer than 25,000 copies per genome. Approximately 7.4% of the genome is single/low-copy DNA (i.e., sequences found in 1 to 5 copies). Sequencing of highly repetitive Cot clones indicates that most Taxodium repeats are highly diverged from previously characterized plant repeat sequences. The bald cypress BAC library consists of 606,336 clones (average insert size of 113 kb) and collectively provides 6.7-fold genome equivalent coverage of the bald cypress genome. Macroarray screening with known genes produced, on average, about 1.5 positive clones per probe per genome-equivalent. Library screening with Cot-1 DNA revealed that approximately 83% of BAC clones contain repetitive sequences iterated 103 to 104 times per genome.ConclusionsThe BAC library for bald cypress is the first to be generated for a conifer species outside of the family Pinaceae. The Taxodium BAC library was shown to be useful in gene isolation and genome characterization and should be an important tool in gymnosperm comparative genomics, physical mapping, genome sequencing, and gene/polymorphism discovery. The single/low-copy (SL) component of bald cypress is 4.6 times the size of the Arabidopsis genome. As suggested for other gymnosperms, the large amount of SL DNA in Taxodium is likely the result of divergence among ancient repeat copies and gene/pseudogene duplication.

Conserved and non-conserved enhancers direct tissue specific transcription in ancient germ layer specific developmental control genes

BackgroundIdentifying DNA sequences (enhancers) that direct the precise spatial and temporal expression of developmental control genes remains a significant challenge in the annotation of vertebrate genomes. Locating these sequences, which in many cases lie at a great distance from the transcription start site, has been a major obstacle in deciphering gene regulation. Coupling of comparative genomics with functional validation to locate such regulatory elements has been a successful method in locating many such regulatory elements. But most of these studies looked either at a single gene only or the whole genome without focusing on any particular process. The pressing need is to integrate the tools of comparative genomics with knowledge of developmental biology to validate enhancers for developmental transcription factors in greater detailResultsOur results show that near four different genes (nkx3.2, pax9, otx1b and foxa2) in zebrafish, only 20-30% of highly conserved DNA sequences can act as developmental enhancers irrespective of the tissue the gene expresses in. We find that some genes also have multiple conserved enhancers expressing in the same tissue at the same or different time points in development. We also located non-conserved enhancers for two of the genes (pax9 and otx1b). Our modified Bacterial artificial chromosome (BACs) studies for these 4 genes revealed that many of these enhancers work in a synergistic fashion, which cannot be captured by individual DNA constructs and are not conserved at the sequence level. Our detailed biochemical and transgenic analysis revealed Foxa1 binds to the otx1b non-conserved enhancer to direct its activity in forebrain and otic vesicle of zebrafish at 24 hpf.ConclusionOur results clearly indicate that high level of functional conservation of genes is not necessarily associated with sequence conservation of its regulatory elements. Moreover certain non conserved DNA elements might have role in gene regulation. The need is to bring together multiple approaches to bear upon individual genes to decipher all its regulatory elements.

Complete Genome Sequence of Human Adenovirus Prototype 17

As one of the first five human adenoviruses (HAdVs) to be sequenced, type 17 was important as a reference tool for comparative genomics of recently isolated HAdV pathogens in species D. HAdV-D17 was the first species D adenovirus to be sequenced and was deposited in GenBank in 1999. These genome data were not of high quality, and a redetermination of the same stock virus provides corrected data; among the differences are a length of 35,139 bp versus 35,100 bp in the original, and 160 mismatches to the original genome were found. Annotation of the coding sequences reveals 39 as opposed to 8, a finding which is important for phylogenomic studies.

Possible chemotaxis in Ruminococcus albus: comparative genomics

Several comparative genomic tools were used to determine the existence of chemiotaxis in ruminal bacteria. Comparative analysis of microbial genomes (database comprehensive microbial resource (CMR)) was used to search for a specific chemiotaxis gene. Then, short sequences of Ruminococcus albus were searched in the database of the National Center for Biotechnology Information, and compared them in the and Concise Microbial Protein. To predict the number of trans-membrane helices of the putative protein, the methyl-accepting chemotaxis proteins (MCP) program was used (TMpred and TopPred). For comparisons of putative protein structures in the MCP and chemotaxis protein, Simple Modular Architecture Research Tool and protein families database were used. Short sequences of the genome of R. albus revealed the presence of chemotaxis genes that could encode a chemoreceptor (MCP) and chemotaxis proteins. Groups of R. albus chemotaxis genes may be responsible for a diverse set of signalling functions, such as the formation of biofilm, adhesion and gene regulation.

Identification of Xanthomonas fragariae, Xanthomonas axonopodis pv. phaseoli, and Xanthomonas fuscans subsp. fuscans with novel markers and using a dot blot platform coupled with automatic data analysis.

Phytosanitary regulations and the provision of plant health certificates still rely mainly on long and laborious culture-based methods of diagnosis, which are frequently inconclusive. DNA-based methods of detection can circumvent many of the limitations of currently used screening methods, allowing a fast and accurate monitoring of samples. The genus Xanthomonas includes 13 phytopathogenic quarantine organisms for which improved methods of diagnosis are needed. In this work, we propose 21 new Xanthomonas-specific molecular markers, within loci coding for Xanthomonas-specific protein domains, useful for DNA-based methods of identification of xanthomonads. The specificity of these markers was assessed by a dot blot hybridization array using 23 non-Xanthomonas species, mostly soil dwelling and/or phytopathogens for the same host plants. In addition, the validation of these markers on 15 Xanthomonas spp. suggested species-specific hybridization patterns, which allowed discrimination among the different Xanthomonas species. Having in mind that DNA-based methods of diagnosis are particularly hampered for unsequenced species, namely, Xanthomonas fragariae, Xanthomonas axonopodis pv. phaseoli, and Xanthomonas fuscans subsp. fuscans, for which comparative genomics tools to search for DNA signatures are not yet applicable, emphasis was given to the selection of informative markers able to identify X. fragariae, X. axonopodis pv. phaseoli, and X. fuscans subsp. fuscans strains. In order to avoid inconsistencies due to operator-dependent interpretation of dot blot data, an image-processing algorithm was developed to analyze automatically the dot blot patterns. Ultimately, the proposed markers and the dot blot platform, coupled with automatic data analyses, have the potential to foster a thorough monitoring of phytopathogenic xanthomonads.

GLAMM: Genome-Linked Application for Metabolic Maps

The Genome-Linked Application for Metabolic Maps (GLAMM) is a unified web interface for visualizing metabolic networks, reconstructing metabolic networks from annotated genome data, visualizing experimental data in the context of metabolic networks and investigating the construction of novel, transgenic pathways. This simple, user-friendly interface is tightly integrated with the comparative genomics tools of MicrobesOnline [Dehal et al. (2010) Nucleic Acids Research, 38, D396–D400]. GLAMM is available for free to the scientific community at glamm.lbl.gov.

Molecular Characterizations of a Novel Putative DNA-Binding Protein LvDBP23 in Marine Shrimp L. vannamei Tissues and Molting Stages

Background Litopenaeus Vannamei, well known as pacific white shrimp, is the most popular shrimp in the world shrimp market. Identification and characterization of shrimp muscle regulatory genes are not only important for shrimp genetic improvement, but also facilitate comparative genomic tools for understanding of muscle development and regeneration.Methodology/Principal FindingsA novel mRNA encoding for a putative DNA-binding protein LvDBP23 was identified from Litopenaeus vannamei abdominal muscle cDNA library. The LvDBP23 cDNA contains 639 nucleotides of protein-coding sequence with deduced 212 amino acids of predicted molecular mass 23.32 kDa with glycine-rich domain at amino acid position 94–130. The mRNA sequence is successfully used for producing LvDBP23 recombinant protein in sf9 insect cell expression system. The expression of LvDBP23 mRNA is presented in abdominal muscle and swimming leg muscle, as well as other tissues including intestine, lymphoid and gill. The mRNA expression has the highest level in abdominal muscle in all tested tissues. LVDBP23 transcript during the molt cycle is highly expressed in the intermolt stage. In vitro nucleic acid-binding assays reveal that LvDBP23 protein can bind to both ssDNA and dsDNA, indicating its possible role of regulation of gene transcription.Conclusions/SignificanceWe are the first to report a DNA-binding protein identified from the abdominal muscle tissue of marine shrimp L. Vannamei. Its high-level specific expression during the intermot stage suggests its role in the regulation of muscle buildup during the growth phase of shrimp molt cycle.

Using Genomic Sequencing for Classical Genetics in E. coli K12

We here develop computational methods to facilitate use of 454 whole genome shotgun sequencing to identify mutations in Escherichia coli K12. We had Roche sequence eight related strains derived as spontaneous mutants in a background without a whole genome sequence. They provided difference tables based on assembling each genome to reference strain E. coli MG1655 (NC_000913). Due to the evolutionary distance to MG1655, these contained a large number of both false negatives and positives. By manual analysis of the dataset, we detected all the known mutations (24 at nine locations) and identified and genetically confirmed new mutations necessary and sufficient for the phenotypes we had selected in four strains. We then had Roche assemble contigs de novo, which we further assembled to full-length pseudomolecules based on synteny with MG1655. This hybrid method facilitated detection of insertion mutations and allowed annotation from MG1655. After removing one genome with less than the optimal 20- to 30-fold sequence coverage, we identified 544 putative polymorphisms that included all of the known and selected mutations apart from insertions. Finally, we detected seven new mutations in a total of only 41 candidates by comparing single genomes to composite data for the remaining six and using a ranking system to penalize homopolymer sequencing and misassembly errors. An additional benefit of the analysis is a table of differences between MG1655 and a physiologically robust E. coli wild-type strain NCM3722. Both projects were greatly facilitated by use of comparative genomics tools in the CoGe software package (http://genomevolution.org/).

A second generation genetic map of the bumblebee Bombus terrestris (Linnaeus, 1758) reveals slow genome and chromosome evolution in the Apidae

BackgroundThe bumblebee Bombus terrestris is an ecologically and economically important pollinator and has become an important biological model system. To study fundamental evolutionary questions at the genomic level, a high resolution genetic linkage map is an essential tool for analyses ranging from quantitative trait loci (QTL) mapping to genome assembly and comparative genomics. We here present a saturated linkage map and match it with the Apis mellifera genome using homologous markers. This genome-wide comparison allows insights into structural conservations and rearrangements and thus the evolution on a chromosomal level.ResultsThe high density linkage map covers ~ 93% of the B. terrestris genome on 18 linkage groups (LGs) and has a length of 2'047 cM with an average marker distance of 4.02 cM. Based on a genome size of ~ 430 Mb, the recombination rate estimate is 4.76 cM/Mb. Sequence homologies of 242 homologous markers allowed to match 15 B. terrestris with A. mellifera LGs, five of them as composites. Comparing marker orders between both genomes we detect over 14% of the genome to be organized in synteny and 21% in rearranged blocks on the same homologous LG.ConclusionsThis study demonstrates that, despite the very high recombination rates of both A. mellifera and B. terrestris and a long divergence time of about 100 million years, the genomes' genetic architecture is highly conserved. This reflects a slow genome evolution in these bees. We show that data on genome organization and conserved molecular markers can be used as a powerful tool for comparative genomics and evolutionary studies, opening up new avenues of research in the Apidae.