High resolution synteny maps allowing direct comparisons between the coffee and tomato genomes

Tomato (Solanum lycopersicum) and coffee (Coffea canephora) belong to the sister families Solanaceae and Rubiaceae, respectively. We report herein the mapping of a common set of 257 Conserved Ortholog Set II genes in the genomes of both species. The mapped markers are well distributed across both genomes allowing the first syntenic comparison between species from these two families. The majority (75%) of the synteny blocks are short (<4 cM); however, some extend up to 50 cM. In an effort to further characterize the synteny between these two genomes, we took advantage of the available sequence for the tomato genome to show that tomato chromosome 7 is syntenic to half of the two coffee linkage groups E and F with the putative break point in tomato localized to the boundary of the heterochromatin and euchromatin on the long arm. In addition to the new insight on genome conservation and evolution between the plant families Solanaceae and Rubiaceae, the comparative maps presented herein provide a translational tool by which coffee researchers may take benefit of DNA sequence and genetic information from tomato and vice versa. It is thus expected that these comparative genome information will help to facilitate and expedite genetic and genomic research in coffee.     

The JAX Synteny Browser for mouse-human comparative genomics

Visualizing regions of conserved synteny between two genomes is supported by numerous software applications. However, none of the current applications allow researchers to select genome features to display or highlight in blocks of synteny based on the annotated biological properties of the features (e.g., type, function, and/or phenotype association). To address this usability gap, we developed an interactive web-based conserved synteny browser, The Jackson Laboratory (JAX) Synteny Browser. The browser allows researchers to highlight or selectively display genome features in the reference and/or the comparison genome according to the biological attributes of the features. Although the current implementation for the browser is limited to the reference genomes for the laboratory mouse and human, the software platform is intentionally genome agnostic. The JAX Synteny Browser software can be deployed for any two genomes where genome coordinates for syntenic blocks are defined and for which biological attributes of the features in one or both genomes are available in widely used standard bioinformatics file formats. The JAX Synteny Browser is available at: http://syntenybrowser.jax.org/. The code base is available from GitHub: https://github.com/TheJacksonLaboratory/syntenybrowser and is distributed under the Creative Commons Attribution license (CC BY).

     

Screening synteny blocks in pairwise genome comparisons through integer programming

Background  

It is difficult to accurately interpret chromosomal correspondences such as true orthology and paralogy due to significant divergence of genomes from a common ancestor. Analyses are particularly problematic among lineages that have repeatedly experienced whole genome duplication (WGD) events. To compare multiple "subgenomes" derived from genome duplications, we need to relax the traditional requirements of "one-to-one" syntenic matchings of genomic regions in order to reflect "one-to-many" or more generally "many-to-many" matchings. However this relaxation may result in the identification of synteny blocks that are derived from ancient shared WGDs that are not of interest. For many downstream analyses, we need to eliminate weak, low scoring alignments from pairwise genome comparisons. Our goal is to objectively select subset of synteny blocks whose total scores are maximized while respecting the duplication history of the genomes in comparison. We call this "quota-based" screening of synteny blocks in order to appropriately fill a quota of syntenic relationships within one genome or between two genomes having WGD events.     

Comparative sequence analysis of the region harboring the hardness locus in barley and its colinear region in rice

The ancestral shared synteny concept has been advocated as an approach to positionally clone genes from complex genomes. However, the unified grass genome model and the study of grasses as a single syntenic genome is a topic of considerable controversy. Hence, more quantitative studies of cereal colinearity at the sequence level are required. This study compared a contiguous 300-kb sequence of the barley (Hordeum vulgare) genome with the colinear region in rice (Oryza sativa). The barley sequence harbors genes involved in endosperm texture, which may be the subject of distinctive evolutionary forces and is located at the extreme telomeric end of the short arm of chromosome 5H. Comparative sequence analysis revealed the presence of five orthologous genes and a complex, postspeciation evolutionary history involving small chromosomal rearrangements, a translocation, numerous gene duplications, and extensive transposon insertion. Discrepancies in gene content and microcolinearity indicate that caution should be exercised in the use of rice as a surrogate for map-based cloning of genes from large genome cereals such as barley.     

Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.     

Fine mapping and comparative genomics integration of two quantitative trait loci controlling resistance to powdery mildew in a Spanish barley landrace

The intervals containing two major quantitative trait loci (QTL) from a Spanish barley landrace conferring broad spectrum resistance to Blumeria graminis were subjected to marker saturation. First, all the available information on recently developed marker resources for barley was exploited. Then, a comparative genomic analysis of the QTL regions with other sequenced grass model species was performed. As a result of the first step, 32 new markers were added to the previous map and new flanking markers closer to both QTL were identified. Next, syntenic integration revealed that the barley target regions showed homology with regions on chromosome 6 of rice (Oryza sativa), chromosome 10 of Sorghum bicolor and chromosome 1 of Brachypodium distachyon. A nested insertion of ancestral syntenic blocks on Brachypodium chromosome 1 was confirmed. Based on sequence information of the most likely candidate orthologous genes, 23 new barley unigene-derived markers were developed and mapped within the barley target regions. The assessment of colinearity revealed an inversion on chromosome 7HL of barley compared to the other three grass species, and nearly perfect colinearity on chromosome 7HS. This two-step marker enrichment allowed for the refinement of the two QTL into much smaller intervals. Inspection of all predicted proteins for the barley unigenes identified within the QTL intervals did not reveal the presence of resistance gene candidates. This study demonstrates the usefulness of sequenced genomes for fine mapping and paves the way for the use of these two loci in barley breeding programs.     

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.     

Cinteny: flexible analysis and visualization of synteny and genome rearrangements in multiple organisms

Background  

Identifying syntenic regions, i.e., blocks of genes or other markers with evolutionary conserved order, and quantifying evolutionary relatedness between genomes in terms of chromosomal rearrangements is one of the central goals in comparative genomics. However, the analysis of synteny and the resulting assessment of genome rearrangements are sensitive to the choice of a number of arbitrary parameters that affect the detection of synteny blocks. In particular, the choice of a set of markers and the effect of different aggregation strategies, which enable coarse graining of synteny blocks and exclusion of micro-rearrangements, need to be assessed. Therefore, existing tools and resources that facilitate identification, visualization and analysis of synteny need to be further improved to provide a flexible platform for such analysis, especially in the context of multiple genomes.     

Development of conserved ortholog set markers linked to the restorer gene Rfp1 in rye

Restoration of male fertility is a prerequisite for hybrid rye breeding and currently the most straightforward approach to minimize ergot infection in hybrid rye varieties. Molecular markers are important tools for the efficient introgression and management of restorer genes like Rfp1 originating from unadapted genetic resources. Furthermore, closely linked markers flanking Rfp1 are indispensible for identifying and selecting individuals with haplotypes showing recombination between Rfp1 and other gene(s) that reside in close proximity and have a negative influence on yield. We identified orthologous gene sets in rice, Brachypodium, and Sorghum and used these gene models as templates to establish conserved ortholog set (COS) markers for the restorer gene Rfp1 on the long arm of rye chromosome 4R. The novel co-dominant markers delimit Rfp1 within a 0.7-cM interval and allow prediction of Rfp1 genotypes with a precision not feasible before. The COS markers enabled an alignment of the improved genetic map of rye chromosome 4R with wheat and barley maps and allowed identification of regions orthologous to Rfp1 in wheat and barley on the short arms of chromosomes 6D and 6H, respectively. Results obtained in this study revealed that micro-collinearity around the Rfp1 locus in rye is affected by rearrangements relative to other grass genomes. The impact of the novel COS markers for practical hybrid rye breeding is discussed.     

High-resolution mapping of the barley leaf rust resistance gene Rph5 using barley expressed sequence tags (ESTs) and synteny with rice

The rapidly growing expressed sequence tag (EST) resources of species representing the Poacea family and availability of comprehensive sequence information for the rice (Oryza sativa) genome create an excellent opportunity for comparative genome analysis. Extensive synteny between rice chromosome 1 and barley (Hordeum vulgare L.) chromosome 3 has proven extremely useful for saturation mapping of chromosomal regions containing target genes of large-genome barley with conserved orthologous genes from the syntenic regions of the rice genome. Rph5 is a gene conferring resistance to the barley leaf rust pathogen Puccinia hordei. It was mapped to chromosome 3HS, which is syntenic with rice chromosome 1S. The objective of this study was to increase marker density within the sub-centimorgan region around Rph5, using sequence-tagged site (STS) markers that were developed based on barley ESTs syntenic to the phage (P1)-derived artificial chromosome (PAC) clones comprising the distal region of rice chromosome 1S. Five rice PAC clones were used as queries in a blastn search to screen 375,187 barley ESTs. Ninety-four non-redundant EST sequences were identified from the EST database and used as templates to design 174 pairs of primer combinations. As a result, 9 barley EST-based STS markers were incorporated into the ‘Bowman’ × ‘Magnif 102’ high-resolution map of the Rph5 region. More importantly, six markers, including five EST-derived STS sequences, were found to co-segregate with Rph5. The results of this study demonstrate the usefulness of rice genomic resources for efficient deployment of barley ESTs for marker saturation of targeted barley genomic regions.     

Towards a whole‐genome sequence for rye (Secale cereale L.)

We report on a whole‐genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole‐genome shotgun sequencing of the 7.9‐Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low‐copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high‐density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high‐density genetic map and establishing a synteny‐based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome‐wide map of selection signals differentiating the divergent gene pools. This rye whole‐genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome‐based breeding in rye.     

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.     

Sequence and structure of Brassica rapa chromosome A3

Background

The species Brassica rapa includes important vegetable and oil crops. It also serves as an excellent model system to study polyploidy-related genome evolution because of its paleohexaploid ancestry and its close evolutionary relationships with Arabidopsis thaliana and other Brassica species with larger genomes. Therefore, its genome sequence will be used to accelerate both basic research on genome evolution and applied research across the cultivated Brassica species.

Results

We have determined and analyzed the sequence of B. rapa chromosome A3. We obtained 31.9 Mb of sequences, organized into nine contigs, which incorporated 348 overlapping BAC clones. Annotation revealed 7,058 protein-coding genes, with an average gene density of 4.6 kb per gene. Analysis of chromosome collinearity with the A. thaliana genome identified conserved synteny blocks encompassing the whole of the B. rapa chromosome A3 and sections of four A. thaliana chromosomes. The frequency of tandem duplication of genes differed between the conserved genome segments in B. rapa and A. thaliana, indicating differential rates of occurrence/retention of such duplicate copies of genes. Analysis of 'ancestral karyotype' genome building blocks enabled the development of a hypothetical model for the derivation of the B. rapa chromosome A3.

Conclusions

We report the near-complete chromosome sequence from a dicotyledonous crop species. This provides an example of the complexity of genome evolution following polyploidy. The high degree of contiguity afforded by the clone-by-clone approach provides a benchmark for the performance of whole genome shotgun approaches presently being applied in B. rapa and other species with complex genomes.     

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross‐species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next‐generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross‐genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter‐ and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.     

First survey of the wheat chromosome 5A composition through a next generation sequencing approach

Wheat is one of the world's most important crops and is characterized by a large polyploid genome. One way to reduce genome complexity is to isolate single chromosomes using flow cytometry. Low coverage DNA sequencing can provide a snapshot of individual chromosomes, allowing a fast characterization of their main features and comparison with other genomes. We used massively parallel 454 pyrosequencing to obtain a 2x coverage of wheat chromosome 5A. The resulting sequence assembly was used to identify TEs, genes and miRNAs, as well as to infer a virtual gene order based on the synteny with other grass genomes. Repetitive elements account for more than 75% of the genome. Gene content was estimated considering non-redundant reads showing at least one match to ESTs or proteins. The results indicate that the coding fraction represents 1.08% and 1.3% of the short and long arm respectively, projecting the number of genes of the whole chromosome to approximately 5,000. 195 candidate miRNA precursors belonging to 16 miRNA families were identified. The 5A genes were used to search for syntenic relationships between grass genomes. The short arm is closely related to Brachypodium chromosome 4, sorghum chromosome 8 and rice chromosome 12; the long arm to regions of Brachypodium chromosomes 4 and 1, sorghum chromosomes 1 and 2 and rice chromosomes 9 and 3. From these similarities it was possible to infer the virtual gene order of 392 (5AS) and 1,480 (5AL) genes of chromosome 5A, which was compared to, and found to be largely congruent with the available physical map of this chromosome.     

Rice-barley synteny and its application to saturation mapping of the barley Rpg1 region.

In order to facilitate the map-based cloning of the barley stem rust resistance gene Rpg1, we have demonstrated a high degree of synteny at a micro level between the telomeric region of barley chromosome 1P and rice chromosome 6. We have also developed and applied a simple and efficient method for selecting useful probes from large insert genomic YAC and cosmid clones. The gene order within the most terminal 6.5 cM of barley chromosome 1P was compared with the most terminal 2.7 cM of rice chromosome 6. Nine rice probes, previously mapped in rice or isolated from YAC or cosmid clones from this region, were mapped in barley. All, except one, were in synteny with the rice gene order. The exception, probe Y617R, was duplicated in barley. One copy was located on a different chromosome and the other in a non-syntenic position on barley chromosome 1P. The barley probes from this region could not be mapped to rice, but two of them were inferred to be in a syntenic location based on their position on a rice YAC. This work demonstrates the utility of applying the results of genetic and physical mapping of the small genome cereal rice to map-based cloning of interesting genes from large genome relatives.