Partial Homology Relations - Satisfiability in terms of Di-Cographs

Directed cographs (di-cographs) play a crucial role in the reconstruction of evolutionary histories of genes based on homology relations which are binary relations between genes. A variety of methods based on pairwise sequence comparisons can be used to infer such homology relations (e.g.\ orthology, paralogy, xenology). They are \emph{satisfiable} if the relations can be explained by an event-labeled gene tree, i.e., they can simultaneously co-exist in an evolutionary history of the underlying genes. Every gene tree is equivalently interpreted as a so-called cotree that entirely encodes the structure of a di-cograph. Thus, satisfiable homology relations must necessarily form a di-cograph. The inferred homology relations might not cover each pair of genes and thus, provide only partial knowledge on the full set of homology relations. Moreover, for particular pairs of genes, it might be known with a high degree of certainty that they are not orthologs (resp.\ paralogs, xenologs) which yields forbidden pairs of genes. Motivated by this observation, we characterize (partial) satisfiable homology relations with or without forbidden gene pairs, provide a quadratic-time algorithm for their recognition and for the computation of a cotree that explains the given relations

Canalization of the evolutionary trajectory of the human influenza virus

Since its emergence in 1968, influenza A (H3N2) has evolved extensively in genotype and antigenic phenotype. Antigenic evolution occurs in the context of a two-dimensional 'antigenic map', while genetic evolution shows a characteristic ladder-like genealogical tree. Here, we use a large-scale individual-based model to show that evolution in a Euclidean antigenic space provides a remarkable correspondence between model behavior and the epidemiological, antigenic, genealogical and geographic patterns observed in influenza virus. We find that evolution away from existing human immunity results in rapid population turnover in the influenza virus and that this population turnover occurs primarily along a single antigenic axis. Thus, selective dynamics induce a canalized evolutionary trajectory, in which the evolutionary fate of the influenza population is surprisingly repeatable and hence, in theory, predictable.Comment: 29 pages, 5 figures, 10 supporting figure

The genome and transcriptome of Trichormus sp NMC-1: insights into adaptation to extreme environments on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP) has the highest biodiversity for an extreme environment worldwide, and provides an ideal natural laboratory to study adaptive evolution. In this study, we generated a draft genome sequence of cyanobacteria Trichormus sp. NMC-1 in the QTP and performed whole transcriptome sequencing under low temperature to investigate the genetic mechanism by which T. sp. NMC-1 adapted to the specific environment. Its genome sequence was 5.9 Mb with a G+C content of 39.2% and encompassed a total of 5362 CDS. A phylogenomic tree indicated that this strain belongs to the Trichormus and Anabaena cluster. Genome comparison between T. sp. NMC-1 and six relatives showed that functionally unknown genes occupied a much higher proportion (28.12%) of the T. sp. NMC-1 genome. In addition, functions of specific, significant positively selected, expanded orthogroups, and differentially expressed genes involved in signal transduction, cell wall/membrane biogenesis, secondary metabolite biosynthesis, and energy production and conversion were analyzed to elucidate specific adaptation traits. Further analyses showed that the CheY-like genes, extracellular polysaccharide and mycosporine-like amino acids might play major roles in adaptation to harsh environments. Our findings indicate that sophisticated genetic mechanisms are involved in cyanobacterial adaptation to the extreme environment of the QTP

Genome landscapes and bacteriophage codon usage

Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonmous codons. Although alternative theories abound, translational selection has been accepted as an important mechanism that shapes the patterns of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns of codon usage across 74 diverse bacteriophages that infect E. coli, P. aeruginosa and L. lactis as their primary host. We introduce the concept of a `genome landscape,' which helps reveal non-trivial, long-range patterns in codon usage across a genome. We develop a series of randomization tests that allow us to interrogate the significance of one aspect of codon usage, such a GC content, while controlling for another aspect, such as adaptation to host-preferred codons. We find that 33 phage genomes exhibit highly non-random patterns in their GC3-content, use of host-preferred codons, or both. We show that the head and tail proteins of these phages exhibit significant bias towards host-preferred codons, relative to the non-structural phage proteins. Our results support the hypothesis of translational selection on viral genes for host-preferred codons, over a broad range of bacteriophages.Comment: 9 Color Figures, 5 Tables, 53 Reference

Genome-wide signatures of convergent evolution in echolocating mammals

Evolution is typically thought to proceed through divergence of genes, proteins, and ultimately phenotypes(1-3). However, similar traits might also evolve convergently in unrelated taxa due to similar selection pressures(4,5). Adaptive phenotypic convergence is widespread in nature, and recent results from a handful of genes have suggested that this phenomenon is powerful enough to also drive recurrent evolution at the sequence level(6-9). Where homoplasious substitutions do occur these have long been considered the result of neutral processes. However, recent studies have demonstrated that adaptive convergent sequence evolution can be detected in vertebrates using statistical methods that model parallel evolution(9,10) although the extent to which sequence convergence between genera occurs across genomes is unknown. Here we analyse genomic sequence data in mammals that have independently evolved echolocation and show for the first time that convergence is not a rare process restricted to a handful of loci but is instead widespread, continuously distributed and commonly driven by natural selection acting on a small number of sites per locus. Systematic analyses of convergent sequence evolution in 805,053 amino acids within 2,326 orthologous coding gene sequences compared across 22 mammals (including four new bat genomes) revealed signatures consistent with convergence in nearly 200 loci. Strong and significant support for convergence among bats and the dolphin was seen in numerous genes linked to hearing or deafness, consistent with an involvement in echolocation. Surprisingly we also found convergence in many genes linked to vision: the convergent signal of many sensory genes was robustly correlated with the strength of natural selection. This first attempt to detect genome-wide convergent sequence evolution across divergent taxa reveals the phenomenon to be much more pervasive than previously recognised

Automatically extracting functionally equivalent proteins from SwissProt

In summary, FOSTA provides an automated analysis of annotations in UniProtKB/Swiss-Prot to enable groups of proteins already annotated as functionally equivalent, to be extracted. Our results demonstrate that the vast majority of UniProtKB/Swiss-Prot functional annotations are of high quality, and that FOSTA can interpret annotations successfully. Where FOSTA is not successful, we are able to highlight inconsistencies in UniProtKB/Swiss-Prot annotation. Most of these would have presented equal difficulties for manual interpretation of annotations. We discuss limitations and possible future extensions to FOSTA, and recommend changes to the UniProtKB/Swiss-Prot format, which would facilitate text-mining of UniProtKB/Swiss-Prot

Sampling and counting genome rearrangement scenarios

Even for moderate size inputs, there are a tremendous number of optimal rearrangement scenarios, regardless what the model is and which specific question is to be answered. Therefore giving one optimal solution might be misleading and cannot be used for statistical inferring. Statistically well funded methods are necessary to sample uniformly from the solution space and then a small number of samples are sufficient for statistical inferring

A single evolutionary innovation drives the deep evolution of symbiotic N₂-fixation in angiosperms

Symbiotic associations occur in every habitat on earth, but we know very little about their evolutionary histories. Current models of trait evolution cannot adequately reconstruct the deep history of symbiotic innovation, because they assume homogenous evolutionary processes across millions of years. Here we use a recently developed, heterogeneous and quantitative phylogenetic framework to study the origin of the symbiosis between angiosperms and nitrogen-fixing (N2) bacterial symbionts housed in nodules. We compile the largest database of global nodulating plant species and reconstruct the symbiosis’ evolution. We identify a single, cryptic evolutionary innovation driving symbiotic N2-fixation evolution, followed by multiple gains and losses of the symbiosis, and the subsequent emergence of ‘stable fixers’ (clades extremely unlikely to lose the symbiosis). Originating over 100 MYA, this innovation suggests deep homology in symbiotic N2-fixation. Identifying cryptic innovations on the tree of life is key to understanding the evolution of complex traits, including symbiotic partnerships

Age related decline in female lar gibbon great call performance suggests that call features correlate with physical condition

Background: White-handed gibbons (Hylobates lar) are small Asian apes known for living in stable territories and producing loud, elaborate vocalizations (songs), often in well-coordinated male/female duets. The female great call, the most conspicuous phrase of the repertoire, has been hypothesized to function in intra-sexual territorial defense. We therefore predicted that characteristics of the great call would correlate with a caller’s physical condition, and thus might honestly reflect resource holding potential (RHP). Because measurement of RHP is virtually impossible for wild animals, we used age as a proxy, hypothesizing that great call climaxes are difficult to produce and maintain over time, and that older adults will therefore perform lower quality great calls than young adults. To test this we analyzed the great call climaxes of 15 wild lar gibbon females at Khao Yai National Park, Thailand and 2 captive females at Leo Conservation Center, Greenwich, CT. Results: Findings show that call climaxes correlate with female age, as young animals (n = 8, mean age: 12.9 years) produced climaxes with a higher frequency range (delta F0), maximum F0 frequency and duty cycle than old animals (n = 9, mean age: 29.6 years). A permuted discriminant function analysis also correctly classified calls by age group. During long song bouts the maximum F0 frequency of great call climaxes’ also decreased. Additional data support the hypothesis that short high notes, associated with rapid inhalation as an individual catches its breath, reflect increased caller effort. Older females produced more high notes than younger females, but the difference only approached statistical significance, suggesting that calling effort may be similar across different ages. Finally, for the first time in this species, we measured peak intensity of calls in captive females. They were capable of producing climaxes in excess of 100 dB at close range (2.7 m). Conclusions: Age and within-bout differences in the lar gibbon great call climax suggest that call features correlate with physical condition and thus the call may have evolved as an honest signal in the context of intra-sexual territorial defense and possibly also in male mate choice via sexual selection, although further testing of these hypotheses is necessary. Results: Findings show that call climaxes correlate with female age, as young animals (n = 8, mean age: 12.9 years) produced climaxes with a higher frequency range (delta F0), maximum F0 frequency and duty cycle than old animals (n = 9, mean age: 29.6 years). A permuted discriminant function analysis also correctly classified calls by age group. During long song bouts the maximum F0 frequency of great call climaxes’ also decreased. Additional data support the hypothesis that short high notes, associated with rapid inhalation as an individual catches its breath, reflect increased caller effort. Older females produced more high notes than younger females, but the difference only approached statistical significance, suggesting that calling effort may be similar across different ages. Finally, for the first time in this species, we measured peak intensity of calls in captive females. They were capable of producing climaxes in excess of 100 dB at close range (2.7 m). Conclusions: Age and within-bout differences in the lar gibbon great call climax suggest that call features correlate with physical condition and thus the call may have evolved as an honest signal in the context of intra-sexual territorial defense and possibly also in male mate choice via sexual selection, although further testing of these hypotheses is necessary