Syntenic relationships between Medicago truncatula and Arabidopsis reveal extensive divergence of genome organization

Arabidopsis and Medicago truncatula represent sister clades within the dicot subclass Rosidae. We used genetic map-based and bacterial artificial chromosome sequence-based approaches to estimate the level of synteny between the genomes of these model plant species. Mapping of 82 tentative orthologous gene pairs reveals a lack of extended macrosynteny between the two genomes, although marker collinearity is frequently observed over small genetic intervals. Divergence estimates based on non-synonymous nucleotide substitutions suggest that a majority of the genes under analysis have experienced duplication in Arabidopsis subsequent to divergence of the two genomes, potentially confounding synteny analysis. Moreover, in cases of localized synteny, genetically linked loci in M. truncatula often share multiple points of synteny with Arabidopsis; this latter observation is consistent with the large number of segmental duplications that compose the Arabidopsis genome. More detailed analysis, based on complete sequencing and annotation of three M. truncatula bacterial artificial chromosome contigs suggests that the two genomes are related by networks of microsynteny that are often highly degenerate. In some cases, the erosion of microsynteny could be ascribed to the selective gene loss from duplicated loci, whereas in other cases, it is due to the absence of close homologs of M. truncatula genes in Arabidopsis.     

Chromosome evolution of MMU16 and RNO11: conserved synteny associated with gene order rearrangements explicable by intrachromosomal recombinations and neocentromere emergence

Comparative mapping between the rat and mouse genomes has shown that some chromosomes are entirely or almost entirely conserved with respect to gene content. Such is the case of rat chromosome 11 (RNO11) and mouse chromosome 16 (MMU16). We determined to what extent such an extensive conservation of synteny is associated with a conserved gene order. Therefore, we regionally localized several genes on RNO11. The comparison of the gene map of RNO11 and MMU16 unambiguously shows that the gene order has not been conserved in the Murinae lineage, thereby implying the occurrence of intrachromosomal evolutionary rearrangements. The transition from one chromosome configuration to the other one can be explained either by two intrachromosomal recombinations or by a single intrachromosomal recombination accompanied by neocentromere emergence.     

A COSII genetic map of the pepper genome provides a detailed picture of synteny with tomato and new insights into recent chromosome evolution in the genus Capsicum

We report herein the development of a pepper genetic linkage map which comprises 299 orthologous markers between the pepper and tomato genomes (including 263 conserved ortholog set II or COSII markers). The expected position of additional 288 COSII markers was inferred in the pepper map via pepper–tomato synteny, bringing the total orthologous markers in the pepper genome to 587. While pepper maps have been previously reported, this is the first complete map in the sense that all markers could be placed in 12 linkage groups corresponding to the 12 chromosomes. The map presented herein is relevant to the genomes of cultivated C. annuum and wild C. annuum (as well as related Capsicum species) which differ by a reciprocal chromosome translocation. This map is also unique in that it is largely based on COSII markers, which permits the inference of a detailed syntenic relationship between the pepper and tomato genomes—shedding new light on chromosome evolution in the Solanaceae. Since divergence from their last common ancestor is approximately 20 million years ago, the two genomes have become differentiated by a minimum number of 19 inversions and 6 chromosome translocations, as well as numerous putative single gene transpositions. Nevertheless, the two genomes share 35 conserved syntenic segments (CSSs) within which gene/marker order is well preserved. The high resolution COSII synteny map described herein provides a platform for cross-reference of genetic and genomic information (including the tomato genome sequence) between pepper and tomato and therefore will facilitate both applied and basic research in pepper. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Cloning and mapping of platypus SOX2 and SOX14: insights into SOX group B evolution

Group B SOX genes, the closest relatives to the sex-determining gene SRY, are thought to have evolved from a single ancestral SOX B by a series of duplications and translocations. The two SOX B genes SOX2 and SOX14 co-localize to chromosome 3q in humans. SOX2 and SOX14 homologues were cloned and characterized in the platypus, a monotreme mammal distantly related to man. The two genes were found to co-localize to chromosome 1q in this species. Proximity of the two related genes has therefore been conserved for 170 Myr, since humans and platypus diverged. The sequence similarity and conserved synteny of these group B genes provide clues to their origin. A simple model of SOX group B gene evolution is proposed.     

Comparative mapping of Z-orthologous genes in vertebrates: implications for the evolution of avian sex chromosomes

Sex chromosomes of birds and mammals are highly differentiated and share several cytological features. However, comparative gene mapping reveals extensive conserved synteny between the chicken Z sex chromosome and human chromosome 9 but not the human X sex chromosome, implying an independent origin of avian and mammalian sex chromosomes. To better understand the evolution of the avian Z chromosome we analysed the synteny of chicken Z-linked genes in zebrafish, which is the best-mapped teleost genome so far. Existing zebrafish maps do not support the existence of an ancestral Z linkage group in the zebrafish genome, whereas mammalian X-linked genes show at least some degree of synteny conservation. This is consistent with in situ hybridisation mapping data in the freshwater pufferfish, Tetraodon nigroviridis where mammalian X-linked genes show a much higher degree of conserved synteny than human chromosome 9 or the avian Z chromosome. Collectively, these data argue in favour of a more recent evolution of the avian Z chromosome, compared with the mammalian X.     

Conservation of the MHC-like region throughout evolution

Identification of conserved regions between the genomes of distant species is a crucial step in the reconstruction of the genomic organization of their last common ancestor. Here we confirm for the first time with robust evidence, the existence of a region of conserved synteny between the human genome and the Drosophila genome. This evolutionarily conserved synteny involves the human MHC and paralogous regions, and we identified 19 conserved genes between these two species in a Drosophila genomic region of less than 2 Mb. The statistical analysis of the distribution of these 19 genes between the Drosophila and human genomes shows that it cannot be explained by chance. Our study constitutes a first step towards the reconstruction of the genome of Urbilateria (the ancestor of all bilaterian) and allows for a better understanding of the evolutionary history of our genome as well as other metazoan genomes.     

The fragile breakage versus random breakage models of chromosome evolution

For many years, studies of chromosome evolution were dominated by the random breakage theory, which implies that there are no rearrangement hot spots in the human genome. In 2003, Pevzner and Tesler argued against the random breakage model and proposed an alternative “fragile breakage” model of chromosome evolution. In 2004, Sankoff and Trinh argued against the fragile breakage model and raised doubts that Pevzner and Tesler provided any evidence of rearrangement hot spots. We investigate whether Sankoff and Trinh indeed revealed a flaw in the arguments of Pevzner and Tesler. We show that Sankoff and Trinh's synteny block identification algorithm makes erroneous identifications even in small toy examples and that their parameters do not reflect the realities of the comparative genomic architecture of human and mouse. We further argue that if Sankoff and Trinh had fixed these problems, their arguments in support of the random breakage model would disappear. Finally, we study the link between rearrangements and regulatory regions and argue that long regulatory regions and inhomogeneity of gene distribution in mammalian genomes may be responsible for the breakpoint reuse phenomenon.     

The Easter Egg Weevil (Pachyrhynchus) genome reveals syntenic patterns in Coleoptera across 200 million years of evolution

Patterns of genomic architecture across insects remain largely undocumented or decoupled from a broader phylogenetic context. For instance, it is unknown whether translocation rates differ between insect orders. We address broad scale patterns of genome architecture across Insecta by examining synteny in a phylogenetic framework from open-source insect genomes. To accomplish this, we add a chromosome level genome to a crucial lineage, Coleoptera. Our assembly of the Pachyrhynchus sulphureomaculatus genome is the first chromosome scale genome for the hyperdiverse Phytophaga lineage and currently the largest insect genome assembled to this scale. The genome is significantly larger than those of other weevils, and this increase in size is caused by repetitive elements. Our results also indicate that, among beetles, there are instances of long-lasting (>200 Ma) localization of genes to a particular chromosome with few translocation events. While some chromosomes have a paucity of translocations, intra-chromosomal synteny was almost absent, with gene order thoroughly shuffled along a chromosome. This large amount of reshuffling within chromosomes with few inter-chromosomal events contrasts with patterns seen in mammals in which the chromosomes tend to exchange larger blocks of material more readily. To place our findings in an evolutionary context, we compared syntenic patterns across Insecta in a phylogenetic framework. For the first time, we find that synteny decays at an exponential rate relative to phylogenetic distance. Additionally, there are significant differences in decay rates between insect orders, this pattern was not driven by Lepidoptera alone which has a substantially different rate.     

Syntenic assignment of dopamine tautomerase (DCT) to bovine chromosome 12

Heterologous primers were used to amplify an exon and intron-containing segment of the bovine homologue of the human dopachrome tautomerase gene. After confirmation of homo-logy by sequence analysis (exon sequence similarity greater than 90%), bovine-specific primers were developed for synteny mapping purposes. The dopachrome tautomerase gene was assigned to bovine chromosome 12 (BTA12) with 97% concordance to the coagulation factor 10 locus. Together with previous synteny mapping of bovine chromosome 12 genes, fms-related tyrosine kinase, esterase D and 5-hydroxytryptamine receptor 2, this assignment further indicates conservation between human chromosome 13q and bovine chromosome 12.     

Sequencing and analysis of a 63 kb bacterial artificial chromosome insert from the Wolbachia endosymbiont of the human filarial parasite Brugia malayi.

Wolbachia endosymbiotic bacteria are widespread in filarial nematodes and are directly involved in the immune response of the host. In addition, antibiotics which disrupt Wolbachia interfere with filarial nematode development thus, Wolbachia provide an excellent target for control of filariasis. A 63.1 kb bacterial artificial chromosome insert, from the Wolbachia endosymbiont of the human filarial parasite Brugia malayi, has been sequenced using the New England Biolabs Inc. Genome Priming System() transposition kit in conjunction with primer walking methods. The bacterial artificial chromosome insert contains approximately 57 potential ORFs which have been compared by individual protein BLAST analysis with the 35 published complete microbial genomes in the Comprehensive Microbial Resource database at The Institute for Genomic Research and in the NCBI GenBank database, as well as to data from 22 incomplete genomes from the DOE Joint Genome Institute. Twenty five of the putative ORFs have significant similarity to genes from the alpha-proteobacteria Rickettsia prowazekii, the most closely related completed genome, as well as to the newly sequenced alpha-proteobacteria endosymbiont Sinorhizobium meliloti. The bacterial artificial chromosome insert sequence however has little conserved synteny with the R. prowazekii and S. meliloti genomes. Significant sequence similarity was also found in comparisons with the currently available sequence data from the Wolbachia endosymbiont of Drosophila melanogaster. Analysis of this bacterial artificial chromosome insert provides useful gene density and comparative genomic data that will contribute to whole genome sequencing of Wolbachia from the B. malayi host. This will also lead to a better understanding of the interactions between the endosymbiont and its host and will offer novel approaches and drug targets for elimination of filarial disease.     

Localization of anchor loci representing five hundred annotated rice genes to wheat chromosomes using PLUG markers

PCR-based Landmark Unique Gene (PLUG) markers are EST-PCR markers developed based on the orthologous gene conservation between rice and wheat, and on the intron polymorphisms among the three orthologous genes derived from the A, B and D genomes of wheat. We designed a total of 960 primer sets from wheat ESTs that showed high similarity with 951 single-copy rice genes. When genomic DNA of Chinese Spring wheat was used as a template, 872 primer sets amplified one to five distinct products. Out of these 872 PLUG markers, 531 were assigned to one or more chromosomes by nullisomic-tetrasomic analysis. For each wheat chromosome, the number of loci detected ranged from 32 for chromosome 6A to 73 for chromosome 7D, with an average of 48 loci per chromosome. Several novel synteny perturbations were identified using deletion bin-mapping of markers. Furthermore, we demonstrated that PLUG markers can be used as probes to simultaneously identify BAC clones that contain homoeologous regions from all three genomes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Gene space dynamics during the evolution of Aegilops tauschii, Brachypodium distachyon, Oryza sativa, and Sorghum bicolor genomes

Nine different regions totaling 9.7 Mb of the 4.02 Gb Aegilops tauschii genome were sequenced using the Sanger sequencing technology and compared with orthologous Brachypodium distachyon, Oryza sativa (rice), and Sorghum bicolor (sorghum) genomic sequences. The ancestral gene content in these regions was inferred and used to estimate gene deletion and gene duplication rates along each branch of the phylogenetic tree relating the four species. The total gene number in the extant Ae. tauschii genome was estimated to be 36,371. The gene deletion and gene duplication rates and total gene numbers in the four genomes were used to estimate the total gene number in each node of the phylogenetic tree. The common ancestor of the Brachypodieae and Triticeae lineages was estimated to have had 28,558 genes, and the common ancestor of the Panicoideae, Ehrhartoideae, and Pooideae subfamilies was estimated to have had 27,152 or 28,350 genes, depending on the ancestral gene scenario. Relative to the Brachypodieae and Triticeae common ancestor, the gene number was reduced in B. distachyon by 3,026 genes and increased in Ae. tauschii by 7,813 genes. The sum of gene deletion and gene duplication rates, which reflects the rate of gene synteny loss, was correlated with the rate of structural chromosome rearrangements and was highest in the Ae. tauschii lineage and lowest in the rice lineage. The high rate of gene space evolution in the Ae. tauschii lineage accounts for the fact that, contrary to the expectations, the level of synteny between the phylogenetically more related Ae. tauschii and B. distachyon genomes is similar to the level of synteny between the Ae. tauschii genome and the genomes of the less related rice and sorghum. The ratio of gene duplication to gene deletion rates in these four grass species closely parallels both the total number of genes in a species and the overall genome size. Because the overall genome size is to a large extent a function of the repeated sequence content in a genome, we suggest that the amount and activity of repeated sequences are important factors determining the number of genes in a genome.     

Mapping of the KHSRP gene to a region of conserved synteny on human chromosome 19p13.3 and mouse chromosome 17

The K homology-type splicing regulatory protein, KSRP, activates splicing through intronic splicing enhancer sequences. It is highly expressed in neural cells and is required for the neural-specific splicing of the c-src N1 exon. In this study, we mapped the gene (gene symbols KHSRP and Khsrp) to human chromosome 19 by using radiation hybrid panels and to mouse chromosome 17 by studying an interspecific backcross panel. Human KHSRP is a positional candidate gene for familial febrile convulsion and Cayman type cerebellar ataxia. Comparative analysis of the human and mouse genomes indicates that the KHSRP gene is located in regions of conserved synteny between the two species.     

Detailed Comparative Mapping of Cereal Chromosome Regions Corresponding to the Ph1 Locus in Wheat

Detailed physical mapping of markers from rice chromosome 9, and from syntenous (at the genetic level) regions of other cereal genomes, has resulted in rice yeast artificial chromosome (YAC) contigs spanning parts of rice 9. This physical mapping, together with comparative genetic mapping, has demonstrated that synteny has been largely maintained between the genomes of several cereals at the level of contiged YACs. Markers located in one region of rice chromosome 9 encompassed by the YAC contigs have exhibited restriction fragment length polymorphism (RFLP) using deletion lines for the Ph1 locus. This has allowed demarcation of the region of rice chromosome 9 syntenous with the ph1b and ph1c deletions in wheat chromosome 5B. A group of probes located in wheat homoeologous group 5 and barley chromosome 5H, however, have synteny with rice chromosomes other than 9. This suggests that the usefulness of comparative trait analysis and of the rice genome as a tool to facilitate gene isolation will differ from one region to the next, and implies that the rice genome is more ancestral in structure than those of the Triticeae.     

Comparative mapping of human claudin-1 (CLDN1) in great apes

The gene encoding claudin-1 (CLDN1) has been mapped to human chromosome 3 (HSA3; 3q28-->q29) using a radiation hybrid panel. Employing fluorescence in situ hybridization (FISH) we here show that a human P1-derived artificial chromosome (PAC) containing CLDN1 detects the orthologous sites in chromosomes of the great apes, chimpanzee, gorilla, and orangutan. Furthermore, the chromosomal position of CLDN1 was determined in mouse chromosomes by FISH. The position of fluorescent signals is confined to a single chromosomal site in both great apes and mouse and in each case maps to the chromosomal region that has conserved synteny with HSA3 (PTR2q28, GGO2q28, PPY2q38 and MMU16B1). Using a gene-specific probe our results are consistent with reports of the striking similarity of great ape and human genomes as illustrated previously by chromosome painting.     

Resolution of conflicting assignments for the bovine casein kinase II alpha (CSNK2A2) gene

The casein kinase II alpha' gene (CSNK2A2), which physically maps to human chromosome 16 (HSA16), has previously been mapped to bovine chromosome 5 (BTA5). Based on these results, a new segment of homology between the human and bovine genomes was suggested. In this paper we demonstrate linkage between CSNK2A2 and several markers on BTA18. Our result is supported by the extensive conservation of synteny between HSA16q and BTA18. Bovine chromosome 18 markers used in this study included several microsatellites, as well as the MC1R gene previously mapped to HSA16q24.3. Sequencing of the PCR-fragment mapped to BTA5 reveals that a CSNK-like retroposon was responsible for the conflicting assignments. The present results further extend the observed conservation of synteny between HSA16q and BTA18.     

Seventy million years of concerted evolution of a homoeologous chromosome pair, in parallel, in major Poaceae lineages

Whole genome duplication ~70 million years ago provided raw material for Poaceae (grass) diversification. Comparison of rice (Oryza sativa), sorghum (Sorghum bicolor), maize (Zea mays), and Brachypodium distachyon genomes revealed that one paleo-duplicated chromosome pair has experienced very different evolution than all the others. For tens of millions of years, the two chromosomes have experienced illegitimate recombination that has been temporally restricted in a stepwise manner, producing structural stratification in the chromosomes. These strata formed independently in different grass lineages, with their similarities (low sequence divergence between paleo-duplicated genes) preserved in parallel for millions of years since the divergence of these lineages. The pericentromeric region of this homeologous chromosome pair accounts for two-thirds of the gene content differences between the modern chromosomes. Both intriguing and perplexing is a distal chromosomal region with the greatest DNA similarity between surviving duplicated genes but also with the highest concentration of lineage-specific gene pairs found anywhere in these genomes and with a significantly elevated gene evolutionary rate. Intragenomic similarity near this chromosomal terminus may be important in hom(e)ologous chromosome pairing. Chromosome structural stratification, together with enrichment of autoimmune response-related (nucleotide binding site-leucine-rich repeat) genes and accelerated DNA rearrangement and gene loss, confer a striking resemblance of this grass chromosome pair to the sex chromosomes of other taxa.     

HSA4 and GGA4: remarkable conservation despite 300-Myr divergence

The chicken (GGA) and human (HSA) genomes diverged around 300-350 Myr ago. Due to this large phylogenetic distance, significant synteny conservation has not been anticipated between the genomes of the two species. However, Zoo-FISH with HSA4 chromosome-specific paint on chicken metaphase chromosomes shows that the human chromosome corresponds largely to the GGA4cen-->q26 region. Comparative gene mapping data in the two species, though limited, provide strong support for these observations. The findings, together with the very recently published data on HSA9-GGAZ and HSA12-GGA1, show that some large chromosomal segments share conserved synteny in the two species. These syntenies are considerably disrupted in the mouse. This makes us believe that despite very early divergence, parts of the human and chicken genomes are more conserved than those of human and mouse, which radiated only 100-120 Myr ago. Moreover, the HSA4-GGA4q correspondence points to a "candidate" chromosome from the karyotype of a mammal-bird ancestor. The findings are thus a small but important step toward understanding the evolution of the two genomes.     

A semi-quantitative, synteny-based method to improve functional predictions for hypothetical and poorly annotated bacterial and archaeal genes

During microbial evolution, genome rearrangement increases with increasing sequence divergence. If the relationship between synteny and sequence divergence can be modeled, gene clusters in genomes of distantly related organisms exhibiting anomalous synteny can be identified and used to infer functional conservation. We applied the phylogenetic pairwise comparison method to establish and model a strong correlation between synteny and sequence divergence in all 634 available Archaeal and Bacterial genomes from the NCBI database and four newly assembled genomes of uncultivated Archaea from an acid mine drainage (AMD) community. In parallel, we established and modeled the trend between synteny and functional relatedness in the 118 genomes available in the STRING database. By combining these models, we developed a gene functional annotation method that weights evolutionary distance to estimate the probability of functional associations of syntenous proteins between genome pairs. The method was applied to the hypothetical proteins and poorly annotated genes in newly assembled acid mine drainage Archaeal genomes to add or improve gene annotations. This is the first method to assign possible functions to poorly annotated genes through quantification of the probability of gene functional relationships based on synteny at a significant evolutionary distance, and has the potential for broad application.