Autobiography, Biography, and Narrative Ethics

Background: Duplicated genes can indefinately persist in genomes if either both copies retain the original function due to dosage benefit (gene conservation), or one of the copies assumes a novel function (neofunctionalization), or both copies become required to perform the function previously accomplished by a single copy (subfunctionalization), or through a combination of these mechanisms. Different models of duplication retention imply different predictions about substitution rates in the coding portion of paralogs and about asymmetry of these rates. Results: We analyse sequence evolution asymmetry in paralogs present in 12 Drosophila genomes using the nearest non-duplicated orthologous outgroup as a reference. Those paralogs present in D. melanogaster are analysed in conjunction with the asymmetry of expression rate and ubiquity and of segregating non-synonymous polymorphisms in the same paralogs. Paralogs accumulate substitutions, on average, faster than their nearest singleton orthologs. The distribution of paralogs\u27 substitution rate asymmetry is overdispersed relative to that of orthologous clades, containing disproportionally more unusually symmetric and unusually asymmetric clades. We show that paralogs are more asymmetric in: a) clades orthologous to highly constrained singleton genes; b) genes with high expression level; c) genes with ubiquitous expression and d) non-tandem duplications. We further demonstrate that, in each asymmetrically evolving pair of paralogs, the faster evolving member of the pair tends to have lower average expression rate, lower expression uniformity and higher frequency of non-synonymous SNPs than its slower evolving counterpart. Conclusions: Our findings are consistent with the hypothesis that many duplications in Drosophila are retained despite stabilising selection being more relaxed in one of the paralogs than in the other, suggesting a widespread unfinished pseudogenization. This phenomenon is likely to make detection of neo- and subfunctionalization signatures difficult, as these models of duplication retention also predict asymmetries in substitution rates and expression profiles. Reviewers: This article has been reviewed by Dr. Jia Zeng (nominated by Dr. I. King Jordan), Dr. Fyodor Kondrashov and Dr. Yuri Wolf

Evolutionary Patterns of Amino Acid Substitutions in 12 Drosophila Genomes

Background: Harnessing 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\u27 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.Results: We 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.Conclusions: We 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

Evolutionary patterns of amino acid substitutions in 12 Drosophila genomes

<p>Abstract</p> <p>Background</p> <p>Harnessing 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 <it>Drosophila</it> 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.</p> <p>Results</p> <p>We demonstrate that amino acid composition of <it>Drosophila</it> 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 <it>Drosophila</it> 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.</p> <p>Conclusions</p> <p>We 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.</p

Widespread Transcriptional Autosomal Dosage Compensation in Drosophila Correlates With Gene Expression Level

Little is known about dosage compensation in autosomal genes. Transcription-level compensation of deletions and other loss-of-function mutations may be a mechanism of dominance of wild-type alleles, a ubiquitous phenomenon whose nature has been a subject of a long debate. We measured gene expression in two isogenic Drosophila lines heterozygous for long deletions and compared our results with previously published gene expression data in a line heterozygous for a long duplication. We find that a majority of genes are at least partially compensated at transcription, both for 1/2-fold dosage (in heterozygotes for deletions) and for 1.5-fold dosage (in heterozygotes for a duplication). The degree of compensation does not vary among functional classes of genes. Compensation for deletions is stronger for highly expressed genes. In contrast, the degree of compensation for duplications is stronger for weakly expressed genes. Thus, partial transcriptional compensation appears to be based on regulatory mechanisms that insure high transcription levels of some genes and low transcription levels of other genes, instead of precise maintenance of a particular homeostatic expression level. Given the ubiquity of transcriptional compensation, dominance of wild-type alleles may be at least partially caused by of the regulation at transcription level

Autobiography, Biography, and Narrative Ethics

Background: Duplicated genes can indefinately persist in genomes if either both copies retain the original function due to dosage benefit (gene conservation), or one of the copies assumes a novel function (neofunctionalization), or both copies become required to perform the function previously accomplished by a single copy (subfunctionalization), or through a combination of these mechanisms. Different models of duplication retention imply different predictions about substitution rates in the coding portion of paralogs and about asymmetry of these rates. Results: We analyse sequence evolution asymmetry in paralogs present in 12 Drosophila genomes using the nearest non-duplicated orthologous outgroup as a reference. Those paralogs present in D. melanogaster are analysed in conjunction with the asymmetry of expression rate and ubiquity and of segregating non-synonymous polymorphisms in the same paralogs. Paralogs accumulate substitutions, on average, faster than their nearest singleton orthologs. The distribution of paralogs\u27 substitution rate asymmetry is overdispersed relative to that of orthologous clades, containing disproportionally more unusually symmetric and unusually asymmetric clades. We show that paralogs are more asymmetric in: a) clades orthologous to highly constrained singleton genes; b) genes with high expression level; c) genes with ubiquitous expression and d) non-tandem duplications. We further demonstrate that, in each asymmetrically evolving pair of paralogs, the faster evolving member of the pair tends to have lower average expression rate, lower expression uniformity and higher frequency of non-synonymous SNPs than its slower evolving counterpart. Conclusions: Our findings are consistent with the hypothesis that many duplications in Drosophila are retained despite stabilising selection being more relaxed in one of the paralogs than in the other, suggesting a widespread unfinished pseudogenization. This phenomenon is likely to make detection of neo- and subfunctionalization signatures difficult, as these models of duplication retention also predict asymmetries in substitution rates and expression profiles. Reviewers: This article has been reviewed by Dr. Jia Zeng (nominated by Dr. I. King Jordan), Dr. Fyodor Kondrashov and Dr. Yuri Wolf

Pou5f3, SoxB1, and Nanog Remodel Chromatin on High Nucleosome Affinity Regions at Zygotic Genome Activation

The zebrafish embryo is transcriptionally mostly quiescent during the first 10 cell cycles, until the main wave of zygotic genome activation (ZGA) occurs, accompanied by fast chromatin remodeling. At ZGA, homologs of the mammalian stem cell transcription factors (TFs) Pou5f3, Nanog, and Sox19b bind to thousands of developmental enhancers to initiate transcription. So far, how these TFs influence chromatin dynamics at ZGA has remained unresolved. To address this question, we analyzed nucleosome positions in wild-type and maternal-zygotic (MZ) mutants for pou5f3 and nanog by MNase-seq. We show that Nanog, Sox19b, and Pou5f3 bind to the high nucleosome affinity regions (HNARs). HNARs are spanning over 600 bp, featuring high in vivo and predicted in vitro nucleosome occupancy and high predicted propeller twist DNA shape value. We suggest a two-step nucleosome destabilization-depletion model, in which the same intrinsic DNA properties of HNAR promote both high nucleosome occupancy and differential binding of TFs. In the first step, already before ZGA, Pou5f3 and Nanog destabilize nucleosomes at HNAR centers genome-wide. In the second step, post-ZGA, Nanog, Pou5f3, and SoxB1 maintain open chromatin state on the subset of HNARs, acting synergistically. Nanog binds to the HNAR center, whereas the Pou5f3 stabilizes the flanks. The HNAR model will provide a useful tool for genome regulatory studies in a variety of biological systems

Warming Overcomes Dispersal-Limitation to Promote Non-native Expansion in Lake Baikal

Non-native species and climate change pose serious threats to global biodiversity. However, the roles of climate, dispersal, and competition are difficult to disentangle in heterogeneous landscapes. We combine empirical data and theory to examine how these forces influence the spread of non-native species in Lake Baikal. We analyze the potential for Daphnia longispina to establish in Lake Baikal, potentially threatening an endemic, cryophillic copepod Epischurella baikalensis. We collected field samples to establish current community composition and compared them to model predictions informed by flow rates, present-day temperatures, and temperature projections. Our data and model agree that expansion is currently limited by dispersal. However, projected increases in temperature reverse this effect, allowing D. longispina to establish in Lake Baikal’s main basin. A strong negative impact emerges from the interaction between climate change and dispersal, outweighing their independent effects. Climate, dispersal, and competition have complex, interactive effects on expansion with important implications for global biodiversity

Functional genomics of acclimation and adaptation in response to thermal stress in Daphnia

BACKGROUND: Gene expression regulation is one of the fundamental mechanisms of phenotypic plasticity and is expected to respond to selection in conditions favoring phenotypic response. The observation that many organisms increase their stress tolerance after acclimation to moderate levels of stress is an example of plasticity which has been long hypothesized to be based on adaptive changes in gene expression. We report genome-wide patterns of gene expression in two heat-tolerant and two heat-sensitive parthenogenetic clones of the zooplankton crustacean Daphnia pulex exposed for three generations to either optimal (18°C) or substressful (28°C) temperature. RESULTS: A large number of genes responded to temperature and many demonstrated a significant genotype-by-environment (GxE) interaction. Among genes with a significant GxE there were approximately equally frequent instances of canalization, i.e. stronger plasticity in heat-sensitive than in heat-tolerant clones, and of enhancement of plasticity along the evolutionary vector toward heat tolerance. The strongest response observed is the across-the-board down-regulation of a variety of genes occurring in heat-tolerant, but not in heat-sensitive clones. This response is particularly obvious among genes involved in core metabolic pathways and those responsible for transcription, translation and DNA repair. CONCLUSIONS: The observed down-regulation of metabolism, consistent with previous findings in yeast and Drosophila, may reflect a general compensatory stress response. The associated down-regulation of DNA repair pathways potentially creates a trade-off between short-term benefits of survival at high temperature and long-term costs of accelerated mutation accumulation. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/1471-2164-15-859) contains supplementary material, which is available to authorized users

Functional Genomics of Acclimation and Adaptation in Response to Thermal Stress in Daphnia

Background: Gene expression regulation is one of the fundamental mechanisms of phenotypic plasticity and is expected to respond to selection in conditions favoring phenotypic response. The observation that many organisms increase their stress tolerance after acclimation to moderate levels of stress is an example of plasticity which has been long hypothesized to be based on adaptive changes in gene expression. We report genome-wide patterns of gene expression in two heat-tolerant and two heat-sensitive parthenogenetic clones of the zooplankton crustacean Daphnia pulex exposed for three generations to either optimal (18°C) or substressful (28°C) temperature. Results: A large number of genes responded to temperature and many demonstrated a significant genotype-byenvironment (GxE) interaction. Among genes with a significant GxE there were approximately equally frequent instances of canalization, i.e. stronger plasticity in heat-sensitive than in heat-tolerant clones, and of enhancement of plasticity along the evolutionary vector toward heat tolerance. The strongest response observed is the across-the-board down-regulation of a variety of genes occurring in heat-tolerant, but not in heat-sensitive clones. This response is particularly obvious among genes involved in core metabolic pathways and those responsible for transcription, translation and DNA repair. Conclusions: The observed down-regulation of metabolism, consistent with previous findings in yeast and Drosophila, may reflect a general compensatory stress response. The associated down-regulation of DNA repair pathways potentially creates a trade-off between short-term benefits of survival at high temperature and long-term costs of accelerated mutation accumulation

Transcriptome‐based phylogeny of endemic Lake Baikal amphipod species flock: fast speciation accompanied by frequent episodes of positive selection

Endemic species flocks inhabiting ancient lakes, oceanic islands and other long‐lived isolated habitats are often interpreted as adaptive radiations. Yet molecular evidence for directional selection during species flocks radiation is scarce. Using partial transcriptomes of 64 species of Lake Baikal (Siberia, Russia) endemic amphipods and two nonendemic outgroups, we report a revised phylogeny of this species flock and analyse evidence for positive selection within the endemic lineages. We confirm two independent invasions of amphipods into Baikal and demonstrate that several morphological features of Baikal amphipods, such as body armour and reduction in appendages and sensory organs, evolved in several lineages in parallel. Radiation of Baikal amphipods has been characterized by short phylogenetic branches and frequent episodes of positive selection which tended to be more frequent in the early phase of the second invasion of amphipods into Baikal when the most intensive diversification occurred. Notably, signatures of positive selection are frequent in genes encoding mitochondrial membrane proteins with electron transfer chain and ATP synthesis functionality. In particular, subunits of both the membrane and substrate‐level ATP synthases show evidence of positive selection in the plankton species Macrohectopus branickii, possibly indicating adaptation to active plankton lifestyle and to survival under conditions of low temperature and high hydrostatic pressures known to affect membranes functioning. Other functional categories represented among genes likely to be under positive selection include Ca‐binding muscle‐related proteins, possibly indicating adaptation to Ca‐deficient low mineralization Baikal waters.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/136009/1/mec13927.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/136009/2/mec13927_am.pd