Evolution and Survival on Eutherian Sex Chromosomes

Since the two eutherian sex chromosomes diverged from an ancestral autosomal pair, the X has remained relatively gene-rich, while the Y has lost most of its genes through the accumulation of deleterious mutations in nonrecombining regions. Presently, it is unclear what is distinctive about genes that remain on the Y chromosome, when the sex chromosomes acquired their unique evolutionary rates, and whether X-Y gene divergence paralleled that of paralogs located on autosomes. To tackle these questions, here we juxtaposed the evolution of X and Y homologous genes (gametologs) in eutherian mammals with their autosomal orthologs in marsupial and monotreme mammals. We discovered that genes on the X and Y acquired distinct evolutionary rates immediately following the suppression of recombination between the two sex chromosomes. The Y-linked genes evolved at higher rates, while the X-linked genes maintained the lower evolutionary rates of the ancestral autosomal genes. These distinct rates have been maintained throughout the evolution of X and Y. Specifically, in humans, most X gametologs and, curiously, also most Y gametologs evolved under stronger purifying selection than similarly aged autosomal paralogs. Finally, after evaluating the current experimental data from the literature, we concluded that unique mRNA/protein expression patterns and functions acquired by Y (versus X) gametologs likely contributed to their retention. Our results also suggest that either the boundary between sex chromosome strata 3 and 4 should be shifted or that stratum 3 should be divided into two strata

Evolution of Parallel Spindles Like genes in plants and highlight of unique domain architecture#

<p>Abstract</p> <p>Background</p> <p>Polyploidy has long been recognized as playing an important role in plant evolution. In flowering plants, the major route of polyploidization is suggested to be sexual through gametes with somatic chromosome number (2<it>n</it>). <it>Parallel Spindle1 </it>gene in <it>Arabidopsis thaliana </it>(<it>AtPS1</it>) was recently demonstrated to control spindle orientation in the 2nd division of meiosis and, when mutated, to induce 2<it>n </it>pollen. Interestingly, <it>AtPS1 </it>encodes a protein with a FHA domain and PINc domain putatively involved in RNA decay (i.e. Nonsense Mediated mRNA Decay). In potato, 2<it>n </it>pollen depending on parallel spindles was described long time ago but the responsible gene has never been isolated. The knowledge derived from <it>AtPS1 </it>as well as the availability of genome sequences makes it possible to isolate potato <it>PSLike </it>(<it>PSL</it>) and to highlight the evolution of <it>PSL </it>family in plants.</p> <p>Results</p> <p>Our work leading to the first characterization of <it>PSLs </it>in potato showed a greater <it>PSL </it>complexity in this species respect to <it>Arabidopsis thaliana</it>. Indeed, a genomic <it>PSL </it>locus and seven cDNAs affected by alternative splicing have been cloned. In addition, the occurrence of at least two other <it>PSL </it>loci in potato was suggested by the sequence comparison of alternatively spliced transcripts.</p> <p>Phylogenetic analysis on 20 <it>Viridaeplantae </it>showed the wide distribution of <it>PSLs </it>throughout the species and the occurrence of multiple copies only in potato and soybean.</p> <p>The analysis of PSL^FHAand PSL^PINcdomains evidenced that, in terms of secondary structure, a major degree of variability occurred in PINc domain respect to FHA. In terms of specific active sites, both domains showed diversification among plant species that could be related to a functional diversification among <it>PSL </it>genes. In addition, some specific active sites were strongly conserved among plants as supported by sequence alignment and by evidence of negative selection evaluated as difference between non-synonymous and synonymous mutations.</p> <p>Conclusions</p> <p>In this study, we highlight the existence of PSLs throughout <it>Viridaeplantae</it>, from mosses to higher plants. We provide evidence that <it>PSLs </it>occur mostly as singleton in the analyzed genomes except in soybean and potato both characterized by a recent whole genome duplication event. In potato, we suggest the candidate <it>PSL </it>gene having a role in 2<it>n </it>pollen that should be deeply investigated.</p> <p>We provide useful insight into evolutionary conservation of FHA and PINc domains throughout plant PSLs which suggest a fundamental role of these domains for PSL function.</p

HIV-Specific Probabilistic Models of Protein Evolution

Comparative sequence analyses, including such fundamental bioinformatics techniques as similarity searching, sequence alignment and phylogenetic inference, have become a mainstay for researchers studying type 1 Human Immunodeficiency Virus (HIV-1) genome structure and evolution. Implicit in comparative analyses is an underlying model of evolution, and the chosen model can significantly affect the results. In general, evolutionary models describe the probabilities of replacing one amino acid character with another over a period of time. Most widely used evolutionary models for protein sequences have been derived from curated alignments of hundreds of proteins, usually based on mammalian genomes. It is unclear to what extent these empirical models are generalizable to a very different organism, such as HIV-1–the most extensively sequenced organism in existence. We developed a maximum likelihood model fitting procedure to a collection of HIV-1 alignments sampled from different viral genes, and inferred two empirical substitution models, suitable for describing between-and within-host evolution. Our procedure pools the information from multiple sequence alignments, and provided software implementation can be run efficiently in parallel on a computer cluster. We describe how the inferred substitution models can be used to generate scoring matrices suitable for alignment and similarity searches. Our models had a consistently superior fit relative to the best existing models and to parameter-rich data-driven models when benchmarked on independent HIV-1 alignments, demonstrating evolutionary biases in amino-acid substitution that are unique to HIV, and that are not captured by the existing models. The scoring matrices derived from the models showed a marked difference from common amino-acid scoring matrices. The use of an appropriate evolutionary model recovered a known viral transmission history, whereas a poorly chosen model introduced phylogenetic error. We argue that our model derivation procedure is immediately applicable to other organisms with extensive sequence data available, such as Hepatitis C and Influenza A viruses

Finding Direction in the Search for Selection.

Tests for positive selection have mostly been developed to look for diversifying selection where change away from the current amino acid is often favorable. However, in many cases we are interested in directional selection where there is a shift toward specific amino acids, resulting in increased fitness in the species. Recently, a few methods have been developed to detect and characterize directional selection on a molecular level. Using the results of evolutionary simulations as well as HIV drug resistance data as models of directional selection, we compare two such methods with each other, as well as against a standard method for detecting diversifying selection. We find that the method to detect diversifying selection also detects directional selection under certain conditions. One method developed for detecting directional selection is powerful and accurate for a wide range of conditions, while the other can generate an excessive number of false positives

Ancient Ancestry of KFDV and AHFV Revealed by Complete Genome Analyses of Viruses Isolated from Ticks and Mammalian Hosts

Alkhurma hemorrhagic fever (AHF) and Kyasanur Forest disease (KFD) viruses both cause serious and sometimes fatal human disease in their respective ranges, Saudi Arabia and India. AHFV was first identified in the mid-1990s and due to its strong genetic similarity to KFDV it has since been considered the result of a recent introduction of KFDV into Saudi Arabia. To gain a better understanding of the evolutionary history of AHFV and KFDV, we sequenced the full-length genomes of 3 KFDV and 16 AHFV. Sequence analyses show a greater genetic diversity within AHFV than previously thought, particularly within the tick population. The phylogeny constructed with these 19 full-length sequences and two AHFV sequences from GenBank indicates AHFV diverged from KFDV almost 700 years ago. Given the presence of competent tick vectors in the regions between and surrounding Saudi Arabia and India and the recent identification of AHFV in Egypt, these results suggest a broader geographic range of AHFV and KFDV, and raise the possibility of other AHFV/KFDV–like viruses circulating in these regions

Assessing constancy of substitution rates in viruses over evolutionary time

<p>Abstract</p> <p>Background</p> <p>Phylogenetic analyses reveal probable patterns of divergence of present day organisms from common ancestors. The points of divergence of lineages can be dated if a corresponding historical or fossil record exists. For many species, in particular viruses, such records are rare. Recently, Bayesian phylogenetic analysis using sequences from closely related organisms isolated at different times have been used to calibrate divergences. Phylogenetic analyses depend on the assumption that the average substitution rates that can be calculated from the data apply throughout the course of evolution. </p> <p>Results</p> <p>The present study tests this crucial assumption by charting the kinds of substitutions observed between pairs of sequences with different levels of total substitutions. Datasets of aligned sequences, both viral and non-viral, were assembled. For each pair of sequences in an aligned set, the distribution of nucleotide interchanges and the total number of changes were calculated. Data were binned according to total numbers of changes and plotted. The accumulation of the six possible interchange types in retroelements as a function of distance followed closely the expected hyperbolic relationship. For other datasets, however, significant deviations from this relationship were noted. A rapid initial accumulation of transition interchanges was frequent among the datasets and anomalous changes occurred at specific divergence levels. </p> <p>Conclusions</p> <p>The accumulation profiles suggested that substantial changes in frequencies of types of substitutions occur over the course of evolution and that such changes should be considered in evaluating and dating viral phylogenies.</p

Genetic Signature of Rapid IHHNV (Infectious Hypodermal and Hematopoietic Necrosis Virus) Expansion in Wild Penaeus Shrimp Populations

Infectious hypodermal and hematopoietic necrosis virus (IHHNV) is a widely distributed single-stranded DNA parvovirus that has been responsible for major losses in wild and farmed penaeid shrimp populations on the northwestern Pacific coast of Mexico since the early 1990's. IHHNV has been considered a slow-evolving, stable virus because shrimp populations in this region have recovered to pre-epizootic levels, and limited nucleotide variation has been found in a small number of IHHNV isolates studied from this region. To gain insight into IHHNV evolutionary and population dynamics, we analyzed IHHNV capsid protein gene sequences from 89 Penaeus shrimp, along with 14 previously published sequences. Using Bayesian coalescent approaches, we calculated a mean rate of nucleotide substitution for IHHNV that was unexpectedly high (1.39×10−4 substitutions/site/year) and comparable to that reported for RNA viruses. We found more genetic diversity than previously reported for IHHNV isolates and highly significant subdivision among the viral populations in Mexican waters. Past changes in effective number of infections that we infer from Bayesian skyline plots closely correspond to IHHNV epizootiological historical records. Given the high evolutionary rate and the observed regional isolation of IHHNV in shrimp populations in the Gulf of California, we suggest regular monitoring of wild and farmed shrimp and restriction of shrimp movement as preventative measures for future viral outbreaks

A new look at the LTR retrotransposon content of the chicken genome

BACKGROUND: LTR retrotransposons contribute approximately 10 % of the mammalian genome, but it has been previously reported that there is a deficit of these elements in the chicken relative to both mammals and other birds. A novel LTR retrotransposon classification pipeline, LocaTR, was developed and subsequently utilised to re-examine the chicken LTR retrotransposon annotation, and determine if the proposed chicken deficit is biologically accurate or simply a technical artefact. RESULTS: Using LocaTR 3.01 % of the chicken galGal4 genome assembly was annotated as LTR retrotransposon-derived elements (nearly double the previous annotation), including 1,073 that were structurally intact. Element distribution is significantly correlated with chromosome size and is non-random within each chromosome. Elements are significantly depleted within coding regions and enriched in gene sparse areas of the genome. Over 40 % of intact elements are found in clusters, unrelated by age or genera, generally in poorly recombining regions. The transcription of most LTR retrotransposons were suppressed or incomplete, but individual domain and full length retroviral transcripts were produced in some cases, although mostly with regularly interspersed stop codons in all reading frames. Furthermore, RNAseq data from 23 diverse tissues enabled greater characterisation of the co-opted endogenous retrovirus Ovex1. This gene was shown to be expressed ubiquitously but at variable levels across different tissues. LTR retrotransposon content was found to be very variable across the avian lineage and did not correlate with either genome size or phylogenetic position. However, the extent of previous, species-specific LTR retrotransposon annotation appears to be a confounding factor. CONCLUSIONS: Use of the novel LocaTR pipeline has nearly doubled the annotated LTR retrotransposon content of the chicken genome compared to previous estimates. Further analysis has described element distribution, clustering patterns and degree of expression in a variety of adult tissues, as well as in three embryonic stages. This study also enabled better characterisation of the co-opted gamma retroviral envelope gene Ovex1. Additionally, this work suggests that there is no deficit of LTR retrotransposons within the Galliformes relative to other birds, or to mammalian genomes when scaled for the three-fold difference in genome size

A Model for Damage Load and Its Implications for the Evolution of Bacterial Aging

Deleterious mutations appearing in a population increase in frequency until stopped by natural selection. The ensuing equilibrium creates a stable frequency of deleterious mutations or the mutational load. Here I develop the comparable concept of a damage load, which is caused by harmful non-heritable changes to the phenotype. A damage load also ensues when the increase of damage is opposed by selection. The presence of a damage load favors the evolution of asymmetrical transmission of damage by a mother to her daughters. The asymmetry is beneficial because it increases fitness variance, but it also leads to aging or senescence. A mathematical model based on microbes reveals that a cell lineage dividing symmetrically is immortal if lifetime damage rates do not exceed a threshold. The evolution of asymmetry allows the lineage to persist above the threshold, but the lineage becomes mortal. In microbes with low genomic mutation rates, it is likely that the damage load is much greater than the mutational load. In metazoans with higher genomic mutation rates, the damage and the mutational load could be of the same magnitude. A fit of the model to experimental data shows that Escherichia coli cells experience a damage rate that is below the threshold and are immortal under the conditions examined. The model estimates the asymmetry level of E. coli to be low but sufficient for persisting at higher damage rates. The model also predicts that increasing asymmetry results in diminishing fitness returns, which may explain why the bacterium has not evolved higher asymmetry