Syntenic Sequences Research Articles

De Novo Assembly and Comparative Analysis of Mitochondrial Genomes of Two Pueraria montana Varieties.

Pueraria montana is a species with important medicinal value and a complex genetic background. In this study, we sequenced and assembled the mitochondrial (mt) genomes of two varieties of P. montana. The mt genome lengths of P. montana var. thomsonii and P. montana var. montana were 457,390 bp and 456,731 bp, respectively. Both P. montana mitogenomes showed a multi-branched structure consisting of two circular molecules, with 56 genes annotated, comprising 33 protein-coding genes, 18 tRNA genes (trnC-GCA and trnM-CAU are multi-copy genes), and 3 rRNA genes. Then, 207 pairs of long repeats and 96 simple sequence repeats (SSRs) were detected in the mt genomes of P. montana, and 484 potential RNA-editing sites were found across the 33 mitochondrial protein-coding genes of each variety. Additionally, a syntenic sequence analysis showed a high collinearity between the two mt genomes. This work is the first to analyze the mt genomes of P. montana. It can provide information that can be used to analyze the structure of mt genomes of higher plants and provide a foundation for future comparative genomic studies and evolutionary biology research in related species.

Evolutionary insights into 3D genome organization and epigenetic landscape of Vigna mungo.

Eukaryotic genomes show an intricate three-dimensional (3D) organization within the nucleus that regulates multiple biological processes including gene expression. Contrary to animals, understanding of 3D genome organization in plants remains at a nascent stage. Here, we investigate the evolution of 3D chromatin architecture in legumes. By using cutting-edge PacBio, Illumina, and Hi-C contact reads, we report a gap-free, chromosome-scale reference genome assembly of Vigna mungo, an important minor legume cultivated in Southeast Asia. We spatially resolved V. mungo chromosomes into euchromatic, transcriptionally active A compartment and heterochromatic, transcriptionally-dormant B compartment. We report the presence of TAD-like-regions throughout the diagonal of the HiC matrix that resembled transcriptional quiescent centers based on their genomic and epigenomic features. We observed high syntenic breakpoints but also high coverage of syntenic sequences and conserved blocks in boundary regions than in the TAD-like region domains. Our findings present unprecedented evolutionary insights into spatial 3D genome organization and epigenetic patterns and their interaction within the V. mungo genome. This will aid future genomics and epigenomics research and breeding programs of V. mungo.

Comparative analysis of mitochondrial genomes of two alpine medicinal plants of Gentiana (Gentianaceae).

Gentiana crassicaulis and G. straminea are alpine plants of Gentiana with important medicinal value and complex genetic backgrounds. In this study, the mitochondrial genomes (mtDNAs) of these two species were sequenced. The mtDNAs of G. crassicaulis and G. straminea are 368,808 and 410,086 bp long, respectively, 52 and 49 unique genes are annotated in the two species, and the gene arrangement varies widely. Compared to G. crassicaulis, G. straminea loses three effective genes, namely atp6, trnG-GCC and trnV-GAC. As a pseudogene, the atp6 gene of G. straminea is incomplete, which is rare in higher plants. We detected 1696 and 1858 pairs of long repeats and 213 SSRs and 250 SSs in the mtDNAs of G. crassicaulis and G. straminea, respectively. There are 392 SNPs and 18 InDels between the two genomes, and syntenic sequence and structural variation analysis show low collinearity between the two genomes. Chloroplast DNA transferring to mtDNA is observed in both species, and 46,511 and 55,043 bp transferred segments containing three tRNA genes are identified, respectively. Comparative analysis of mtDNAs of G. crassicaulis, G. straminea and four species of Gentianales determined 18 core genes, and there is no specific gene in G. crassicaulis and G. straminea. The phylogenetic tree based on mtDNAs places Gentianaceae in a branch of Gentianales. This study is the first to analyze the mtDNAs of Gentianaceae, which could provide information for analysis of the structure of mtDNAs of higher plants and phylogenetic research of Gentianaceae and Gentianales.

Reactivation of transposable elements following hybridization in fission yeast.

Hybridization is thought to reactivate transposable elements (TEs) that were efficiently suppressed in the genomes of the parental hosts. Here, we provide evidence for this “genomic shock hypothesis” in the fission yeast Schizosaccharomyces pombe. In this species, two divergent lineages (Sp and Sk) have experienced recent, likely human-induced, hybridization. We used long-read sequencing data to assemble genomes of 37 samples derived from 31 S. pombe strains spanning a wide range of ancestral admixture proportions. A comprehensive TE inventory revealed exclusive presence of long terminal repeat (LTR) retrotransposons. Sequence analysis of active full-length elements, as well as solo LTRs, revealed a complex history of homologous recombination. Population genetic analyses of syntenic sequences placed insertion of many solo LTRs before the split of the Sp and Sk lineages. Most full-length elements were inserted more recently, after hybridization. With the exception of a single full-length element with signs of positive selection, both solo LTRs and, in particular, full-length elements carry signatures of purifying selection indicating effective removal by the host. Consistent with reactivation upon hybridization, the number of full-length LTR retrotransposons, varying extensively from zero to 87 among strains, significantly increases with the degree of genomic admixture. This study gives a detailed account of global TE diversity in S. pombe, documents complex recombination histories within TE elements, and provides evidence for the “genomic shock hypothesis.”

Chromosome-level genome assembly of a regenerable maize inbred line A188

BackgroundThe maize inbred line A188 is an attractive model for elucidation of gene function and improvement due to its high embryogenic capacity and many contrasting traits to the first maize reference genome, B73, and other elite lines. The lack of a genome assembly of A188 limits its use as a model for functional studies.ResultsHere, we present a chromosome-level genome assembly of A188 using long reads and optical maps. Comparison of A188 with B73 using both whole-genome alignments and read depths from sequencing reads identify approximately 1.1 Gb of syntenic sequences as well as extensive structural variation, including a 1.8-Mb duplication containing the Gametophyte factor1 locus for unilateral cross-incompatibility, and six inversions of 0.7 Mb or greater. Increased copy number of carotenoid cleavage dioxygenase 1 (ccd1) in A188 is associated with elevated expression during seed development. High ccd1 expression in seeds together with low expression of yellow endosperm 1 (y1) reduces carotenoid accumulation, accounting for the white seed phenotype of A188. Furthermore, transcriptome and epigenome analyses reveal enhanced expression of defense pathways and altered DNA methylation patterns of the embryonic callus.ConclusionsThe A188 genome assembly provides a high-resolution sequence for a complex genome species and a foundational resource for analyses of genome variation and gene function in maize. The genome, in comparison to B73, contains extensive intra-species structural variations and other genetic differences. Expression and network analyses identify discrete profiles for embryonic callus and other tissues.

Reconstructing the Clonal and Mutational Architecture of Myeloma through Avian Leukosis Virus (ALV)-Mediated Genome Editing

Background and Aims. Despite the fact that the myeloma mutational landscape was elucidated several years ago, there are yet no approved precision therapies targeting specific driver mutations. The role of common, but somewhat enigmatic, mutational events remains elusive, for example, the loss-of-function mutation of the RNA-processing exosome ribonuclease, DIS3, the third most-common mutation after KRAS and NRAS mutations. Furthermore, complex clonal architectural patterns (such as subclonal BRAF mutations nested within RAS-mutant clones) complicate the design of targeted therapy trials.To generate precision in vivo models that recapitulate the complex clonal and mutational architecture of myeloma, we integrated two powerful technological platforms: the first RAS-driven myeloma in vivo model (VQ model) created through a collaboration between the Zhang and Asimakopoulos labs coupled with flexible avian leukosis virus (ALV)-mediated gene transfer that permits myeloma genome editing in a temporally- and spatially- controlled fashion (Prdm1: TVb-mRFP model)(Asimakopoulos and Varmus, J Virol, 2009; 83(10):4835-43).VQ is based on the state-of-art Vk*MYC model (Chesi et al., Cancer Cell 2008;13(2):167-80). VQ was generated by conditionally activating an NrasQ61R allele in Vk*MYC mice. The resultant VQ model recapitulates an advanced, proliferative and drug-resistant phase of disease characterized by RAS mutations. The transducibility of VQ restored the functionality of the ALV-gene transfer technology. This technology exploits the fact that avian retroviruses cannot transduce mammalian cells unless they are engineered to ectopically express an avian retroviral receptor (in this case, the receptor TVb, driven by a plasma-cell-specific regulatory elements). Recently, the ALV gene transfer technology was adapted to achieve CRISPR-Cas9 mediated gene editing in the tumor niche in vivo (Oldrini et al., Nat Comm 2018;9(1):1466).Approach. We generated TVb-expressing VQ lines (CD138+ B220neg) that home to the bone marrow upon transplantation into Prdm1:TVb-mRFP fully immunocompetent syngeneic recipients in the C57BL6/J background. These VQ cells were engineered to express the nuclease Cas9. To determine the efficiency of gene transfer through ALV vectors (as a prelude to sgRNA delivery), pilot experiments were conducted using two types of ALV vector carrying a fluorescent reporter, GFP (green fluorescent protein): a) an RCAS-type vector that drives reporter expression from the native ALV LTR and b) an RCAN-type vector that drives reporter expression from an internal CMV promoter.Results. ALV titers achieved were routinely in excess of 107-108 infectious units/mL. Parental Vk*MYC cells expressing the avian receptor TVb were poorly transducible with ALV vectors. By contrast, VQ cells expressing TVb were readily transducible. At multiplicity of infection (MOI)=2, 68% of VQ myeloma cells expressed reporter fluorescence after a single round of transduction. At MOI=5 the transduction rate rose to 80%. The mean fluorescence intensity (MFI) positively correlated with MOI, consistent with higher-copy proviral integration at higher MOI. At equivalent MOI, GFP fluorescence from RCAS vectors was more intense than GFP fluorescence driven by an internal CMV promoter in RCAN vectors, suggesting that the ALV LTR is a stronger promoter than CMV in VQ myeloma cells. 100% transduced VQ cells were CD138+B220neg at all conditions tested.Conclusions and Current Focus. We are using ALV-CRISPR technology to generate precision myeloma models to support the next wave of personalized therapy and immunotherapy for this recalcitrant and incurable cancer.Progress on three in vivo models under construction will be presented at the meeting:•Subclonal Braf mutation nested within clonal Nras mutation. These models will be treated with next-generation RAF inhibitors (Karoulia Z et al.Nat Rev Cancer 2017;17(11):676-91).•Dis3 loss-of function mutation in constellation with co-clonal Ras mutation. The pathogenic consequences of Dis3 disruption remain poorly understood as no parallels exist in other cancers.•The murine homolog of the t(4;14) translocation. We are targeting breakpoint cluster regions within perfectly syntenic sequences on mouse chr. 5 (containing Fgfr3 and Mmset loci) and mouse chr.12 (containing the murine Igh locus). DisclosuresNo relevant conflicts of interest to declare.

Identification of a gene responsible for cytoplasmic male-sterility in onions (Allium cepa L.) using comparative analysis of mitochondrial genome sequences of two recently diverged cytoplasms.

Almost identical mitochondrial genome sequences of two recently diverged male-fertile normal and male-sterile CMS-T-like cytoplasms were obtained in onions. A chimeric gene, orf725 , was found to be a CMS-inducing gene. In onions (Allium cepa L.), cytoplasmic male-sterility (CMS) has been widely used in hybrid seed production. Two types of CMS (CMS-S and CMS-T) have been reported in onions. A complete mitochondrial genome sequence of the CMS-S cytoplasm has been reported in our previous study. Draft mitochondrial genome sequences of male-fertile normal and CMS-T-like cytoplasms are reported in this study. Raw reads obtained from normal and CMS-T-like cytoplasms were assembled into eight and nine almost identical contigs, respectively. After connection and reorganization of contigs by PCR amplification and genome walking, four scaffold sequences with total length of 339 and 180bp were produced for the normal cytoplasm. A mitochondrial genome sequence of the CMS-T-like cytoplasm was obtained by mapping trimmed reads of CMS-T onto scaffold sequences of the normal cytoplasm. Compared with the CMS-S mitochondrial genome, the normal mitochondrial genome was highly rearranged with 31 syntenic blocks. A total of 499 single nucleotide polymorphisms (SNPs) or insertions/deletions were identified in these syntenic regions. On the other hand, normal and CMS-T-like mitochondrial genome sequences were almost identical except for orf725, a chimeric gene consisting of cox1 with other sequences. Only three SNPs were identified between normal and CMS-T-like syntenic sequences. These results indicate that orf725 is likely to be the casual gene for CMS induction in onions and that CMS-T-like cytoplasm has recently diverged from the normal cytoplasm by introduction of orf725.

Construction of Red Fox Chromosomal Fragments from the Short-Read Genome Assembly.

The genome of a red fox (Vulpes vulpes) was recently sequenced and assembled using next-generation sequencing (NGS). The assembly is of high quality, with 94X coverage and a scaffold N50 of 11.8 Mbp, but is split into 676,878 scaffolds, some of which are likely to contain assembly errors. Fragmentation and misassembly hinder accurate gene prediction and downstream analysis such as the identification of loci under selection. Therefore, assembly of the genome into chromosome-scale fragments was an important step towards developing this genomic model. Scaffolds from the assembly were aligned to the dog reference genome and compared to the alignment of an outgroup genome (cat) against the dog to identify syntenic sequences among species. The program Reference-Assisted Chromosome Assembly (RACA) then integrated the comparative alignment with the mapping of the raw sequencing reads generated during assembly against the fox scaffolds. The 128 sequence fragments RACA assembled were compared to the fox meiotic linkage map to guide the construction of 40 chromosomal fragments. This computational approach to assembly was facilitated by prior research in comparative mammalian genomics, and the continued improvement of the red fox genome can in turn offer insight into canid and carnivore chromosome evolution. This assembly is also necessary for advancing genetic research in foxes and other canids.

Functionally conserved enhancers with divergent sequences in distant vertebrates.

BackgroundTo examine the contributions of sequence and function conservation in the evolution of enhancers, we systematically identified enhancers whose sequences are not conserved among distant groups of vertebrate species, but have homologous function and are likely to be derived from a common ancestral sequence. Our approach combined comparative genomics and epigenomics to identify potential enhancer sequences in the genomes of three groups of distantly related vertebrate species.ResultsWe searched for sequences that were conserved within groups of closely related species but not between groups of more distant species, and were associated with an epigenetic mark of enhancer activity. To facilitate inferring orthology between non-conserved sequences, we limited our search to introns whose orthology could be unambiguously established by mapping the bracketing exons. We show that a subset of these non-conserved but syntenic sequences from the mouse and zebrafish genomes have homologous functions in a zebrafish transgenic enhancer assay. The conserved expression patterns driven by these enhancers are probably associated with short transcription factor-binding motifs present in the divergent sequences.ConclusionsWe have identified numerous potential enhancers with divergent sequences but a conserved function. These results indicate that selection on function, rather than sequence, may be a common mode of enhancer evolution; evidence for selection at the sequence level is not a necessary criterion to define a gene regulatory element.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2070-7) contains supplementary material, which is available to authorized users.

A disease resistance locus on potato and tomato chromosome 4 exhibits a conserved multipartite structure displaying different rates of evolution in different lineages.

BackgroundIn plant genomes, NB-LRR based resistance (R) genes tend to occur in clusters of variable size in a relatively small number of genomic regions. R-gene sequences mostly differentiate by accumulating point mutations and gene conversion events. Potato and tomato chromosome 4 harbours a syntenic R-gene locus (known as the R2 locus in potato) that has mainly been examined in central American/Mexican wild potato species on the basis of its contribution to resistance to late blight, caused by the oomycete pathogen Phytophthora infestans. Evidence to date indicates the occurrence of a fast evolutionary mode characterized by gene conversion events at the locus in these genotypes.ResultsA physical map of the R2 locus was developed for three Solanum tuberosum genotypes and used to identify the tomato syntenic sequence. Functional annotation of the locus revealed the presence of numerous resistance gene homologs (RGHs) belonging to the R2 gene family (R2GHs) organized into a total of 4 discrete physical clusters, three of which were conserved across S. tuberosum and tomato. Phylogenetic analysis showed clear orthology/paralogy relationships between S. tuberosum R2GHs but not in R2GHs cloned from Solanum wild species. This study confirmed that, in contrast to the wild species R2GHs, which have evolved through extensive sequence exchanges between paralogs, gene conversion was not a major force for differentiation in S. tuberosum R2GHs, and orthology/paralogy relationships have been maintained via a slow accumulation of point mutations in these genotypes.ConclusionsS. tuberosum and Solanum lycopersicum R2GHs evolved mostly through duplication and deletion events, followed by gradual accumulation of mutations. Conversely, widespread gene conversion is the major evolutionary force that has shaped the locus in Mexican wild potato species. We conclude that different selective forces shaped the evolution of the R2 locus in these lineages and that co-evolution with a pathogen steered selection on different evolutionary paths.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-015-0645-8) contains supplementary material, which is available to authorized users.

Interspecific and intraspecific gene variability in a 1-Mb region containing the highest density of NBS-LRR genes found in the melon genome.

BackgroundPlant NBS-LRR -resistance genes tend to be found in clusters, which have been shown to be hot spots of genome variability. In melon, half of the 81 predicted NBS-LRR genes group in nine clusters, and a 1 Mb region on linkage group V contains the highest density of R-genes and presence/absence gene polymorphisms found in the melon genome. This region is known to contain the locus of Vat, an agronomically important gene that confers resistance to aphids. However, the presence of duplications makes the sequencing and annotation of R-gene clusters difficult, usually resulting in multi-gapped sequences with higher than average errors.ResultsA 1-Mb sequence that contains the largest NBS-LRR gene cluster found in melon was improved using a strategy that combines Illumina paired-end mapping and PCR-based gap closing. Unknown sequence was decreased by 70% while about 3,000 SNPs and small indels were corrected. As a result, the annotations of 18 of a total of 23 NBS-LRR genes found in this region were modified, including additional coding sequences, amino acid changes, correction of splicing boundaries, or fussion of ORFs in common transcription units. A phylogeny analysis of the R-genes and their comparison with syntenic sequences in other cucurbits point to a pattern of local gene amplifications since the diversification of cucurbits from other families, and through speciation within the family. A candidate Vat gene is proposed based on the sequence similarity between a reported Vat gene from a Korean melon cultivar and a sequence fragment previously absent in the unrefined sequence.ConclusionsA sequence refinement strategy allowed substantial improvement of a 1 Mb fragment of the melon genome and the re-annotation of the largest cluster of NBS-LRR gene homologues found in melon. Analysis of the cluster revealed that resistance genes have been produced by sequence duplication in adjacent genome locations since the divergence of cucurbits from other close families, and through the process of speciation within the family a candidate Vat gene was also identified using sequence previously unavailable, which demonstrates the advantages of genome assembly refinements when analyzing complex regions such as those containing clusters of highly similar genes.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-1131) contains supplementary material, which is available to authorized users.

A genetic linkage map of kelp grouper (Epinephelus bruneus) based on microsatellite markers

Kelp grouper (Epinephelus bruneus) is an important aquaculture species in Japan, Korea and China. Mariculture production of the species has increased due to its high demand and market price. However, some problems affect kelp grouper aquaculture such as the low growth rate, high mortality due to diseases and low ability to survive the larval stage. To analyze economically important traits, genetic linkage maps are an effective tool. We constructed sex-specific linkage maps of kelp grouper using 222 microsatellite markers. The male map consisted of 23 linkage groups with 161 markers and the female map consisted of 25 linkage groups with 173 markers. The total lengths of the male and female maps were 650.5cM and 944.4cM, respectively, and the average intervals were 5.0cM and 6.7cM, respectively. The average ratio of recombination between males and females was 1:1.5. Moreover, syntenic sequence comparisons provided basic information of several potential candidate genes affecting organism physiological and biochemical reactions. Based on the linkage map, further quantitative trait loci (QTL) or candidate gene(s) detection can be anticipated to contribute to assist breeding programs of kelp grouper. In addition, by providing basic genome information of kelp grouper, the map provides a first step towards comparative QTL as well as comparative genome analyses with other groupers in the future.

Complete mitochondrial genome sequence and identification of a candidate gene responsible for cytoplasmic male sterility in radish (Raphanus sativus L.) containing DCGMS cytoplasm

A novel cytoplasmic male sterility (CMS) conferred by Dongbu cytoplasmic and genic male-sterility (DCGMS) cytoplasm and its restorer-of-fertility gene (Rfd1) was previously reported in radish (Raphanus sativus L.). Its inheritance of fertility restoration and profiles of mitochondrial DNA (mtDNA)-based molecular markers were reported to be different from those of Ogura CMS, the first reported CMS in radish. The complete mitochondrial genome sequence (239,186bp; GenBank accession No. KC193578) of DCGMS mitotype is reported in this study. Thirty-four protein-coding genes and three ribosomal RNA genes were identified. Comparative analysis of a mitochondrial genome sequence of DCGMS and previously reported complete sequences of normal and Ogura CMS mitotypes revealed various recombined structures of seventeen syntenic sequence blocks. Short-repeat sequences were identified in almost all junctions between syntenic sequence blocks. Phylogenetic analysis of three radish mitotypes showed that DCGMS was more closely related to the normal mitotype than to the Ogura mitotype. A single 1,551-bp unique region was identified in DCGMS mtDNA sequences and a novel chimeric gene, designated orf463, consisting of 128-bp partial sequences of cox1 gene and 1,261-bp unidentified sequences were found in the unique region. No other genes with a chimeric structure, a major feature of most characterized CMS-associated genes in other plant species, were found in rearranged junctions of syntenic sequence blocks. Like other known CMS-associated mitochondrial genes, the predicted gene product of orf463 contained 12 transmembrane domains. Thus, this gene product might be integrated into the mitochondrial membrane. In total, the results indicate that orf463 is likely to be a casual factor for CMS induction in radish containing the DCGMS cytoplasm.

Just One Copy

Most model plant species have undergone at least one whole-genome duplication event relative to their nearest well-studied relative. The occurrence of a duplicated genome in the history of the species, even if not retained in the extant species, can confound estimates of molecular evolution. Wang et al. have investigated the relative rates of molecular evolution in a conserved syntenic sequence in members of the Solanaceae family (eggplant, pepper, petunia, potato, and tomato) that have undergone no genome-wide duplication events since they diverged, which from this study was estimated to have been ∼30 million years ago. Although some regions were duplicated within individual species, the authors were able to analyze genes and regions that have, more or less, been evolving as a single copy within the genome. They identified small-scale differences in the location of insertions and deletions, lineage-specific selection, gene content, order, and orientation, and estimated that approximately one-third of the examined regions is under selection and two-thirds of the sequences under selection are outside of genes and not associated with domestication.— LMZ

Pericentromeric euchromatin is conserved in minute human supernumerary chromosomes: a study using cross-species colour segmenting (RxFISH).

Investigation of marker chromosomes is one of the most challenging areas of clinical cytogenetics, especially in the prenatal scenario. A range of techniques including microdissection/reverse painting, SKY and M-FISH are available for the investigation of larger markers (>3 Mb). All these techniques rely on hybridization of unique, homologous sequences with simultaneous suppression of repeat sequences. In contrast, RxFISH is based on hybridization of cross-species syntenic sequences; repeat sequences do not hybridize due to species divergence. We have used RxFISH to analyse a group of the smallest, i.e. minute, supernumerary marker chromosomes. Our results suggest that even the smallest marker chromosomes often contain conserved pericentric euchromatin. More detailed characterization of pericentric genetic content is needed to assess the clinical significance of minute supernumerary markers.

The TIGR rice genome annotation resource: annotating the rice genome and creating resources for plant biologists.

Rice is not only a major food staple for the world's population but it also is a model species for a major group of flowering plants, the monocotyledonous plants. Draft genomic sequence of two subspecies of rice, Oryza sativa spp. japonica and indica ssp. are publicly available. To provide the community with a resource to data-mine the rice genome, we have constructed an annotation resource for rice (http://www.tigr.org/tdb/e2k1/osa1/). In this resource, we have annotated the rice genome for gene content, identified motifs/domains within the predicted genes, constructed a rice repeat database, identified related sequences in other plant species, and identified syntenic sequences between rice and maize. All of the data is available through web-based interfaces, FTP downloads, and a Distributed Annotation System.

Alignment of whole genomes.

A new system for aligning whole genome sequences is described. Using an efficient data structure called a suffix tree, the system is able to rapidly align sequences containing millions of nucleotides. Its use is demonstrated on two strains of Mycoplasma tuberculosis, on two less similar species of Mycoplasma bacteria and on two syntenic sequences from human chromosome 12 and mouse chromosome 6. In each case it found an alignment of the input sequences, using between 30 s and 2 min of computation time. From the system output, information on single nucleotide changes, translocations and homologous genes can easily be extracted. Use of the algorithm should facilitate analysis of syntenic chromosomal regions, strain-to-strain comparisons, evolutionary comparisons and genomic duplications.

On the complexity and approximation of syntenic distance

The paper studies the computational complexity and approximation algorithms for a new evolutionary distance between multi-chromosomal genomes introduced recently by Ferretti, Nadeau and Sankoff. Here, a chromosome is represented as a set of genes and a genome is a collections of chromosomes. The syntenic distance between two genomes is defined as the minimum number of translocations, fusions and fissions required to transform one genome into the other. We prove that computing the syntenic distance is NP-hard and give a simple approximation algorithm with performance ratio 2. For the case when an upper bound d on the syntenic distance is known, we show that an optimal syntenic sequence can be found in O( nk + 2 o( d 2) ) time, where n and k are the number of chromosomes in the two given genomes. Next, we show that if the set of operations for transforming a genome is significantly restricted, we can nevertheless find a solution that performs at most O(log d) additional moves, where d is the number of moves performed by the unrestricted optimum. This result should help in the design of approximation algorithms. Finally, we investigate the median problem: Given three genomes, construct a genome minimizing the total syntenic distance to the three given genomes and compute the corresponding median distance. The problem has application in the inference of phytogenies based on the syntenic distance. We prove that the problem is NP-hard and design a polynomial time approximation algorithm with a performance ratio of 4+ ε for any constant ε > 0.