Dataset Information

Performance and scalability of discriminative metrics for comparative gene identification in 12 Drosophila genomes.

ABSTRACT: Comparative genomics of multiple related species is a powerful methodology for the discovery of functional genomic elements, and its power should increase with the number of species compared. Here, we use 12 Drosophila genomes to study the power of comparative genomics metrics to distinguish between protein-coding and non-coding regions. First, we study the relative power of different comparative metrics and their relationship to single-species metrics. We find that even relatively simple multi-species metrics robustly outperform advanced single-species metrics, especially for shorter exons (< or =240 nt), which are common in animal genomes. Moreover, the two capture largely independent features of protein-coding genes, with different sensitivity/specificity trade-offs, such that their combinations lead to even greater discriminatory power. In addition, we study how discovery power scales with the number and phylogenetic distance of the genomes compared. We find that species at a broad range of distances are comparably effective informants for pairwise comparative gene identification, but that these are surpassed by multi-species comparisons at similar evolutionary divergence. In particular, while pairwise discovery power plateaued at larger distances and never outperformed the most advanced single-species metrics, multi-species comparisons continued to benefit even from the most distant species with no apparent saturation. Last, we find that genes in functional categories typically considered fast-evolving can nonetheless be recovered at very high rates using comparative methods. Our results have implications for comparative genomics analyses in any species, including the human.

SUBMITTER: Lin MF

PROVIDER: S-EPMC2291194 | biostudies-literature | 2008 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Performance and scalability of discriminative metrics for comparative gene identification in 12 Drosophila genomes.

Lin Michael F MF Deoras Ameya N AN Rasmussen Matthew D MD Kellis Manolis M

PLoS computational biology 20080418 4

Comparative genomics of multiple related species is a powerful methodology for the discovery of functional genomic elements, and its power should increase with the number of species compared. Here, we use 12 Drosophila genomes to study the power of comparative genomics metrics to distinguish between protein-coding and non-coding regions. First, we study the relative power of different comparative metrics and their relationship to single-species metrics. We find that even relatively simple multi- ...[more]

PMID: 18421375

Similar Datasets

Project description:BackgroundHarnessing vast amounts of genomic data in phylogenetic context stemming from massive sequencing of multiple closely related genomes requires new tools and approaches. We present a tool for the genome-wide analysis of frequencies and patterns of amino acid substitutions in multiple alignments of genes' coding regions, and a database of amino acid substitutions in the phylogeny of 12 Drosophila genomes. We illustrate the use of these resources to address three types of evolutionary genomics questions: about fluxes in amino acid composition in proteins, about asymmetries in amino acid substitutions and about patterns of molecular evolution in duplicated genes.ResultsWe demonstrate that amino acid composition of Drosophila proteins underwent a significant shift over the last 70 million years encompassed by the studied phylogeny, with less common amino acids (Cys, Met, His) increasing in frequency and more common ones (Ala, Leu, Glu) becoming less frequent. These fluxes are strongly correlated with polarity of source and destination amino acids, resulting in overall systematic decrease of mean polarity of amino acids found in Drosophila proteins. Frequency and radicality of amino acid substitutions are higher in paralogs than in orthologous single-copy genes and are higher in gene families with paralogs than in gene families without surviving duplications. Rate and radicality of substitutions, as expected, are negatively correlated with overall level and uniformity of gene expression. However, these correlations are not observed for substitutions occurring in duplicated genes, indicating a different selective constraint on the evolution of paralogous sequences. Clades resulting from duplications show a marked asymmetry in rate and radicality of amino acid substitutions, possibly a signal of widespread neofunctionalization. These patterns differ among protein families of different functionality, with genes coding for RNA-binding proteins differing from most other functional groups in terms of amino acid substitution patterns in duplicated and single-copy genes.ConclusionsWe demonstrate that deep phylogenetic analysis of amino acid substitutions can reveal interesting genome-wide patterns. Amino acid composition of drosophilid proteins is shaped by fluxes similar to those previously observed in prokaryotic, yeast and mammalian genomes, indicating globally present patterns. Increased frequency and radicality of amino acid substitutions in duplicated genes and the presence of asymmetry of these parameters between paralogous clades indicate widespread neofunctionalization among paralogs as the mechanism of duplication retention.

Project description:BackgroundMiniature inverted-repeat transposable elements (MITEs) are non-autonomous DNA-mediated transposable elements (TEs) derived from autonomous TEs. Unlike in many plants or animals, MITEs and other types of DNA-mediated TEs were previously thought to be either rare or absent in Drosophila. Most other TE families in Drosophila exist at low or intermediate copy number (around < 100 per genome).ResultsWe present evidence here that the dispersed repeat Drosophila interspersed element 1 (DINE-1; also named INE-1 and DNAREP1) is a highly abundant DNA-mediated TE containing inverted repeats found in all 12 sequenced Drosophila genomes. All DINE-1s share a similar sequence structure, but are more homogeneous within species than they are among species. The inferred phylogenetic relationship of the DINE-1 consensus sequence from each species is generally consistent with the known species phylogeny, suggesting vertical transmission as the major mechanism for DINE-1 propagation. Exceptions observed in D. willistoni and D. ananassae could be due to either horizontal transfer or reactivation of ancestral copies. Our analysis of pairwise percentage identity of DINE-1 copies within species suggests that the transpositional activity of DINE-1 is extremely dynamic, with some lineages showing evidence for recent transpositional bursts and other lineages appearing to have silenced their DINE-1s for long periods of time. We also find that all species have many DINE-1 insertions in introns and adjacent to protein-coding genes. Finally, we discuss our results in light of a recent proposal that DINE-1s belong to the Helitron family of TEs.ConclusionWe find that all 12 Drosophila species with whole-genome sequence contain the high copy element DINE-1. Although all DINE-1s share a similar structure, species-specific variation in the distribution of average pairwise divergence suggests that DINE-1 has gone through multiple independent cycles of activation and suppression. DINE-1 also has had a significant impact on gene structure evolution.

Dataset Information

Performance and scalability of discriminative metrics for comparative gene identification in 12 Drosophila genomes.

Publications

Performance and scalability of discriminative metrics for comparative gene identification in 12 Drosophila genomes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets