Structural Genomic Analysis Research Articles

G-SAIP: Graphical Sequence Alignment Through Parallel Programming in the Post-Genomic Era.

A common task in bioinformatics is to compare DNA sequences to identify similarities between organisms at the sequence level. An approach to such comparison is the dot-plots, a 2-dimensional graphical representation to analyze DNA or protein alignments. Dot-plots alignment software existed before the sequencing revolution, and now there is an ongoing limitation when dealing with large-size sequences, resulting in very long execution times. High-Performance Computing (HPC) techniques have been successfully used in many applications to reduce computing times, but so far, very few applications for graphical sequence alignment using HPC have been reported. Here, we present G-SAIP (Graphical Sequence Alignment in Parallel), a software capable of spawning multiple distributed processes on CPUs, over a supercomputing infrastructure to speed up the execution time for dot-plot generation up to 1.68× compared with other current fastest tools, improve the efficiency for comparative structural genomic analysis, phylogenetics because the benefits of pairwise alignments for comparison between genomes, repetitive structure identification, and assembly quality checking.

Structural Genomic Analysis of SARS-CoV-2 and Other Coronaviruses.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is the causative agent of the coronavirus disease 2019 (COVID-19) pandemic. In this study, we conducted a comparative analysis of the structural genes of SARS-CoV-2 and other CoVs. We found that the sequence of the E gene was the most evolutionarily conserved across 200 SARS-CoV-2 isolates. The E gene and M gene sequences of SARS-CoV-2 and NC014470 CoV were closely related and fell within the same branch of a phylogenetic tree. The absolute diversity of E gene and M gene sequences of SARS-CoV-2 isolates was similar to that of common CoVs (C-CoVs) infecting other organisms. The absolute diversity of the M gene sequence of the KJ481931 CoV that can infect humans was similar to that of SARS-CoV-2 and C-CoVs infecting other organisms. The M gene sequence of KJ481931 CoV (infecting humans), SARS-CoV-2 and NC014470 CoV (infecting other organisms) were closely related, falling within the same branch of a phylogenetic tree. Patterns of variation and evolutionary characteristics of the N gene and S gene were very similar. These data may be of value for understanding the origins and intermediate hosts of SARS-CoV-2.

Structural Bioinformatic and Genomic Analyses of RORg Interacting with Cognate Response Elements

Structural Bioinformatic and Genomic Analyses of RORg Interacting with Cognate Response Elements

MITOCHONDRIAL GENOME STRUCTURE OF MACROSIPHUM ALBIFRONS ESSIG, 1911 — A DANGEROUS PEST OF LUPINE

This work provides information on the assembly, annotation, and results of the structural analysis of the Macrosiphum albifrons mitochondrial genome. The mitochondrial genome of M. albifrons has an ancestral ability and a deviation towards the content of adenine bases in their saturation region of repeats. The presence of a complete set of tRNAs was shown, in 14 of which the paired bases are aberrant.

Structural genome analysis in cultivated potato taxa

Key messageTwelve potato accessions were selected to represent two principal views on potato taxonomy. The genomes were sequenced and analyzed for structural variation (copy number variation) against three published potato genomes.The common potato (Solanum tuberosum L.) is an important staple crop with a highly heterozygous and complex tetraploid genome. The other taxa of cultivated potato contain varying ploidy levels (2X–5X), and structural variations are common in the genomes of these species, likely contributing to the diversification or agronomic traits during domestication. Increased understanding of the genomes and genomic variation will aid in the exploration of novel agronomic traits. Thus, sequencing data from twelve potato landraces, representing the four ploidy levels, were used to identify structural genomic variation compared to the two currently available reference genomes, a double monoploid potato genome and a diploid inbred clone of S. chacoense. The results of a copy number variation analysis showed that in the majority of the genomes, while the number of deletions is greater than the number of duplications, the number of duplicated genes is greater than the number of deleted ones. Specific regions in the twelve potato genomes have a high density of CNV events. Further, the auxin-induced SAUR genes (involved in abiotic stress), disease resistance genes and the 2-oxoglutarate/Fe(II)-dependent oxygenase superfamily proteins, among others, had increased copy numbers in these sequenced genomes relative to the references.

Structural and Comparative Analysis of the Complete Chloroplast Genome of Pyrus hopeiensis-"Wild Plants with a Tiny Population"-and Three Other Pyrus Species.

Pyrus hopeiensis is a valuable wild resource of Pyrus in the Rosaceae. Due to its limited distribution and population decline, it has been listed as one of the “wild plants with a tiny population” in China. To date, few studies have been conducted on P. hopeiensis. This paper offers a systematic review of P. hopeiensis, providing a basis for the conservation and restoration of P. hopeiensis resources. In this study, the chloroplast genomes of two different genotypes of P. hopeiensis, P. ussuriensis Maxin. cv. Jingbaili, P. communis L. cv. Early Red Comice, and P. betulifolia were sequenced, compared and analyzed. The two P. hopeiensis genotypes showed a typical tetrad chloroplast genome, including a pair of inverted repeats encoding the same but opposite direction sequences, a large single copy (LSC) region, and a small single copy (SSC) region. The length of the chloroplast genome of P. hopeiensis HB-1 was 159,935 bp, 46 bp longer than that of the chloroplast genome of P. hopeiensis HB-2. The lengths of the SSC and IR regions of the two Pyrus genotypes were identical, with the only difference present in the LSC region. The GC content was only 0.02% higher in P. hopeiensis HB-1. The structure and size of the chloroplast genome, the gene species, gene number, and GC content of P. hopeiensis were similar to those of the other three Pyrus species. The IR boundary of the two genotypes of P. hopeiensis showed a similar degree of expansion. To determine the evolutionary history of P. hopeiensis within the genus Pyrus and the Rosaceae, 57 common protein-coding genes from 36 Rosaceae species were analyzed. The phylogenetic tree showed a close relationship between the genera Pyrus and Malus, and the relationship between P. hopeiensis HB-1 and P. hopeiensis HB-2 was the closest.

Evolutionary acquired robustness and vulnerability in cancer genome: negligible negative selection in carcinogenesis

Cancer is a disease involving the accumulation of somatic genetic alterations accompanied by clonal evolution (1,2). A wide range of oncogenes and tumor suppressor genes are known to be mutated, deleted, rearranged, and/or amplified in most cancers. These genes have been identified by functional analyses (transformation of NIH3T3), genetic linkage analyses, and structural genomic analyses such as cytogenetics.

Crystal Structure of Shigella flexneri SF173 Reveals a Dimeric Helical Bundle Conformation

We report the crystal structure and bioinformatic analysis of SF173, a functionally uncharacterized protein from the human enteropathogenic bacteria Shigella flexneri. The structure shows a tightly interlinked dimer formed by adimeric core comprising α2 and α3 helices from both subunits and swapping the N-terminal α1 helix of each monomer. Structural inspection and genomic analysis results suggest that the SF173 might play its putative function by binding to SF172, the partially overlapped upstream product in the operon. As YaeO (an SF172 orthologue) has been identified to be an inhibitor of the bacterial transcription terminator Rho protein, SF173 is suggested to be involved in the regulation of Rho-dependent transcription termination, by inhibiting the Rho protein binding to SF172/YaeO.

Adult high-grade B-cell lymphoma with Burkitt lymphoma signature: genomic features and potential therapeutic targets

AbstractThe adult high-grade B-cell lymphomas sharing molecular features with Burkitt lymphoma (BL) are highly aggressive lymphomas with poor clinical outcome. High-resolution structural and functional genomic analysis of adult Burkitt lymphoma (BL) and high-grade B-cell lymphoma with BL gene signature (adult-molecularly defined BL [mBL]) revealed the MYC-ARF-p53 axis as the primary deregulated pathway. Adult-mBL had either unique or more frequent genomic aberrations (del13q14, del17p, gain8q24, and gain18q21) compared with pediatric-mBL, but shared commonly mutated genes. Mutations in genes promoting the tonic B-cell receptor (BCR)→PI3K pathway (TCF3 and ID3) did not differ by age, whereas effectors of chronic BCR→NF-κB signaling were associated with adult-mBL. A subset of adult-mBL had BCL2 translocation and mutation and elevated BCL2 mRNA and protein expression, but had a mutation profile similar to mBL. These double-hit lymphomas may have arisen from a tumor precursor that acquired both BCL2 and MYC translocations and/or KMT2D (MLL2) mutation. Gain/amplification of MIR17HG and its paralogue loci was observed in 50% of adult-mBL. In vitro studies suggested miR-17∼92's role in constitutive activation of BCR signaling and sensitivity to ibrutinib. Overall integrative analysis identified an interrelated gene network affected by copy number and mutation, leading to disruption of the p53 pathway and the BCR→PI3K or NF-κB activation, which can be further exploited in vivo by small-molecule inhibitors for effective therapy in adult-mBL.

Forward genetic screen of human transposase genomic rearrangements.

BackgroundNumerous human genes encode potentially active DNA transposases or recombinases, but our understanding of their functions remains limited due to shortage of methods to profile their activities on endogenous genomic substrates.ResultsTo enable functional analysis of human transposase-derived genes, we combined forward chemical genetic hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1) screening with massively parallel paired-end DNA sequencing and structural variant genome assembly and analysis. Here, we report the HPRT1 mutational spectrum induced by the human transposase PGBD5, including PGBD5-specific signal sequences (PSS) that serve as potential genomic rearrangement substrates.ConclusionsThe discovered PSS motifs and high-throughput forward chemical genomic screening approach should prove useful for the elucidation of endogenous genome remodeling activities of PGBD5 and other domesticated human DNA transposases and recombinases.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2877-x) contains supplementary material, which is available to authorized users.

Structural organization of fatty acid desaturase loci in linseed lines with contrasting linolenic acid contents.

Flax (Linum usitatissimum L.), the richest crop source of omega-3 fatty acids (FAs), is a diploid plant with an estimated genome size of ~370Mb and is well suited for studying genomic organization of agronomically important traits. In this study, 12 bacterial artificial chromosome clones harbouring the six FA desaturase loci sad1, sad2, fad2a, fad2b, fad3a and fad3b from the conventional variety CDC Bethune and the high linolenic acid line M5791 were sequenced, analysed and compared to determine the structural organization of these loci and to gain insights into the genetic mechanisms underlying FA composition in flax. With one gene every 3.2-4.6kb, the desaturase loci have a higher gene density than the genome's average of one gene per 7.8-8.2kb. The gene order and orientation across the two genotypes were generally conserved with the exception of the sad1 locus that was predicted to have additional genes in CDC Bethune. High sequence conservation in both genic and intergenic regions of the sad and fad2b loci contrasted with the significant level of variation of the fad2a and fad3 loci, with SNPs being the most frequently observed mutation type. The fad2a locus had 297 SNPs and 36 indels over ~95kb contrasting with the fad2b locus that had a mere seven SNPs and four indels in ~110kb. Annotation of the gene-rich loci revealed other genes of known role in lipid or carbohydrate metabolic/catabolic pathways. The organization of the fad2b locus was particularly complex with seven copies of the fad2b gene in both genotypes. The presence of Gypsy, Copia, MITE, Mutator, hAT and other novel repeat elements at the desaturase loci was similar to that of the whole genome. This structural genomic analysis provided some insights into the genomic organization and composition of the main desaturase loci of linseed and of their complex evolution through both tandem and whole genome duplications.

Genomic structural variants are linked with intellectual disability.

Mutations in more than 500 genes have been associated with intellectual disability (ID) and related disorders of cognitive function, such as autism and schizophrenia. Here we aimed to unravel the molecular epidemiology of non-specific ID in a genetic isolate using a combination of population and molecular genetic approaches. A large multigenerational pedigree was ascertained within a Dagestan Genetic Heritage research program in a genetic isolate of indigenous ethnics. Clinical characteristics of the affected members were based on combining diagnoses from regional psychiatric hospitals with our own clinical assessment, using a Russian translation of the structured psychiatric interviews, the Diagnostic Interview for Genetic Studies and the Family Interview for Genetic Studies, based on DSM-IV criteria. Weber/CHLC 9.0 STRs set was used for multipoint parametric linkage analyses (Simwalk2.91). Next, we checked CNVs and LOH (based on Affymetrix SNP 5.0 data) in the linked with ID genomic regions with the aim to identify candidate genes associated with mutations in linked regions. The number of statistically significant (p ≤ 0.05) suggestive linkage peaks with 1.3 < LOD < 3.0 we detected in a total of 10 genomic regions: 1q41, 2p25.3-p24.2, 3p13-p12.1, 4q13.3, 10p11, 11q23, 12q24.22-q24.31, 17q24.2-q25.1, 21q22.13 and 22q12.3-q13.1. Three significant linkage signals with LOD >3 were obtained at 2p25.3-p24.2 under the dominant model, with a peak at 21 cM flanked by loci D2S2976 and D2S2952; at 12q24.22-q24.31 under the recessive model, with a peak at -120 cM flanked by marker D12S2070 and D12S395 and at 22q12.3 under the dominant model, with a peak at 32 cM flanked by marker D22S683 and D22S445. After a set of genes had been designated as possible candidates in these specific chromosomal regions,we conducted an exploratory search for LOH and CNV based on microarray data to detect structural genomic variants within five ID-linked regions with LOD scores between 2.0 and 3.9. In these selected regions we obtained 173 ROH segments and 98 CN segments. Further analysis of region 2p25.3-p24.2 revealed deletions within genes encoding MYTL, SNTG2 and TPO among five of 21 affected cases at 2p25.3-p24.2. In the ID-linked region at 12q24.22-12q24.31 19 out of 21 ID cases carried segmental CNV and 20 of 21 them displayed ROH segments with mean size lengths for ID cases 2512 kb (500-6,472 kb) and for healthy control 682 kb (531-986 kb), including the genes MED13L, HRK, FBXW8, TESC, CDK2AP1 and SBNO1. Seven of 21 affected pedigree members displayed segmental deletions at 22q12.3 that includes the gene LARGE. Eight affected pedigree members carried ROH segments and 6 CN segments at 10p11.23-p11.21 containing the genes ZEB1, c10orf68 and EPC1. Our linkage and structural genomic variation analyses in a remote highland genetic isolate with aggregation of ID demonstrated that even highly isolated single kindred ID has oligo/polygenic pathogenesis. The results obtained implicate 10 genomic regions linked with ID that contain some of previously reported candidate genes, including HRK, FBXW8, TESC, CDK2AP1 and SBNO1 at 12q24 that were shown in recent studies as associated with brain measures derived from MRI scans.

Mapping contigs using CONTIGuator.

Obtaining bacterial genomic sequences has become a routine task in today's biology. The emergence of the comparative genomics approach has led to an increasing number of bacterial species having more than one strain sequenced, thus facilitating the annotation process. On the other hand, many genomic sequences are now left in the "draft" status, as a series of contigs, mainly for the labor-intensive finishing task. As a result, many genomic analyses are incomplete (e.g., in their annotation) or impossible to be performed (e.g., structural genomics analysis). Many approaches have been recently developed to facilitate the finishing process or at least to produce higher quality scaffolds; taking advantage of the comparative genomics paradigm, closely related genomes are used to align the contigs and determine their relative orientation, thus facilitating the finishing process, but also producing higher quality scaffolds. In this chapter we present the use of the CONTIGuator algorithm, which aligns the contigs from a draft genome to a closely related closed genome and resolves their relative orientation based on this alignment, producing a scaffold and a series of PCR primer pairs for the finishing process. The CONTIGuator algorithm is also capable of handling multipartite genomes (i.e., genomes having chromosomes and other plasmids), univocally mapping contigs to the most similar replicon. The program also produces a series of contig maps that allow to perform structural genomics analysis on the draft genome. The functionalities of the web interface, as well as the command line version, are presented.

Analysis of the small RNA spf in the plant pathogen Pseudomonas syringae pv. tomato strain DC3000

Bacteria contain small non-coding RNAs (ncRNAs) that are typically responsible for altering transcription, translation or mRNA stability. ncRNAs are important because they often regulate virulence factors and susceptibility to various stresses. Here, the regulation of a recently described ncRNA of Pseudomonas syringae DC3000, spot 42 (now referred to as spf), was investigated. A putative RpoE binding site was identified upstream of spf in strain DC3000. RpoE is shown to regulate the expression of spf. Also, deletion of spf results in increased sensitivity to hydrogen peroxide compared with the wild-type strain, suggesting that spf plays a role in susceptibility to oxidative stress. Furthermore, expression of alg8 is shown to be influenced by spf, suggesting that this ncRNA plays a role in alginate biosynthesis. Structural and comparative genomic analyses show this ncRNA is well conserved among the pseudomonads. The findings provide new information on the regulation and role of this ncRNA in P. syringae.

Structural genomics analysis of uncharacterized protein families overrepresented in human gut bacteria identifies a novel glycoside hydrolase.

BackgroundBacteroides spp. form a significant part of our gut microbiome and are well known for optimized metabolism of diverse polysaccharides. Initial analysis of the archetypal Bacteroides thetaiotaomicron genome identified 172 glycosyl hydrolases and a large number of uncharacterized proteins associated with polysaccharide metabolism.ResultsBT_1012 from Bacteroides thetaiotaomicron VPI-5482 is a protein of unknown function and a member of a large protein family consisting entirely of uncharacterized proteins. Initial sequence analysis predicted that this protein has two domains, one on the N- and one on the C-terminal. A PSI-BLAST search found over 150 full length and over 90 half size homologs consisting only of the N-terminal domain. The experimentally determined three-dimensional structure of the BT_1012 protein confirms its two-domain architecture and structural analysis of both domains suggests their specific functions. The N-terminal domain is a putative catalytic domain with significant similarity to known glycoside hydrolases, the C-terminal domain has a beta-sandwich fold typically found in C-terminal domains of other glycosyl hydrolases, however these domains are typically involved in substrate binding. We describe the structure of the BT_1012 protein and discuss its sequence-structure relationship and their possible functional implications.ConclusionsStructural and sequence analyses of the BT_1012 protein identifies it as a glycosyl hydrolase, expanding an already impressive catalog of enzymes involved in polysaccharide metabolism in Bacteroides spp. Based on this we have renamed the Pfam families representing the two domains found in the BT_1012 protein, PF13204 and PF12904, as putative glycoside hydrolase and glycoside hydrolase-associated C-terminal domain respectively.

Integrated genetic and epigenetic analysis of childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the commonest childhood malignancy and is characterized by recurring structural genetic alterations. Previous studies of DNA methylation suggest epigenetic alterations may also be important, but an integrated genome-wide analysis of genetic and epigenetic alterations in ALL has not been performed. We analyzed 137 B-lineage and 30 T-lineage childhood ALL cases using microarray analysis of DNA copy number alterations and gene expression, and genome-wide cytosine methylation profiling using the HpaII tiny fragment enrichment by ligation-mediated PCR (HELP) assay. We found that the different genetic subtypes of ALL are characterized by distinct DNA methylation signatures that exhibit significant correlation with gene expression profiles. We also identified an epigenetic signature common to all cases, with correlation to gene expression in 65% of these genes, suggesting that a core set of epigenetically deregulated genes is central to the initiation or maintenance of lymphoid transformation. Finally, we identified aberrant methylation in multiple genes also targeted by recurring DNA copy number alterations in ALL, suggesting that these genes are inactivated far more frequently than suggested by structural genomic analyses alone. Together, these results demonstrate subtype- and disease-specific alterations in cytosine methylation in ALL that influence transcriptional activity, and are likely to exert a key role in leukemogenesis.

Structural analysis of the genome of breast cancer cell line ZR-75-30 identifies twelve expressed fusion genes

BackgroundIt has recently emerged that common epithelial cancers such as breast cancers have fusion genes like those in leukaemias. In a representative breast cancer cell line, ZR-75-30, we searched for fusion genes, by analysing genome rearrangements.ResultsWe first analysed rearrangements of the ZR-75-30 genome, to around 10kb resolution, by molecular cytogenetic approaches, combining array painting and array CGH. We then compared this map with genomic junctions determined by paired-end sequencing. Most of the breakpoints found by array painting and array CGH were identified in the paired end sequencing—55% of the unamplified breakpoints and 97% of the amplified breakpoints (as these are represented by more sequence reads). From this analysis we identified 9 expressed fusion genes: APPBP2-PHF20L1, BCAS3-HOXB9, COL14A1-SKAP1, TAOK1-PCGF2, TIAM1-NRIP1, TIMM23-ARHGAP32, TRPS1-LASP1, USP32-CCDC49 and ZMYM4-OPRD1. We also determined the genomic junctions of a further three expressed fusion genes that had been described by others, BCAS3-ERBB2, DDX5-DEPDC6/DEPTOR and PLEC1-ENPP2. Of this total of 12 expressed fusion genes, 9 were in the coamplification. Due to the sensitivity of the technologies used, we estimate these 12 fusion genes to be around two-thirds of the true total. Many of the fusions seem likely to be driver mutations. For example, PHF20L1, BCAS3, TAOK1, PCGF2, and TRPS1 are fused in other breast cancers. HOXB9 and PHF20L1 are members of gene families that are fused in other neoplasms. Several of the other genes are relevant to cancer—in addition to ERBB2, SKAP1 is an adaptor for Src, DEPTOR regulates the mTOR pathway and NRIP1 is an estrogen-receptor coregulator.ConclusionsThis is the first structural analysis of a breast cancer genome that combines classical molecular cytogenetic approaches with sequencing. Paired-end sequencing was able to detect almost all breakpoints, where there was adequate read depth. It supports the view that gene breakage and gene fusion are important classes of mutation in breast cancer, with a typical breast cancer expressing many fusion genes.

CONTIGuator: a bacterial genomes finishing tool for structural insights on draft genomes

Recent developments in sequencing technologies have given the opportunity to sequence many bacterial genomes with limited cost and labor, compared to previous techniques. However, a limiting step of genome sequencing is the finishing process, needed to infer the relative position of each contig and close sequencing gaps. An additional degree of complexity is given by bacterial species harboring more than one replicon, which are not contemplated by the currently available programs. The availability of a large number of bacterial genomes allows geneticists to use complete genomes (possibly from the same species) as templates for contigs mapping.Here we present CONTIGuator, a software tool for contigs mapping over a reference genome which allows the visualization of a map of contigs, underlining loss and/or gain of genetic elements and permitting to finish multipartite genomes. The functionality of CONTIGuator was tested using four genomes, demonstrating its improved performances compared to currently available programs.Our approach appears efficient, with a clear visualization, allowing the user to perform comparative structural genomics analysis on draft genomes. CONTIGuator is a Python script for Linux environments and can be used on normal desktop machines and can be downloaded from http://contiguator.sourceforge.net.

Reconstructing the History of Large-Scale Genomic Changes: Biological Questions and Computational Challenges

In addition to point mutations, larger-scale structural changes (including rearrangements, duplications, insertions, and deletions) are also prevalent between different mammalian genomes. Capturing these large-scale changes is critical to unraveling the history of mammalian evolution in order to better understand the human genome. It also has profound biomedical significance, because many human diseases are associated with structural genomic aberrations. The increasing number of mammalian genomes being sequenced as well as the recent advancement in DNA sequencing technologies are allowing us to identify these structural genomic changes with vastly greater accuracy. However, there are a considerable number of computational challenges related to these problems. In this article, we introduce the ancestral genome reconstruction problem, which enables us to explain the large-scale genomic changes between species in an evolutionary context. The application of these methods to within-species structural variation and disease genome analysis is also discussed. The target audience of this article is advanced undergraduate students in biology.

Genomic insight into the common carp (Cyprinus carpio) genome by sequencing analysis of BAC-end sequences

BackgroundCommon carp is one of the most important aquaculture teleost fish in the world. Common carp and other closely related Cyprinidae species provide over 30% aquaculture production in the world. However, common carp genomic resources are still relatively underdeveloped. BAC end sequences (BES) are important resources for genome research on BAC-anchored genetic marker development, linkage map and physical map integration, and whole genome sequence assembling and scaffolding.ResultTo develop such valuable resources in common carp (Cyprinus carpio), a total of 40,224 BAC clones were sequenced on both ends, generating 65,720 clean BES with an average read length of 647 bp after sequence processing, representing 42,522,168 bp or 2.5% of common carp genome. The first survey of common carp genome was conducted with various bioinformatics tools. The common carp genome contains over 17.3% of repetitive elements with GC content of 36.8% and 518 transposon ORFs. To identify and develop BAC-anchored microsatellite markers, a total of 13,581 microsatellites were detected from 10,355 BES. The coding region of 7,127 genes were recognized from 9,443 BES on 7,453 BACs, with 1,990 BACs have genes on both ends. To evaluate the similarity to the genome of closely related zebrafish, BES of common carp were aligned against zebrafish genome. A total of 39,335 BES of common carp have conserved homologs on zebrafish genome which demonstrated the high similarity between zebrafish and common carp genomes, indicating the feasibility of comparative mapping between zebrafish and common carp once we have physical map of common carp.ConclusionBAC end sequences are great resources for the first genome wide survey of common carp. The repetitive DNA was estimated to be approximate 28% of common carp genome, indicating the higher complexity of the genome. Comparative analysis had mapped around 40,000 BES to zebrafish genome and established over 3,100 microsyntenies, covering over 50% of the zebrafish genome. BES of common carp are tremendous tools for comparative mapping between the two closely related species, zebrafish and common carp, which should facilitate both structural and functional genome analysis in common carp.