Explore the vast range of titles available.

Proteins in the slow lane Protein evolution occurs only slowly, since most amino-acid substitutions are likely to be deleterious and selection will favour conservation of function. Inna Povolotskaya and Fyodor Kondrashov set out to find out just how slowly by asking the question, are ancient extant proteins — those that were present in the last universal common ancestor (LUCA) — continuing to diverge from the ancestral sequence? Their calculations, based on the approach used by Edwin Hubble in his study of the recession of galaxies in the physical Universe, suggests that extant protein sequences are still expanding from each other, and, therefore, from their common ancestor. Divergence is very slow: the 3.5 billion years or so since LUCA's time has not been long enough for the limit of sequence divergence to be reached. This tardiness is a consequence of the sparseness of functional protein sequences in sequence space and the ruggedness of the protein fitness landscape: 98% of sites cannot accept an amino-acid substitution at a given moment, yet most sites may eventually be permitted to evolve when other, compensatory, changes occur. The need to maintain the structural and functional integrity of an evolving protein limits the range of acceptable amino-acid substitutions — but to what extent does this constrain how far homologous protein sequences can diverge? Here, sequence divergence data are used to explore the limits of protein evolution, and to conclude that ancient proteins are continuing to diverge from one another, indicating that the protein sequence universe is slowly expanding. The need to maintain the structural and functional integrity of an evolving protein severely restricts the repertoire of acceptable amino-acid substitutions 1 , 2 , 3 , 4 . However, it is not known whether these restrictions impose a global limit on how far homologous protein sequences can diverge from each other. Here we explore the limits of protein evolution using sequence divergence data. We formulate a computational approach to study the rate of divergence of distant protein sequences and measure this rate for ancient proteins, those that were present in the last universal common ancestor. We show that ancient proteins are still diverging from each other, indicating an ongoing expansion of the protein sequence universe. The slow rate of this divergence is imposed by the sparseness of functional protein sequences in sequence space and the ruggedness of the protein fitness landscape: ∼98 per cent of sites cannot accept an amino-acid substitution at any given moment but a vast majority of all sites may eventually be permitted to evolve when other, compensatory, changes occur. Thus, ∼3.5 × 10 9  yr has not been enough to reach the limit of divergent evolution of proteins, and for most proteins the limit of sequence similarity imposed by common function may not exceed that of random sequences.

Characterizing the fitness landscape, a representation of fitness for a large set of genotypes, is key to understanding how genetic information is interpreted to create functional organisms. Here we determined the evolutionarily-relevant segment of the fitness landscape of His3, a gene coding for an enzyme in the histidine synthesis pathway, focusing on combinations of amino acid states found at orthologous sites of extant species. Just 15% of amino acids found in yeast His3 orthologues were always neutral while the impact on fitness of the remaining 85% depended on the genetic background. Furthermore, at 67% of sites, amino acid replacements were under sign epistasis, having both strongly positive and negative effect in different genetic backgrounds. 46% of sites were under reciprocal sign epistasis. The fitness impact of amino acid replacements was influenced by only a few genetic backgrounds but involved interaction of multiple sites, shaping a rugged fitness landscape in which many of the shortest paths between highly fit genotypes are inaccessible.

The rhesus macaque is an abundant species of Old World monkeys and a valuable model organism for biomedical research due to its close phylogenetic relationship to humans. Copy number variation is one of the main sources of genomic diversity within and between species and a widely recognized cause of inter-individual differences in disease risk. However, copy number differences among rhesus macaques and between the human and macaque genomes, as well as the relevance of this diversity to research involving this nonhuman primate, remain understudied. Here we present a high-resolution map of sequence copy number for the rhesus macaque genome constructed from a dataset of 198 individuals. Our results show that about one-eighth of the rhesus macaque reference genome is composed of recently duplicated regions, either copy number variable regions or fixed duplications. Comparison with human genomic copy number maps based on previously published data shows that, despite overall similarities in the genome-wide distribution of these regions, there are specific differences at the chromosome level. Some of these create differences in the copy number profile between human disease genes and their rhesus macaque orthologs. Our results highlight the importance of addressing the number of copies of target genes in the design of experiments and cautions against human-centered assumptions in research conducted with model organisms. Overall, we present a genome-wide copy number map from a large sample of rhesus macaque individuals representing an important novel contribution concerning the evolution of copy number in primate genomes.

Background Whole genome re-sequencing data from dogs and wolves are now commonly used to study how natural and artificial selection have shaped the patterns of genetic diversity. Single nucleotide polymorphisms, microsatellites and variants in mitochondrial DNA have been interrogated for links to specific phenotypes or signals of domestication. However, copy number variation (CNV), despite its increasingly recognized importance as a contributor to phenotypic diversity, has not been extensively explored in canids. Results Here, we develop a new accurate probabilistic framework to create fine-scale genomic maps of segmental duplications (SDs), compare patterns of CNV across groups and investigate their role in the evolution of the domestic dog by using information from 34 canine genomes. Our analyses show that duplicated regions are enriched in genes and hence likely possess functional importance. We identify 86 loci with large CNV differences between dogs and wolves, enriched in genes responsible for sensory perception, immune response, metabolic processes, etc. In striking contrast to the observed loss of nucleotide diversity in domestic dogs following the population bottlenecks that occurred during domestication and breed creation, we find a similar proportion of CNV loci in dogs and wolves, suggesting that other dynamics are acting to particularly select for CNVs with potentially functional impacts. Conclusions This work is the first comparison of genome wide CNV patterns in domestic and wild canids using whole-genome sequencing data and our findings contribute to study the impact of novel kinds of genetic changes on the evolution of the domestic dog.

The origins of neural systems remain unresolved. In contrast to other basal metazoans, ctenophores (comb jellies) have both complex nervous and mesoderm-derived muscular systems. These holoplanktonic predators also have sophisticated ciliated locomotion, behaviour and distinct development. Here we present the draft genome of Pleurobrachia bachei , Pacific sea gooseberry, together with ten other ctenophore transcriptomes, and show that they are remarkably distinct from other animal genomes in their content of neurogenic, immune and developmental genes. Our integrative analyses place Ctenophora as the earliest lineage within Metazoa. This hypothesis is supported by comparative analysis of multiple gene families, including the apparent absence of HOX genes, canonical microRNA machinery, and reduced immune complement in ctenophores. Although two distinct nervous systems are well recognized in ctenophores, many bilaterian neuron-specific genes and genes of ‘classical’ neurotransmitter pathways either are absent or, if present, are not expressed in neurons. Our metabolomic and physiological data are consistent with the hypothesis that ctenophore neural systems, and possibly muscle specification, evolved independently from those in other animals. The draft genome of the ctenophore Pleurobrachia bachei (Pacific sea gooseberry) is presented, together with ten other ctenophore transcriptomes — these genomes have a very different neurogenic, immune and developmental gene content when compared with other animal genomes, and it is proposed that ctenophore neural systems, and possibly muscle specification, evolved independently from those in other animals. Life, not quite as we know it Ctenophores (comb jellies) are enigmatic animals that combine two distinct nervous systems with an elementary brain-like centre and possess mesoderm-derived muscles appropriate to their predatory lifestyle. Leonid Moroz et al . present the draft genome of the ctenophore Pleurobrachia bachei (Pacific sea gooseberry), together with ten other ctenophore transcriptomes. These genomes have a neurogenic, immune and developmental gene content that differs markedly from other animal genomes: HOX genes and canonical microRNA machinery are absent, and the immune gene complement is reduced. Many bilaterian neuron-specific genes and genes of 'classical' neurotransmitter pathways either are absent or not expressed in neurons. The authors propose that ctenophore neural systems, and possibly muscle specification, evolved independently from those in other animals.

Bioluminescence is found across the entire tree of life, conferring a spectacular set of visually oriented functions from attracting mates to scaring off predators. Half a dozen different luciferins, molecules that emit light when enzymatically oxidized, are known. However, just one biochemical pathway for luciferin biosynthesis has been described in full, which is found only in bacteria. Here, we report identification of the fungal luciferase and three other key enzymes that together form the biosynthetic cycle of the fungal luciferin from caffeic acid, a simple and widespread metabolite. Introduction of the identified genes into the genome of the yeast Pichia pastoris along with caffeic acid biosynthesis genes resulted in a strain that is autoluminescent in standard media. We analyzed evolution of the enzymes of the luciferin biosynthesis cycle and found that fungal bioluminescence emerged through a series of events that included two independent gene duplications. The retention of the duplicated enzymes of the luciferin pathway in nonluminescent fungi shows that the gene duplication was followed by functional sequence divergence of enzymes of at least one gene in the biosynthetic pathway and suggests that the evolution of fungal bioluminescence proceeded through several closely related stepping stone nonluminescent biochemical reactions with adaptive roles. The availability of a complete eukaryotic luciferin biosynthesis pathway provides several applications in biomedicine and bioengineering.

Background The evolution and genomic stop codon frequencies have not been rigorously studied with the exception of coding of non-canonical amino acids. Here we study the rate of evolution and frequency distribution of stop codons in bacterial genomes. Results We show that in bacteria stop codons evolve slower than synonymous sites, suggesting the action of weak negative selection. However, the frequency of stop codons relative to genomic nucleotide content indicated that this selection regime is not straightforward. The frequency of TAA and TGA stop codons is GC-content dependent, with TAA decreasing and TGA increasing with GC-content, while TAG frequency is independent of GC-content. Applying a formal, analytical model to these data we found that the relationship between stop codon frequencies and nucleotide content cannot be explained by mutational biases or selection on nucleotide content. However, with weak nucleotide content-dependent selection on TAG, -0.5 < Nes < 1.5, the model fits all of the data and recapitulates the relationship between TAG and nucleotide content. For biologically plausible rates of mutations we show that, in bacteria, TAG stop codon is universally associated with lower fitness, with TAA being the optimal for G-content < 16% while for G-content > 16% TGA has a higher fitness than TAG. Conclusions Our data indicate that TAG codon is universally suboptimal in the bacterial lineage, such that TAA is likely to be the preferred stop codon for low GC content while the TGA is the preferred stop codon for high GC content. The optimization of stop codon usage may therefore be useful in genome engineering or gene expression optimization applications. Reviewers This article was reviewed by Michail Gelfand, Arcady Mushegian and Shamil Sunyaev. For the full reviews, please go to the Reviewers’ Comments section.

The genomic changes underlying both early and late stages of horse domestication remain largely unknown. We examined the genomes of 14 early domestic horses from the Bronze and Iron Ages, dating to between ~4.1 and 2.3 thousand years before present. We find early domestication selection patterns supporting the neural crest hypothesis, which provides a unified developmental origin for common domestic traits. Within the past 2.3 thousand years, horses lost genetic diversity and archaic DNA tracts introgressed from a now-extinct lineage. They accumulated deleterious mutations later than expected under the cost-of-domestication hypothesis, probably because of breeding from limited numbers of stallions. We also reveal that Iron Age Scythian steppe nomads implemented breeding strategies involving no detectable inbreeding and selection for coat-color variation and robust forelimbs.

Background Divergence of two independently evolving sequences that originated from a common ancestor can be described by two parameters, the asymptotic level of divergence E and the rate r at which this level of divergence is approached. Constant negative selection impedes allele replacements and, therefore, is routinely assumed to decelerate sequence divergence. However, its impact on E and on r has not been formally investigated. Results Strong selection that favors only one allele can make E arbitrarily small and r arbitrarily large. In contrast, in the case of 4 possible alleles and equal mutation rates, the lowest value of r , attained when two alleles confer equal fitnesses and the other two are strongly deleterious, is only two times lower than its value under selective neutrality. Conclusions Constant selection can strongly constrain the level of sequence divergence, but cannot reduce substantially the rate at which this level is approached. In particular, under any constant selection the divergence of sequences that accumulated one substitution per neutral site since their origin from the common ancestor must already constitute at least one half of the asymptotic divergence at sites under such selection. Reviewers This article was reviewed by Drs. Nicolas Galtier, Sergei Maslov, and Nick Grishin.

Bi-allelic pathogenic variations within POLR3A have been associated with a spectrum of hereditary disorders. Among these, a less frequently observed condition is Wiedemann–Rautenstrauch syndrome (WRS), also known as neonatal progeroid syndrome. This syndrome typically manifests neonatally and is characterized by growth retardation, evident generalized lipodystrophy with distinctively localized fat accumulations, sparse scalp hair, and atypical facial features. Our objective was to elucidate the underlying molecular mechanisms of Wiedemann–Rautenstrauch syndrome (WRS). In this study, we present a clinical case of a 7-year-old female patient diagnosed with WRS. Utilizing whole-exome sequencing (WES), we identified a novel missense variant c.3677T>C (p.Leu1226Pro) in the POLR3A gene (NM_007055.4) alongside two cis intronic variants c.1909+22G>A and c.3337-11T>C. Via the analysis of mRNA derived from fibroblasts, we reconfirmed the splicing-affecting nature of the c.3337-11T>C variant. Furthermore, our investigation led to the reclassification of the c.3677T>C (p.Leu1226Pro) variant as a likely pathogenic variant. Therefore, this is the first case demonstrating the molecular genetics of a patient with Wiedemann–Rautenstrauch syndrome from the Russian Federation. A limited number of clinical cases have been documented until this moment; therefore, broadening the linkage between phenotype and molecular changes in the POLR3A gene will significantly contribute to the comprehensive understanding of the molecular basis of POLR3A-related disorders.