Bootstrapping phylogenies inferred from rearrangement data

ABSTRACT: Large-scale sequencing of genomes has enabled the inference of phylogenies based on the evolution of genomic architecture, under such events as rearrangements, duplications, and losses. Many evolutionary models and associated algorithms have been designed over the last few years and have found use in comparative genomics and phylogenetic inference. However, the assessment of phylogenies built from such data has not been properly addressed to date. The standard method used in sequence-based phylogenetic inference is the bootstrap, but it relies on a large number of homologous characters that can be resampled; yet in the case of rearrangements, the entire genome is a single character. Alternatives such as the jackknife suffer from the same problem, while likelihood tests cannot be applied in the absence of well established probabilistic models. We present a new approach to the assessment of distance-based phylogenetic inference from whole-genome data; our approach combines features of the jackknife and the bootstrap and remains nonparametric. For each feature of our method, we give an equivalent feature in the sequence-based framework; we also present the results of extensive experimental testing, in both sequence-based and genome-based frameworks. Through the feature-by-feature comparison and the experimental results, we show that our bootstrapping approach is on par with the classic phylogenetic bootstrap used in sequence-based reconstruction, and we establish the clear superiority of the classic bootstrap for sequence data and of our corresponding new approach for rearrangement data over proposed variants. Finally, we test our approach on a small dataset of mammalian genomes, verifying that the support values match current thinking about the respective branches. Our method is the first to provide a standard of assessment to match that of the classic phylogenetic bootstrap for aligned sequences. Its support values follow a similar scale and its receiver-operating characteristics are nearly identical, indicating that it provides similar levels of sensitivity and specificity. Thus our assessment method makes it possible to conduct phylogenetic analyses on whole genomes with the same degree of confidence as for analyses on aligned sequences. Extensions to search-based inference methods such as maximum parsimony and maximum likelihood are possible, but remain to be thoroughly tested.     

Chromosome number and sex determination coevolve in turtles

Although much progress has been achieved in understanding the genetic basis of adaptation, the drivers of genome evolution remain obscure. For instance, extensive variation among reptilian genomes continues largely unexplained, yet reptiles hold critical clues about vertebrate evolution. Turtles posses diverse chromosome numbers (2N = 28-66) derived from extensive genomic rearrangements, plus varied sex-determining mechanisms (genotypic and temperature-dependent). Here, we show that rates of evolution in turtle chromosome number are ~20-fold higher along phylogenetic branches where transitions between sex-determining mechanisms also occur, revealing a strong coevolution of these traits and making drift a less likely driver. Directional tests indicate that both traits evolved effectively in synchrony. These events occurred near global extremes in temperature shifts over the last 200 million years, although the role of climate change remains unknown at this point. Two alternative testable explanations for these patterns are proposed. First, selection for sex determination turnover may co-opt mechanisms (e.g., chromatin remodeling) favoring genomic rearrangements. Alternatively, chromosomal rearrangements underlying diploid number evolution may alter gene regulation enabling transitions in sex-determining mechanisms. Our data indicate that the evolution of sex determination is intimately linked to profound genomic changes underlying diploid number evolution, the ecological context of which remains intriguing.     

Estimation of the Extent of Synteny Between Tetraodon nigroviridis and Homo sapiens Genomes

This paper presents a genomic comparison between 20 sequenced BACs (or fragments of BACs) from Tetraodon nigroviridis and the human genome. A total of 199 fish genes were identified by informatics resources, together with their putative human orthologues. Comparisons of the localizations in both species led to the identification of 32 syntenic regions and a minimum of 131 rearrangements in these regions that occurred during independent evolution of these species. This made it possible to estimate the rate of genomic rearrangements that occurred per million years (and per megabase). This rate is comparable to that obtained by comparison of the Fugu rubripes shotgun sequence data to human data but is significantly higher that those obtained by comparing the human genome to mammalian genomes. Overall, it suggests that genomic evolution by rearrangement is not uniform within the vertebrate group.Sequence data for the genomic BAC clones have been deposited with the DDBJ/EMBL/GenBank Data Libraries under accession numbers BX629360, BX629354, BX629355, BX629356, BX629357, BX629358, BX629359, and BX629360.     

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.     

Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish

ABSTRACT: BACKGROUND: The availability of a large number of recently sequenced vertebrate genomes opens new avenues to integrate cytogenetics and genomics in comparative and evolutionary studies. Cytogenetic mapping can offer alternative means to identify conserved synteny shared by distinct genomes and also to define genome regions that are still not fine characterized even after wide-ranging nucleotide sequence efforts. An efficient way to perform comparative cytogenetic mapping is based on BAC clones mapping by fluorescence in situ hybridization. In this report, to address the knowledge gap on the genome evolution in cichlid fishes, BAC clones of an Oreochromis niloticus library covering the linkage groups (LG) 1, 3, 5, and 7 were mapped onto the chromosomes of 9 African cichlid species. The cytogenetic mapping data were also integrated with BAC-end sequences information of O. niloticus and comparatively analyzed against the genome of other fish species and vertebrates. RESULTS: The location of BACs from LG1, 3, 5, and 7 revealed a strong chromosomal conservation among the analyzed cichlid species genomes, which evidenced a synteny of the markers of each LG. Comparative in silico analysis also identified large genomic blocks that were conserved in distantly related fish groups and also in other vertebrates. CONCLUSIONS: Although it has been suggested that fishes contain plastic genomes with high rates of chromosomal rearrangements and probably low rates of synteny conservation, our results evidence that large syntenic chromosome segments have been maintained conserved during evolution, at least for the considered markers. Additionally, our current cytogenetic mapping efforts integrated with genomic approaches conduct to a new perspective to address important questions involving chromosome evolution in fishes.     

Comparative organization of wheat homoeologous group 3S and 7L using wheat-rice synteny and identification of potential markers for genes controlling xanthophyll content in wheat

EST and genomic DNA sequencing efforts for rice and wheat have provided the basis for interpreting genome organization and evolution. In this study we have used EST and genomic sequencing information and a bioinformatic approach in a two-step strategy to align portions of the wheat and rice genomes. In the first step, wheat ESTs were used to identify rice orthologs and it was shown that wheat 3S and rice 1 contain syntenic units with intrachromosomal rearrangements. Further analysis using anchored rice contiguous sequences and TBLASTX alignments in a second alignment step showed interruptions by orthologous genes that map elsewhere in the wheat genome. This indicates that gene content and order is not as conserved as large chromosomal blocks as previously predicted. Similarly, chromosome 7L contains syntenic units with rice 6 and 8 but is interrupted by combinations of intrachromosomal and interchromosomal rearrangements involving syntenic units and single gene orthologs from other rice chromosome groups. We have used the rice sequence annotations to identify genes that can be used to develop markers linked to biosynthetic pathways on 3BS controlling xanthophyll production in wheat and thus involved in determining flour colour.Electronic Supplementary Material Supplementary material is available in the online version of this article at .     

Phylogenomics of the Maverick Virus-Like Mobile Genetic Elements of Vertebrates

Mavericks are virus-like mobile genetic elements found in the genomes of eukaryotes. Although Mavericks encode capsid morphogenesis homologs, their viral particles have not been observed. Here, we provide new evidence supporting the viral nature of Mavericks and the potential existence of virions. To this end, we conducted a phylogenomic analysis of Mavericks in hundreds of vertebrate genomes, discovering 134 elements with an intact coding capacity in 17 host species. We reveal an extensive genomic fossil record in 143 species and date three groups of elements to the Late Cretaceous. Bayesian phylogenetic analysis using genomic fossil orthologs suggests that Mavericks have infected osteichthyans for ∼419 My. They have undergone frequent cross-species transmissions in cyprinid fish and all core genes are subject to strong purifying selection. We conclude that vertebrate Mavericks form an ancient lineage of aquatic dsDNA viruses which are probably still functional in some vertebrate lineages.     

Molecular cytogenetic and genomic insights into chromosomal evolution

This review summarizes aspects of the extensive literature on the patterns and processes underpinning chromosomal evolution in vertebrates and especially placental mammals. It highlights the growing synergy between molecular cytogenetics and comparative genomics, particularly with respect to fully or partially sequenced genomes, and provides novel insights into changes in chromosome number and structure across deep division of the vertebrate tree of life. The examination of basal numbers in the deeper branches of the vertebrate tree suggest a haploid (n) chromosome number of 10-13 in an ancestral vertebrate, with modest increases in tetrapods and amniotes most probably by chromosomal fissioning. Information drawn largely from cross-species chromosome painting in the data-dense Placentalia permits the confident reconstruction of an ancestral karyotype comprising n=23 chromosomes that is similarly retained in Boreoeutheria. Using in silico genome-wide scans that include the newly released frog genome we show that of the nine ancient syntenies detected in conserved karyotypes of extant placentals (thought likely to reflect the structure of ancestral chromosomes), the human syntenic segmental associations 3p/21, 4pq/8p, 7a/16p, 14/15, 12qt/22q and 12pq/22qt predate the divergence of tetrapods. These findings underscore the enhanced quality of ancestral reconstructions based on the integrative molecular cytogenetic and comparative genomic approaches that collectively highlight a pattern of conserved syntenic associations that extends back ～360 million years ago.     

A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

Synorth: exploring the evolution of synteny and long-range regulatory interactions in vertebrate genomes

Genomic regulatory blocks are chromosomal regions spanned by long clusters of highly conserved noncoding elements devoted to long-range regulation of developmental genes, often immobilizing other, unrelated genes into long-lasting syntenic arrangements. Synorth is a web resource for exploring and categorizing the syntenic relationships in genomic regulatory blocks across multiple genomes, tracing their evolutionary fate after teleost whole genome duplication at the level of genomic regulatory block loci, individual genes, and their phylogenetic context.     

New approaches for reconstructing phylogenies from gene order data

We report on new techniques we have developed for reconstructing phylogenies on whole genomes. Our mathematical techniques include new polynomial-time methods for bounding the inversion length of a candidate tree and new polynomial-time methods for estimating genomic distances which greatly improve the accuracy of neighbor-joining analyses. We demonstrate the power of these techniques through an extensive performance study based on simulating genome evolution under a wide range of model conditions. Combining these new tools with standard approaches (fast reconstruction with neighbor-joining, exploration of all possible refinements of strict consensus trees, etc.) has allowed us to analyze datasets that were previously considered computationally impractical. In particular, we have conducted a complete phylogenetic analysis of a subset of the Campanulaceae family, confirming various conjectures about the relationships among members of the subset and about the principal mechanism of evolution for their chloroplast genome. We give representative results of the extensive experimentation we conducted on both real and simulated datasets in order to validate and characterize our approaches. We find that our techniques provide very accurate reconstructions of the true tree topology even when the data are generated by processes that include a significant fraction of transpositions and when the data are close to saturation.     

Details of Retropositional Genome Dynamics That Provide a Rationale for a Generic Division: The Distinct Branching of All the Pacific Salmon and Trout (Oncorhynchus) from the Atlantic Salmon and Trout (Salmo)

Salmonid species contain numerous short interspersed repetitive elements (SINEs), known collectively as the HpaI family, in their genomes. Amplification and successive integration of individual SINEs into the genomes have occurred during the evolution of salmonids. We reported previously a strategy for determining the phylogenetic relationships among the Pacific salmonids in which these SINEs were used as temporal landmarks of evolution. Here, we provide evidence for extensive genomic rearrangements that involved retropositions and deletions in a common ancestor of all the Pacific salmon and trout. Our results provide genetic support for the recent phylogenetic reassignment of steelhead and related species from the genus Salmo to the genus Oncorhynchus. Several other informative loci identified by insertions of HpaI SINEs have been isolated, and previously proposed branching orders of the Oncorhynchus species have been confirmed. The authenticity of our phylogenetic tree is supported both by the isolation of more than two informative loci per branching point and by the congruence of all our data, which suggest that the period between succesive speciations was sufficiently long for each SINE that had been amplified in the original species to become fixed in all individuals of that species.     

Comparison of Pax1/9 locus reveals 500-Myr-old syntenic block and evolutionary conserved noncoding regions

Identification of conserved genomic regions within and between different genomes is crucial when studying genome evolution. Here, we described regions of strong synteny conservation between vertebrate deuterostomes (tetrapods and teleosts) and invertebrate deuterostomes (amphioxus and sea urchin). The shared gene contents across phylogenetically distant species demonstrate that the conservation of the regions stemmed from an ancestral segment instead of a series of independent convergent events. Comparison of the syntenic regions allows us to postulate the primitive gene organization in the last common ancestor of deuterostomes and the evolutionary events that occurred to the 3 distinct lineages of sea urchin, amphioxus, and vertebrates after their separation. In addition, alignment of the syntenic regions led to the identification of 8 noncoding evolutionarily conserved regions shared between amphioxus and vertebrates. To our knowledge, this is the first report of conserved noncoding sequences shared by vertebrates and nonvertebrates. These noncoding sequences have high possibility of being elements that regulate neighboring genes. They are likely to be a factor in the maintenance of conserved synteny over long phylogenetic distance in different deuterostome lineages.     

Sequencing wheat chromosome arm 7BS delimits the 7BS/4AL translocation and reveals homoeologous gene conservation

Complex Triticeae genomes pose a challenge to genome sequencing efforts due to their size and repetitive nature. Genome sequencing can reveal details of conservation and rearrangements between related genomes. We have applied Illumina second generation sequencing technology to sequence and assemble the low copy and unique regions of Triticum aestivum chromosome arm 7BS, followed by the construction of a syntenic build based on gene order in Brachypodium. We have delimited the position of a previously reported translocation between 7BS and 4AL with a resolution of one or a few genes and report approximately 13% genes from 7BS having been translocated to 4AL. An additional 13 genes are found on 7BS which appear to have originated from 4AL. The gene content of the 7DS and 7BS syntenic builds indicate a total of ~77,000 genes in wheat. Within wheat syntenic regions, 7BS and 7DS share 740 genes and a common gene conservation rate of ~39% of the genes from the corresponding regions in Brachypodium, as well as a common rate of colinearity with Brachypodium of ~60%. Comparison of wheat homoeologues revealed ~84% of genes previously identified in 7DS have a homoeologue on 7BS or 4AL. The conservation rates we have identified among wheat homoeologues and with Brachypodium provide a benchmark of homoeologous gene conservation levels for future comparative genomic analysis. The syntenic build of 7BS is publicly available at http://www.wheatgenome.info.     

Rearrangement rate following the whole-genome duplication in teleosts

It is now clear that a whole-genome duplication (WGD) occurred at the base of the teleost fish lineage. Like the other anciently polyploid genomes investigated so far, teleost genomes now behave like diploids with chromosomes forming pairs at meiosis. The diploidization process is currently poorly understood. It is associated with many gene deletions, such that one of the duplicates is lost at most loci and has also been proposed to coincide with an increase in genomic instability. Here we ask whether WGD is a determinant of the genomic rearrangement rate in teleosts. We study variability of the rates of rearrangement along a vertebrate phylogenetic tree, composed of 3 tetrapods (human, chicken, and mouse) and 3 teleost fishes (zebrafish, Tetraodon, and Takifugu), whose complete genome sequences are available. We devise a simple parsimony method for counting rearrangements, which takes into account various methodological complications caused by the WGD and the subsequent gene losses. We show that there does appear to be an increase in rearrangement rate after WGD, but that there is also a great deal of additional variability in rearrangement rates across species.     

QuartetS: a fast and accurate algorithm for large-scale orthology detection

The unparalleled growth in the availability of genomic data offers both a challenge to develop orthology detection methods that are simultaneously accurate and high throughput and an opportunity to improve orthology detection by leveraging evolutionary evidence in the accumulated sequenced genomes. Here, we report a novel orthology detection method, termed QuartetS, that exploits evolutionary evidence in a computationally efficient manner. Based on the well-established evolutionary concept that gene duplication events can be used to discriminate homologous genes, QuartetS uses an approximate phylogenetic analysis of quartet gene trees to infer the occurrence of duplication events and discriminate paralogous from orthologous genes. We used function- and phylogeny-based metrics to perform a large-scale, systematic comparison of the orthology predictions of QuartetS with those of four other methods [bi-directional best hit (BBH), outgroup, OMA and QuartetS-C (QuartetS followed by clustering)], involving 624 bacterial genomes and >2 million genes. We found that QuartetS slightly, but consistently, outperformed the highly specific OMA method and that, while consuming only 0.5% additional computational time, QuartetS predicted 50% more orthologs with a 50% lower false positive rate than the widely used BBH method. We conclude that, for large-scale phylogenetic and functional analysis, QuartetS and QuartetS-C should be preferred, respectively, in applications where high accuracy and high throughput are required.     

Large-scale analysis of microRNA evolution

ABSTRACT: BACKGROUND: In animals, microRNAs (miRNA) are important genetic regulators. Animal miRNAs appear to have expanded in conjunction with an escalation in complexity during early bilaterian evolution. Their small size and high-degree of similarity makes them challenging for phylogenetic approaches. Furthermore, genomic locations encoding miRNAs are not clearly defined in many species. A number of studies have looked at the evolution of individual miRNA families. However, we currently lack resources for large-scale analysis of miRNA evolution. RESULTS: We addressed some of these issues in order to analyse the evolution of miRNAs. We perform syntenic and phylogenetic analysis for miRNAs from 80 animal species. We present synteny maps, phylogenies and functional data for miRNAs across these species. These data represent the basis of our analyses and also act as a resource for the community. CONCLUSIONS: We use these data to explore the distribution of miRNAs across phylogenetic space, characterise their birth and death, and examine functional relationships between miRNAs and other genes. These data confirm a number of previously reported findings on a larger scale and also offer novel insights into the evolution of the miRNA repertoire in animals, and it's genomic organization.     

Comparative genomic organization of the human and bovine PRNP locus

We sequenced a 208-kb BAC clone spanning the bovine prion protein (PRNP) locus, and compared the genomic structure with that of human. As a result, we determined the precise breakpoint between the two syntenic genomes, located on the 5' UTR of the PRNP gene, and discovered two highly repetitive sequences near the breakpoint. Further analysis demonstrated that the genomic structure of three genes, PRNP, PRND, and RASSF2, within the syntenic region of the bovine genome is highly conserved in order and orientation. The PRNT locus was not found in bovine but is conserved in several primates, including human. Moreover, we confirmed that the bovine RASSF2 is composed of 10 exons, as is the human gene, showing some difference from a previous report. Our findings may provide useful clues for understanding the evolutional process in the PRNP locus and also the mechanism that allows TSE from cattle to infect humans.     

Unraveling the effect of genomic structural changes in the rhesus macaque - implications for the adaptive role of inversions

Background

By reshuffling genomes, structural genomic reorganizations provide genetic variation on which natural selection can work. Understanding the mechanisms underlying this process has been a long-standing question in evolutionary biology. In this context, our purpose in this study is to characterize the genomic regions involved in structural rearrangements between human and macaque genomes and determine their influence on meiotic recombination as a way to explore the adaptive role of genome shuffling in mammalian evolution.

Results

We first constructed a highly refined map of the structural rearrangements and evolutionary breakpoint regions in the human and rhesus macaque genomes based on orthologous genes and whole-genome sequence alignments. Using two different algorithms, we refined the genomic position of known rearrangements previously reported by cytogenetic approaches and described new putative micro-rearrangements (inversions and indels) in both genomes. A detailed analysis of the rhesus macaque genome showed that evolutionary breakpoints are in gene-rich regions, being enriched in GO terms related to immune system. We also identified defense-response genes within a chromosome inversion fixed in the macaque lineage, underlying the relevance of structural genomic changes in evolutionary and/or adaptation processes. Moreover, by combining in silico and experimental approaches, we studied the recombination pattern of specific chromosomes that have suffered rearrangements between human and macaque lineages.

Conclusions

Our data suggest that adaptive alleles – in this case, genes involved in the immune response – might have been favored by genome rearrangements in the macaque lineage.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-530) contains supplementary material, which is available to authorized users.