Comparative Genomics Tools Research Articles

Evolutionary Dynamics of Chloroplast Genomes in Low Light: A Case Study of the Endolithic Green Alga Ostreobium quekettii.

Some photosynthetic organisms live in extremely low light environments. Light limitation is associated with selective forces as well as reduced exposure to mutagens, and over evolutionary timescales it can leave a footprint on species’ genomes. Here, we present the chloroplast genomes of four green algae (Bryopsidales, Ulvophyceae), including the endolithic (limestone-boring) alga Ostreobium quekettii, which is a low light specialist. We use phylogenetic models and comparative genomic tools to investigate whether the chloroplast genome of Ostreobium corresponds to our expectations of how low light would affect genome evolution. Ostreobium has the smallest and most gene-dense chloroplast genome among Ulvophyceae reported to date, matching our expectation that light limitation would impose resource constraints reflected in the chloroplast genome architecture. Rates of molecular evolution are significantly slower along the phylogenetic branch leading to Ostreobium, in agreement with the expected effects of low light and energy levels on molecular evolution. We expected the ability of Ostreobium to perform photosynthesis in very low light to be associated with positive selection in genes related to the photosynthetic machinery, but instead, we observed that these genes may be under stronger purifying selection. Besides shedding light on the genome dynamics associated with a low light lifestyle, this study helps to resolve the role of environmental factors in shaping the diversity of genome architectures observed in nature.

Aligning the unalignable: bacteriophage whole genome alignments.

BackgroundIn recent years, many studies focused on the description and comparison of large sets of related bacteriophage genomes. Due to the peculiar mosaic structure of these genomes, few informative approaches for comparing whole genomes exist: dot plots diagrams give a mostly qualitative assessment of the similarity/dissimilarity between two or more genomes, and clustering techniques are used to classify genomes. Multiple alignments are conspicuously absent from this scene. Indeed, whole genome aligners interpret lack of similarity between sequences as an indication of rearrangements, insertions, or losses. This behavior makes them ill-prepared to align bacteriophage genomes, where even closely related strains can accomplish the same biological function with highly dissimilar sequences.ResultsIn this paper, we propose a multiple alignment strategy that exploits functional collinearity shared by related strains of bacteriophages, and uses partial orders to capture mosaicism of sets of genomes. As classical alignments do, the computed alignments can be used to predict that genes have the same biological function, even in the absence of detectable similarity. The Alpha aligner implements these ideas in visual interactive displays, and is used to compute several examples of alignments of Staphylococcus aureus and Mycobacterium bacteriophages, involving up to 29 genomes. Using these datasets, we prove that Alpha alignments are at least as good as those computed by standard aligners. Comparison with the progressiveMauve aligner – which implements a partial order strategy, but whose alignments are linearized – shows a greatly improved interactive graphic display, while avoiding misalignments.ConclusionsMultiple alignments of whole bacteriophage genomes work, and will become an important conceptual and visual tool in comparative genomics of sets of related strains.A python implementation of Alpha, along with installation instructions for Ubuntu and OSX, is available on bitbucket (https://bitbucket.org/thekswenson/alpha).Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-015-0869-5) contains supplementary material, which is available to authorized users.

Bioinformatics Tools for Exploring the SUMO Gene Network.

Plant sumoylation research has seen significant advances in recent years, particularly since high-throughput proteomics strategies have enabled the discovery of hundreds of potential SUMO targets and interactors of SUMO pathway components. In the present chapter, we introduce the SUMO Gene Network (SGN), a curated assembly of Arabidopsis thaliana genes that have been functionally associated with sumoylation, from SUMO pathway components to targets and interactors. The enclosed tutorial helps interpret and manage these datasets, and details bioinformatics tools that can be used for in silico-based hypothesis generation. The latter include tools for sumoylation site prediction, comparative genomics, and gene network analysis.

Genome-wide characterization and expression analysis of common bean bHLH transcription factors in response to excess salt concentration.

Members of basic helix-loop-helix (bHLH) gene family found in all eukaryotes play crucial roles in response to stress. Though, most eukaryotes carry the proteins of this family, biological functions of the most bHLH family members are not deeply evaluated in plants. In this study, we conducted a comprehensive genome-wide analysis of bHLH transcription factors in salt tolerant common bean. We identified 155 bHLH protein-encoding genes (PvbHLH) by using in silico comparative genomics tools. Based on the phylogenetic tree, PvbHLH genes were classified into 8 main groups with 21 subfamilies. Exon-intron analysis indicated that proteins belonging to same main groups exhibited a closely related gene structure. While, the PvbHLH gene family has been mainly expanded through segmental duplications, a total of 11 tandem duplication were detected. Genome-wide expression analysis of bHLH genes showed that 63 PvbHLH genes were differentially expressed in at least one tissue. Three of them displayed higher expression values in both leaf and root tissues. The in silico micro-RNA target transcript analyses revealed that totally 100 PvHLH genes targeted by 86 plant miRNAs. The most abundant transcripts, which were targeted by all 18 plant miRNA, were belonging to PvHLH-22 and PvHLH-44 genes. The expression of 16 PvbHLH genes in the root and leaf tissues of salt-stressed common bean was evaluated using qRT-PCR. Among them, two of PvbHLHs, PvbHLH-54, PvbHLH-148, were found to be up-regulated in both tissues in correlation with RNA-seq measurements. The results of this study could help improve understanding of biological functions of common bean bHLH family under salt stress. Additionally, it may provide basic resources for analyzing bHLH protein function for improving economic, agronomic and ecological benefit in common bean and other species.

High resolution physical mapping of single gene fragments on pachytene chromosome 4 and 7 of Rosa.

BackgroundRosaceae is a family containing many economically important fruit and ornamental species. Although fluorescence in situ hybridization (FISH)-based physical mapping of plant genomes is a valuable tool for map-based cloning, comparative genomics and evolutionary studies, no studies using high resolution physical mapping have been performed in this family. Previously we proved that physical mapping of single-copy genes as small as 1.1 kb is possible on mitotic metaphase chromosomes of Rosa wichurana using Tyramide-FISH. In this study we aimed to further improve the physical map of Rosa wichurana by applying high resolution FISH to pachytene chromosomes.ResultsUsing high resolution Tyramide-FISH and multicolor Tyramide-FISH, 7 genes (1.7–3 kb) were successfully mapped on pachytene chromosomes 4 and 7 of Rosa wichurana. Additionally, by using multicolor Tyramide-FISH three closely located genes were simultaneously visualized on chromosome 7. A detailed map of heterochromatine/euchromatine patterns of chromosome 4 and 7 was developed with indication of the physical position of these 7 genes. Comparison of the gene order between Rosa wichurana and Fragaria vesca revealed a poor collinearity for chromosome 7, but a perfect collinearity for chromosome 4.ConclusionsHigh resolution physical mapping of short probes on pachytene chromosomes of Rosa wichurana was successfully performed for the first time. Application of Tyramide-FISH on pachytene chromosomes allowed the mapping resolution to be increased up to 20 times compared to mitotic metaphase chromosomes. High resolution Tyramide-FISH and multicolor Tyramide-FISH might become useful tools for further physical mapping of single-copy genes and for the integration of physical and genetic maps of Rosa wichurana and other members of the Rosaceae.Electronic supplementary materialThe online version of this article (doi:10.1186/s12863-015-0233-9) contains supplementary material, which is available to authorized users.

A quantitative hybridization approach using 17 DNA markers for identification and clustering analysis of Ralstonia solanacearum

Ralstonia solanacearum (Rs) is a quarantine phytopathogenic bacterium accountable for heavy economic losses worldwide. Monitoring and eradication programmes required for this pathogen are dependent on the availability of time‐ and cost‐efficient detection and typing methods. However, members of the Rs species complex are characterized by a high phenotypic and genetic diversity, which requires improved diagnostics methods. The currently available full genome sequences of several Rs strains allow for the selection of novel specific DNA markers using comparative genomics tools. In this work, 17 novel markers were selected based on Rs‐specific protein domains and thoroughly validated for specificity and stability, both in silico and using ‘wet lab’ assays. Polymerase chain reaction‐ and hybridization‐based validation assays revealed that the DNA regions selected as markers were unevenly distributed amongst the tested strains, with nine markers present throughout the species complex. The distribution of the remaining eight markers was highly variable between the different analysed strains and enabled the attainment of strain‐specific dot blot hybridization patterns, particularly informative for typing. The average probability value of each strain being positive for each of the 17 markers was calculated by an algorithm and used to obtain a dendrogram representing hierarchical clustering analysis of Rs, according to the similarity of their hybridization patterns. This method should prove to be a robust and straightforward procedure for genotyping members of the Rs species complex. Furthermore, this quantitative hybridization approach will allow the construction of informative databases to determine new Rs genotypes and infer epidemiological patterns.

The Alternaria genomes database: a comprehensive resource for a fungal genus comprised of saprophytes, plant pathogens, and allergenic species.

BackgroundAlternaria is considered one of the most common saprophytic fungal genera on the planet. It is comprised of many species that exhibit a necrotrophic phytopathogenic lifestyle. Several species are clinically associated with allergic respiratory disorders although rarely found to cause invasive infections in humans. Finally, Alternaria spp. are among the most well known producers of diverse fungal secondary metabolites, especially toxins.DescriptionWe have recently sequenced and annotated the genomes of 25 Alternaria spp. including but not limited to many necrotrophic plant pathogens such as A. brassicicola (a pathogen of Brassicaceous crops like cabbage and canola) and A. solani (a major pathogen of Solanaceous plants like potato and tomato), and several saprophytes that cause allergy in human such as A. alternata isolates. These genomes were annotated and compared. Multiple genetic differences were found in the context of plant and human pathogenicity, notably the pro-inflammatory potential of A. alternata. The Alternaria genomes database was built to provide a public platform to access the whole genome sequences, genome annotations, and comparative genomics data of these species. Genome annotation and comparison were performed using a pipeline that integrated multiple computational and comparative genomics tools. Alternaria genome sequences together with their annotation and comparison data were ported to Ensembl database schemas using a self-developed tool (EnsImport). Collectively, data are currently hosted using a customized installation of the Ensembl genome browser platform.ConclusionRecent efforts in fungal genome sequencing have facilitated the studies of the molecular basis of fungal pathogenicity as a whole system. The Alternaria genomes database provides a comprehensive resource of genomics and comparative data of an important saprophytic and plant/human pathogenic fungal genus. The database will be updated regularly with new genomes when they become available. The Alternaria genomes database is freely available for non-profit use at http://alternaria.vbi.vt.edu.

Unlimited Thirst for Genome Sequencing, Data Interpretation, and Database Usage in Genomic Era: The Road towards Fast-Track Crop Plant Improvement.

The number of sequenced crop genomes and associated genomic resources is growing rapidly with the advent of inexpensive next generation sequencing methods. Databases have become an integral part of all aspects of science research, including basic and applied plant and animal sciences. The importance of databases keeps increasing as the volume of datasets from direct and indirect genomics, as well as other omics approaches, keeps expanding in recent years. The databases and associated web portals provide at a minimum a uniform set of tools and automated analysis across a wide range of crop plant genomes. This paper reviews some basic terms and considerations in dealing with crop plant databases utilization in advancing genomic era. The utilization of databases for variation analysis with other comparative genomics tools, and data interpretation platforms are well described. The major focus of this review is to provide knowledge on platforms and databases for genome-based investigations of agriculturally important crop plants. The utilization of these databases in applied crop improvement program is still being achieved widely; otherwise, the end for sequencing is not far away.

GenomeVISTA--an integrated software package for whole-genome alignment and visualization.

With the ubiquitous generation of complete genome assemblies for a variety of species, efficient tools for whole-genome alignment along with user-friendly visualization are critically important. Our VISTA family of tools for comparative genomics, based on algorithms for pairwise and multiple alignments of genomic sequences and whole-genome assemblies, has become one of the standard techniques for comparative analysis. Most of the VISTA programs have been implemented as Web-accessible servers and are extensively used by the biomedical community. In this manuscript, we introduce GenomeVISTA: a novel implementation that incorporates most features of the VISTA family--fast and accurate alignment, visualization capabilities, GUI and analytical tools within a stand-alone software package. GenomeVISTA thus provides flexibility and security for users who need to conduct whole-genome comparisons on their own computers. Implemented in Perl, C/C++ and Java, the source code is freely available for download at the VISTA Web site: http://genome.lbl.gov/vista/.

CoryneBase: Corynebacterium Genomic Resources and Analysis Tools at Your Fingertips

Corynebacteria are used for a wide variety of industrial purposes but some species are associated with human diseases. With increasing number of corynebacterial genomes having been sequenced, comparative analysis of these strains may provide better understanding of their biology, phylogeny, virulence and taxonomy that may lead to the discoveries of beneficial industrial strains or contribute to better management of diseases. To facilitate the ongoing research of corynebacteria, a specialized central repository and analysis platform for the corynebacterial research community is needed to host the fast-growing amount of genomic data and facilitate the analysis of these data. Here we present CoryneBase, a genomic database for Corynebacterium with diverse functionality for the analysis of genomes aimed to provide: (1) annotated genome sequences of Corynebacterium where 165,918 coding sequences and 4,180 RNAs can be found in 27 species; (2) access to comprehensive Corynebacterium data through the use of advanced web technologies for interactive web interfaces; and (3) advanced bioinformatic analysis tools consisting of standard BLAST for homology search, VFDB BLAST for sequence homology search against the Virulence Factor Database (VFDB), Pairwise Genome Comparison (PGC) tool for comparative genomic analysis, and a newly designed Pathogenomics Profiling Tool (PathoProT) for comparative pathogenomic analysis. CoryneBase offers the access of a range of Corynebacterium genomic resources as well as analysis tools for comparative genomics and pathogenomics. It is publicly available at http://corynebacterium.um.edu.my/.

Comparative genomic analysis of the R2R3 MYB secondary cell wall regulators of Arabidopsis, poplar, rice, maize, and switchgrass

BackgroundR2R3 MYB proteins constitute one of the largest plant transcription factor clades and regulate diverse plant-specific processes. Several R2R3 MYB proteins act as regulators of secondary cell wall (SCW) biosynthesis in Arabidopsis thaliana (At), a dicotyledenous plant. Relatively few studies have examined SCW R2R3 MYB function in grasses, which may have diverged from dicots in terms of SCW regulatory mechanisms, as they have in cell wall composition and patterning. Understanding cell wall regulation is especially important for improving lignocellulosic bioenergy crops, such as switchgrass.ResultsHere, we describe the results of applying phylogenic, OrthoMCL, and sequence identity analyses to classify the R2R3 MYB family proteins from the annotated proteomes of Arabidposis, poplar, rice, maize and the initial genome (v0.0) and translated transcriptome of switchgrass (Panicum virgatum). We find that the R2R3 MYB proteins of the five species fall into 48 subgroups, including three dicot-specific, six grass-specific, and two panicoid grass-expanded subgroups. We observe four classes of phylogenetic relationships within the subgroups of known SCW-regulating MYB proteins between Arabidopsis and rice, ranging from likely one-to-one orthology (for AtMYB26, AtMYB103, AtMYB69) to no homologs identifiable (for AtMYB75). Microarray data for putative switchgrass SCW MYBs indicate that many maintain similar expression patterns with the Arabidopsis SCW regulators. However, some of the switchgrass-expanded candidate SCW MYBs exhibit differences in gene expression patterns among paralogs, consistent with subfunctionalization. Furthermore, some switchgrass representatives of grass-expanded clades have gene expression patterns consistent with regulating SCW development.ConclusionsOur analysis suggests that no single comparative genomics tool is able to provide a complete picture of the R2R3 MYB protein family without leaving ambiguities, and establishing likely false-negative and -positive relationships, but that used together a relatively clear view emerges. Generally, we find that most R2R3 MYBs that regulate SCW in Arabidopsis are likely conserved in the grasses. This comparative analysis of the R2R3 MYB family will facilitate transfer of understanding of regulatory mechanisms among species and enable control of SCW biosynthesis in switchgrass toward improving its biomass quality.

Hormonal modulators of glial ABCA1 and apoE levels

Apolipoprotein E (apoE) is the major lipid carrier in the central nervous system. As apoE plays a major role in the pathogenesis of Alzheimer disease (AD) and also mediates repair pathways after several forms of acute brain injury, modulating the expression, secretion, or function of apoE may provide potential therapeutic approaches for several neurological disorders. Here we show that progesterone and a synthetic progestin, lynestrenol, significantly induce apoE secretion from human CCF-STTG1 astrocytoma cells, whereas estrogens and the progesterone metabolite allopregnanolone have negligible effects. Intriguingly, lynestrenol also increases expression of the cholesterol transporter ABCA1 in CCF-STTG1 astrocytoma cells, primary murine glia, and immortalized murine astrocytes that express human apoE3. The progesterone receptor inhibitor RU486 attenuates the effect of progestins on apoE expression in CCF-STTG1 astrocytoma cells but has no effect on ABCA1 expression in all glial cell models tested, suggesting that the progesterone receptor (PR) may participate in apoE but does not affect ABCA1 regulation. These results suggest that selective reproductive steroid hormones have the potential to influence glial lipid homeostasis through liver X receptor-dependent and progesterone receptor-dependent pathways.

Phylogeography and population structure of the biologically invasive phytopathogen Erwinia amylovora inferred using minisatellites

Erwinia amylovora causes a major disease of pome fruit trees worldwide, and is regulated as a quarantine organism in many countries. While some diversity of isolates has been observed, molecular epidemiology of this bacterium is hindered by a lack of simple molecular typing techniques with sufficiently high resolution. We report a molecular typing system of E. amylovora based on variable number of tandem repeats (VNTR) analysis. Repeats in the E. amylovora genome were identified with comparative genomic tools, and VNTR markers were developed and validated. A Multiple-Locus VNTR Analysis (MLVA) was applied to E. amylovora isolates from bacterial collections representing global and regional distribution of the pathogen. Based on six repeats, MLVA allowed the distinction of 227 haplotypes among a collection of 833 isolates of worldwide origin. Three geographically separated groups were recognized among global isolates using Bayesian clustering methods. Analysis of regional outbreaks confirmed presence of diverse haplotypes but also high representation of certain haplotypes during outbreaks. MLVA analysis is a practical method for epidemiological studies of E. amylovora, identifying previously unresolved population structure within outbreaks. Knowledge of such structure can increase our understanding on how plant diseases emerge and spread over a given geographical region.

BlastGraph: a comparative genomics tool based on BLAST and graph algorithms

BlastGraph is an interactive Java program for comparative genome analysis based on Basic Local Alignment Search Tool (BLAST), graph clustering and data visualization. The software generates clusters of sequences of multiple genomes from all-to-all BLAST results and visualizes the results in graph plots together with related information such as sequence features, gene conservation and similarity relationships. Pruning algorithms are used to reduce results to more meaningful subclusters. Subsequent analyses can then be conducted based on the predicted clusters, including gene content, genome phylogenetics and gene gain and loss. https://github.com/bigwiv/BlastGraph.

A Whole Genome Pairwise Comparative and Functional Analysis of Geobacter sulfurreducens PCA

Geobacter species are involved in electricity production, bioremediations, and various environmental friendly activities. Whole genome comparative analyses of Geobacter sulfurreducens PCA, Geobacter bemidjiensis Bem, Geobacter sp. FRC-32, Geobacter lovleyi SZ, Geobacter sp. M21, Geobacter metallireducens GS-15, Geobacter uraniireducens Rf4 have been made to find out similarities and dissimilarities among them. For whole genome comparison of Geobacter species, an in-house tool, Geobacter Comparative Genomics Tool (GCGT) has been developed using BLASTALL program, and these whole genome analyses yielded conserved genes and they are used for functional prediction. The conserved genes identified are about 2184 genes, and these genes are classified into 14 groups based on the pathway information. Functions for 74 hypothetical proteins have been predicted based on the conserved genes. The predicted functions include pilus type proteins, flagellar proteins, ABC transporters, and other proteins which are involved in electron transfer. A phylogenetic tree from 16S rRNA of seven Geobacter species showed that G. sulfurreducens PCA is closely related to G. metallireducens GS-15 and G. lovleyi SZ. For evolutionary study, acetate kinase protein is used, which showed closeness to Pelobacter propionicus, Pelobacter carbinolicus, and Deferribacteraceae family bacterial species. These results will be useful to enhance electricity production by using biotechnological approaches.

A high-resolution cucumber cytogenetic map integrated with the genome assembly

BackgroundHigh-resolution cytogenetic map can provide not only important biological information on genome organization but also solid foundation for genetic and genomic research. The progress in the molecular and cytogenetic studies has created the basis for developing the cytogenetic map in cucumber (Cucumis sativus L.).ResultsHere, the cytogenetic maps of four cucumber chromosomes (chromosomes 1, 3–5) were constructed by fluorescence in situ hybridization (FISH) analysis on cucumber pachytene chromosomes. Together with our previously constructed cytogenetic maps of three cucumber chromosomes (chromosomes 2, 6–7), cucumber has a complete cytogenetic map with 76 anchoring points between the genetic, the cytogenetic and the draft genome assembly maps. To compare our pachytene FISH map directly to the genetic linkage and draft genome assembly maps, we used a standardized map unit—relative map position (RMP) to produce the comparative map alignments. The alignments allowed a global view of the relationship of genetic and physical distances along each cucumber chromosome, and accuracy and coverage of the draft genome assembly map.ConclusionsWe demonstrated a good correlation between positions of the markers in the linkage and physical maps, and essentially complete coverage of chromosome arms by the draft genome assembly. Our study not only provides essential information for the improvement of sequence assembly but also offers molecular tools for cucumber genomics research, comparative genomics and evolutionary study.

An automated graphics tool for comparative genomics: the Coulson plot generator

BackgroundComparative analysis is an essential component to biology. When applied to genomics for example, analysis may require comparisons between the predicted presence and absence of genes in a group of genomes under consideration. Frequently, genes can be grouped into small categories based on functional criteria, for example membership of a multimeric complex, participation in a metabolic or signaling pathway or shared sequence features and/or paralogy. These patterns of retention and loss are highly informative for the prediction of function, and hence possible biological context, and can provide great insights into the evolutionary history of cellular functions. However, representation of such information in a standard spreadsheet is a poor visual means from which to extract patterns within a dataset.ResultsWe devised the Coulson Plot, a new graphical representation that exploits a matrix of pie charts to display comparative genomics data. Each pie is used to describe a complex or process from a separate taxon, and is divided into sectors corresponding to the number of proteins (subunits) in a complex/process. The predicted presence or absence of proteins in each complex are delineated by occupancy of a given sector; this format is visually highly accessible and makes pattern recognition rapid and reliable. A key to the identity of each subunit, plus hierarchical naming of taxa and coloring are included. A java-based application, the Coulson plot generator (CPG) automates graphic production, with a tab or comma-delineated text file as input and generating an editable portable document format or svg file.ConclusionsCPG software may be used to rapidly convert spreadsheet data to a graphical matrix pie chart format. The representation essentially retains all of the information from the spreadsheet but presents a graphically rich format making comparisons and identification of patterns significantly clearer. While the Coulson plot format is highly useful in comparative genomics, its original purpose, the software can be used to visualize any dataset where entity occupancy is compared between different classes.AvailabilityCPG software is available at sourceforge http://sourceforge.net/projects/coulson and http://dl.dropbox.com/u/6701906/Web/Sites/Labsite/CPG.html

Toward systems metabolic engineering of Aspergillus and Pichia species for the production of chemicals and biofuels

Recently genome sequence data have become available for Aspergillus and Pichia species of industrial interest. This has stimulated the use of systems biology approaches for large-scale analysis of the molecular and metabolic responses of Aspergillus and Pichia under defined conditions, which has resulted in much new biological information. Case-specific contextualization of this information has been performed using comparative and functional genomic tools. Genomics data are also the basis for constructing genome-scale metabolic models, and these models have helped in the contextualization of knowledge on the fundamental biology of Aspergillus and Pichia species. Furthermore, with the availability of these models, the engineering of Aspergillus and Pichia is moving from traditional approaches, such as random mutagenesis, to a systems metabolic engineering approach. Here we review the recent trends in systems biology of Aspergillus and Pichia species, highlighting the relevance of these developments for systems metabolic engineering of these organisms for the production of hydrolytic enzymes, biofuels and chemicals from biomass.

A radiation hybrid map of chromosome 1D reveals synteny conservation at a wheat speciation locus

The species cytoplasm specific (scs) genes affect nuclear-cytoplasmic interactions in interspecific hybrids. A radiation hybrid (RH) mapping population of 188 individuals was employed to refine the location of the scs (ae) locus on Triticum aestivum chromosome 1D. "Wheat Zapper," a comparative genomics tool, was used to predict synteny between wheat chromosome 1D, Oryza sativa, Brachypodium distachyon, and Sorghum bicolor. A total of 57 markers were developed based on synteny or literature and genotyped to produce a RH map spanning 205.2cR. A test-cross methodology was devised for phenotyping of RH progenies, and through forward genetic, the scs (ae) locus was pinpointed to a 1.1Mb-segment containing eight genes. Further, the high resolution provided by RH mapping, combined with chromosome-wise synteny analysis, located the ancestral point of fusion between the telomeric and centromeric repeats of two paleochromosomes that originated chromosome 1D. Also, it indicated that the centromere of this chromosome is likely the result of a neocentromerization event, rather than the conservation of an ancestral centromere as previously believed. Interestingly, location of scs locus in the vicinity of paleofusion is not associated with the expected disruption of synteny, but rather with a good degree of conservation across grass species. Indeed, these observations advocate the evolutionary importance of this locus as suggested by "Maan's scs hypothesis."

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

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