Intracellular Selection, Conversion Bias, and the Expected Substitution Rate of Organelle Genes

A key step in the substitution of a new organelle mutant throughout a population is the generation of germ-line cells homoplasmic for that mutant. Given that each cell typically contains multiple copies of organelles, each of which in turn contains multiple copies of the organelle genome, processes akin to drift and selection in a population are responsible for producing homoplasmic cells. This paper examines the expected substitution rate of new mutants by obtaining the probability that a new mutant is fixed throughout a cell, allowing for arbitrary rates of genome turnover within an organelle and organelle turnover within the cell, as well as (possibly biased) gene conversion and genetic differences in genome and/or organelle replication rates. Analysis is based on a variation of Moran's model for drift in a haploid population. One interesting result is that if the rate of unbiased conversion is sufficiently strong, it creates enough intracellular drift to overcome even strong differences in the replication rates of wild-type and mutant genomes. Thus, organelles with very high conversion rates are more resistant to intracellular selection based on differences in genome replication and/or degradation rates. It is found that the amount of genetic exchange between organelles within the cell greatly influences the probability of fixation. In the absence of exchange, biased gene conversion and/or differences in genome replication rates do not influence the probability of fixation beyond the initial fixation within a single organelle. With exchange, both these processes influence the probability of fixation throughout the entire cell. Generally speaking, exchange between organelles accentuates the effects of directional intracellular forces.(ABSTRACT TRUNCATED AT 250 WORDS)     

Molecular Clock of Neutral Mutations in a Fitness-Increasing Evolutionary Process

The molecular clock of neutral mutations, which represents linear mutation fixation over generations, is theoretically explained by genetic drift in fitness-steady evolution or hitchhiking in adaptive evolution. The present study is the first experimental demonstration for the molecular clock of neutral mutations in a fitness-increasing evolutionary process. The dynamics of genome mutation fixation in the thermal adaptive evolution of Escherichia coli were evaluated in a prolonged evolution experiment in duplicated lineages. The cells from the continuously fitness-increasing evolutionary process were subjected to genome sequencing and analyzed at both the population and single-colony levels. Although the dynamics of genome mutation fixation were complicated by the combination of the stochastic appearance of adaptive mutations and clonal interference, the mutation fixation in the population was simply linear over generations. Each genome in the population accumulated 1.6 synonymous and 3.1 non-synonymous neutral mutations, on average, by the spontaneous mutation accumulation rate, while only a single genome in the population occasionally acquired an adaptive mutation. The neutral mutations that preexisted on the single genome hitchhiked on the domination of the adaptive mutation. The successive fixation processes of the 128 mutations demonstrated that hitchhiking and not genetic drift were responsible for the coincidence of the spontaneous mutation accumulation rate in the genome with the fixation rate of neutral mutations in the population. The molecular clock of neutral mutations to the fitness-increasing evolution suggests that the numerous neutral mutations observed in molecular phylogenetic trees may not always have been fixed in fitness-steady evolution but in adaptive evolution.     

Three new Drosophila markers of intracellular membranes

The need for cellular markers that permit a quick and accurate evaluation of a protein's subcellular localization has increased with the surge of new data generated by the Drosophila genome project. In this report, we present three ubiquitously expressed Drosophila transgenes that expressed a green fluorescent protein variant (enhanced yellow fluorescent protein) that has been targeted to different intracellular membrane targets: the Golgi apparatus, mitochondria, and endoplasmic reticulum. These markers serve as an internal standard for characterizing a protein's subcellular localization or as a means of tracking the dynamics of intracellular organelles during normal or abnormal cellular or developmental processes. We have also examined fixation artifacts using these constructs to illustrate the effects that fixation and permeabilization have on intracellular membrane organization.     

Human endogenous retroviruses in the primate lineage and their influence on host genomes

Primates emerged about 60 million years ago. Since that time various primate-targeting retroviruses have integrated in the germ line of primate species, and some drifted to fixation. After germ line fixation, continued activity of proviruses resulted in intragenomic spread of so-called endogenous retroviruses (ERVs). Variant ERVs emerged, amplified in the genome and profoundly altered genome structures and potentially functionality. Importantly, ERVs are genome modifiers of exogenous origin. The human genome contains about 8% of sequences of retroviral origin. The human ERVs (HERVs) comprise many distinct families that amplified to copy numbers of up to several thousand. We review here the evolution of several well-characterized HERV families in the human lineage since initial germ line fixation. It is apparent that endogenous retroviruses profoundly affected the genomes of species in the evolutionary lineage leading to Homo sapiens.     

A Model of Genome Size Dynamics During Speciation

A model developed for the evolving size of the repetitive part of the eukaryote genome during speciation was subjected to analytical and computer treatment. The basic assumption of the model was that two classes of repetitive DNA contribute mainly to macroevolutionary changes in genome size: arrays of tandem repeats (ATR) changing through unequal crossover and mobile genetic elements (MGE) changing presumably through an integration mechanism of the Tn- and Is-kind operating in bacteria. Within the framework of this model, the macroevolution of the MGE size is formally equivalent to that of the ATR in the particular case when shifts of chromatids have only one repeat out of register. This allowed us to consider genome size as a large set of various ATRs. The results obtained are as follows. If the duplication and deletion of repeats have unequal fixation probabilities during each speciation act, the predicted species distributions of genome size significantly deviate from the real ones; if they have equal fixation probabilities, there is a conformance between calculated and real distributions. In the latter case, the model reproduces the salient features of real distributions upon acceptance of 1) upper selective boundary nonspecifically limiting increase in genome size within the evolving taxonomic group and 2) non-neutrality of variability in genome size with respect to speciation.     

Co-ordination and fine-tuning of nitrogen fixation in Azotobacter vinelandii

Nitrogen fixation by the free-living organism Azotobacter vinelandii can occur through the activity of three different systems that are genetically distinct but mechanistically related. A combination of bioinformatic and biochemical-genetic studies has revealed that at least 82 different genes are likely to be associated with the formation and regulation of these systems. Studies performed over many years have established that cross-talk occurs between the various nitrogen fixation systems, and that expression and fine-tuning of their activities are integrated with overall cellular physiology. Martinez-Noel and co-workers now report another newly discovered aspect of the process. Evidence is presented to suggest that a nitrogen fixation-specific paralogue of ClpX is used to control the accumulation of proteins involved in formation of a metal-sulphur cluster that provides a nitrogenase active site. The intriguing aspect of this work is that it indicates that the nitrogen fixation-associated ClpX must recruit ClpP, for which a paralogue is not duplicated within any of the nitrogen fixation regions of the genome, to achieve its function related to nitrogen fixation. Inspection of the A. vinelandii genome indicates that such recruitment of cellular housekeeping components is a common feature used to integrate nitrogen fixation with global cellular physiology.     

Hybrid genome evolution by transposition

Species hybridization is reviewed focusing on its role as a source of evolutionary novelties. Contrary to the view that hybrids are lineages devoid of evolutionary value, a number of case studies are given that show how hybrids are responsible for reticulate evolution that may lead to the origin of new species. Hybrid evolution is mediated by extensive genome repatterning followed by rapid stabilization and fixation of highly adapted genotypes. Some well-documented cases demonstrate that bursts of transposition follow hybridization and may contribute to the genetic instability observed after hybridization. The mechanism that triggers transposition in hybrids is largely unknown, but coupling of hybrid transposition and demethylation has been observed in mammals and plants. A natural scenario is proposed in which marginal small hybrid populations undergo transposition mediated genome reorganizations accompanied by exogenous and endogenous selection that, in concert with drift, lead to rapid fixation of high fitness hybrid genotypes. These genotypes may represent parental introgressed species or be entirely new species.     

Analysis of human accelerated DNA regions using archaic hominin genomes

Several previous comparisons of the human genome with other primate and vertebrate genomes identified genomic regions that are highly conserved in vertebrate evolution but fast-evolving on the human lineage. These human accelerated regions (HARs) may be regions of past adaptive evolution in humans. Alternatively, they may be the result of non-adaptive processes, such as biased gene conversion. We captured and sequenced DNA from a collection of previously published HARs using DNA from an Iberian Neandertal. Combining these new data with shotgun sequence from the Neandertal and Denisova draft genomes, we determine at least one archaic hominin allele for 84% of all positions within HARs. We find that 8% of HAR substitutions are not observed in the archaic hominins and are thus recent in the sense that the derived allele had not come to fixation in the common ancestor of modern humans and archaic hominins. Further, we find that recent substitutions in HARs tend to have come to fixation faster than substitutions elsewhere in the genome and that substitutions in HARs tend to cluster in time, consistent with an episodic rather than a clock-like process underlying HAR evolution. Our catalog of sequence changes in HARs will help prioritize them for functional studies of genomic elements potentially responsible for modern human adaptations.     

Two genome sequences of the same bacterial strain, Gluconacetobacter diazotrophicus PAl 5, suggest a new standard in genome sequence submission

Gluconacetobacter diazotrophicus PAl 5 is of agricultural significance due to its ability to provide fixed nitrogen to plants. Consequently, its genome sequence has been eagerly anticipated to enhance understanding of endophytic nitrogen fixation. Two groups have sequenced the PAl 5 genome from the same source (ATCC 49037), though the resulting sequences contain a surprisingly high number of differences. Therefore, an optical map of PAl 5 was constructed in order to determine which genome assembly more closely resembles the chromosomal DNA by aligning each sequence against a physical map of the genome. While one sequence aligned very well, over 98% of the second sequence contained numerous rearrangements. The many differences observed between these two genome sequences could be owing to either assembly errors or rapid evolutionary divergence. The extent of the differences derived from sequence assembly errors could be assessed if the raw sequencing reads were provided by both genome centers at the time of genome sequence submission. Hence, a new genome sequence standard is proposed whereby the investigator supplies the raw reads along with the closed sequence so that the community can make more accurate judgments on whether differences observed in a single stain may be of biological origin or are simply caused by differences in genome assembly procedures.     

A Model-Based Analysis of GC-Biased Gene Conversion in the Human and Chimpanzee Genomes

GC-biased gene conversion (gBGC) is a recombination-associated process that favors the fixation of G/C alleles over A/T alleles. In mammals, gBGC is hypothesized to contribute to variation in GC content, rapidly evolving sequences, and the fixation of deleterious mutations, but its prevalence and general functional consequences remain poorly understood. gBGC is difficult to incorporate into models of molecular evolution and so far has primarily been studied using summary statistics from genomic comparisons. Here, we introduce a new probabilistic model that captures the joint effects of natural selection and gBGC on nucleotide substitution patterns, while allowing for correlations along the genome in these effects. We implemented our model in a computer program, called phastBias, that can accurately detect gBGC tracts about 1 kilobase or longer in simulated sequence alignments. When applied to real primate genome sequences, phastBias predicts gBGC tracts that cover roughly 0.3% of the human and chimpanzee genomes and account for 1.2% of human-chimpanzee nucleotide differences. These tracts fall in clusters, particularly in subtelomeric regions; they are enriched for recombination hotspots and fast-evolving sequences; and they display an ongoing fixation preference for G and C alleles. They are also significantly enriched for disease-associated polymorphisms, suggesting that they contribute to the fixation of deleterious alleles. The gBGC tracts provide a unique window into historical recombination processes along the human and chimpanzee lineages. They supply additional evidence of long-term conservation of megabase-scale recombination rates accompanied by rapid turnover of hotspots. Together, these findings shed new light on the evolutionary, functional, and disease implications of gBGC. The phastBias program and our predicted tracts are freely available.     

Genomics insights into symbiotic nitrogen fixation

Following an interaction with rhizobial soil bacteria, legume plants are able to form a novel organ, termed the root nodule. This organ houses the rhizobial microsymbionts, which perform the biological nitrogen fixation process resulting in the incorporation of ammonia into plant organic molecules. Recent advances in genomics have opened exciting new perspectives in this field by providing the complete gene inventory of two rhizobial microsymbionts. The complete genome sequences of Mesorhizobium loti, the symbiont of several Lotus species, and Sinorhizobium meliloti, the symbiont of alfalfa, were determined and annotated in detail. For legume macrosymbionts, expressed sequence tag projects and expression analyses using DNA arrays in conjunction with proteomics approaches have identified numerous genes involved in root nodule formation and nitrogen fixation. The isolation of legume genes by tagging or positional cloning recently allowed the identification of genes that control the very early steps of root nodule organogenesis.     

A neutral theory predicts multigenic aging and increased concentrations of deleterious mutations on the mitochondrial and Y chromosomes

Population genetic forces have molded the constitution of the human genome over evolutionary time, and some of the most important parameters are the initial frequency of the allele, p, the effective population size, Ne, and the selection coefficient, s. There is considerable agreement among evolutionary gerontologists that the amplitude of -s is small for alleles that are Deleterious In Late Life (DILL), and thus DILL traits are effectively neutral and should be fixed in the human population in relationship to Ne and p. Even higher rates of fixation of deleterious mutations are predicted to occur in the two nonrecombinant genomes in humans, i.e., the Y chromosome and the mitochondrial genome, as a consequence of their lower Ne than autosomes, and the predicted higher rate of fixation of deleterious alleles on the Y may explain the reduced average life span of males vs. females. The high probability of fixation of neutral and mildly deleterious mutations in the mitochondrial genome explains in part its fast rate of evolution, the high observed frequency of mitochondrial disease in relationship to this genome's small size, and may be the underlying reason for the transfer of mitochondrial genes over evolutionary time to the nucleus. The predicted higher concentration of deleterious mutations on the mitochondrial genome could have some leverage to cause more dysfunction than that predicted by mitochondrial gene number alone, because of the essential role of mitochondrial gene function in multisubunit complexes, the coupling of mitochondrial functions, the observation that some mtDNA sequences facilitate somatic mutation, and the likelihood of deleterious mutations either increasing the production of or the sensitivity to mitochondrial ROS.     

A new approach for using genome scans to detect recent positive selection in the human genome

Genome-wide scanning for signals of recent positive selection is essential for a comprehensive and systematic understanding of human adaptation. Here, we present a genomic survey of recent local selective sweeps, especially aimed at those nearly or recently completed. A novel approach was developed for such signals, based on contrasting the extended haplotype homozygosity (EHH) profiles between populations. We applied this method to the genome single nucleotide polymorphism (SNP) data of both the International HapMap Project and Perlegen Sciences, and detected widespread signals of recent local selection across the genome, consisting of both complete and partial sweeps. A challenging problem of genomic scans of recent positive selection is to clearly distinguish selection from neutral effects, given the high sensitivity of the test statistics to departures from neutral demographic assumptions and the lack of a single, accurate neutral model of human history. We therefore developed a new procedure that is robust across a wide range of demographic and ascertainment models, one that indicates that certain portions of the genome clearly depart from neutrality. Simulations of positive selection showed that our tests have high power towards strong selection sweeps that have undergone fixation. Gene ontology analysis of the candidate regions revealed several new functional groups that might help explain some important interpopulation differences in phenotypic traits.     

Genome sequence of Rhizobium etli CNPAF512, a nitrogen-fixing symbiont isolated from bean root nodules in Brazil

Rhizobium etli is a Gram-negative soil-dwelling alphaproteobacterium that carries out symbiotic biological nitrogen fixation in close association with legume hosts. R. etli strains exhibit high sequence divergence and are geographically structured, with a potentially dramatic influence on the outcome of symbiosis. Here, we present the genome sequence of R. etli CNPAF512, a Brazilian isolate from bean nodules. We anticipate that the availability of genome sequences of R. etli strains from distinctly different areas will provide valuable new insights into the geographic mosaic of the R. etli pangenome and the evolutionary dynamics that shape it.     

Life in hot acid: pathway analyses in extremely thermoacidophilic archaea

The extremely thermoacidophilic archaea are a particularly intriguing group of microorganisms that must simultaneously cope with biologically extreme pHs (< or = 4) and temperatures (Topt > or = 60 degrees C) in their natural environments. Their expanding biotechnological significance relates to their role in biomining of base and precious metals and their unique mechanisms of survival in hot acid, at both the cellular and biomolecular levels. Recent developments, such as advances in understanding of heavy metal tolerance mechanisms, implementation of a genetic system, and discovery of a new carbon fixation pathway, have been facilitated by the availability of genome sequence data and molecular genetic systems. As a result, new insights into the metabolic pathways and physiological features that define extreme thermoacidophily have been obtained, in some cases suggesting prospects for biotechnological opportunities.     

Benchmarking of single-virus genomics: a new tool for uncovering the virosphere

Metagenomics and single-cell genomics have enabled the discovery of relevant uncultured microbes. Recently, single-virus genomics (SVG), although still in an incipient stage, has opened new avenues in viral ecology by allowing the sequencing of one single virus at a time. The investigation of methodological alternatives and optimization of existing procedures for SVG is paramount to deliver high-quality genomic data. We report a sequencing dataset of viral single-amplified genomes (vSAGs) from cultured and uncultured viruses obtained by applying different conditions in each SVG step, from viral preservation and novel whole-genome amplification (WGA) to sequencing platforms and genome assembly. Sequencing data showed that cryopreservation and mild fixation were compatible with WGA, although fresh samples delivered better genome quality data. The novel TruPrime WGA, based on primase-polymerase features, and WGA-X employing a thermostable phi29 polymerase, were proven to be with sufficient sensitivity in SVG. The Oxford Nanopore (ON) sequencing platform did not provide a significant improvement of vSAG assembly compared to Illumina alone. Finally, the SPAdes assembler performed the best. Overall, our results represent a valuable genomic dataset that will help to standardized and advance new tools in viral ecology.     

Recombination and hitchhiking of deleterious alleles

When new advantageous alleles arise and spread within a population, deleterious alleles at neighboring loci can hitchhike alongside them and spread to fixation in areas of low recombination, introducing a fixed mutation load. We use branching processes and diffusion equations to calculate the probability that a deleterious allele hitchhikes and fixes alongside an advantageous mutant. As expected, the probability of fixation of a deleterious hitchhiker rises with the selective advantage of the sweeping allele and declines with the selective disadvantage of the deleterious hitchhiker. We then use computer simulations of a genome with an infinite number of loci to investigate the increase in load after an advantageous mutant is introduced. We show that the appearance of advantageous alleles on genetic backgrounds loaded with deleterious alleles has two potential effects: it can fix deleterious alleles, and it can facilitate the persistence of recombinant lineages that happen to occur. The latter is expected to reduce the signals of selection in the surrounding region. We consider these results in light of human genetic data to infer how likely it is that such deleterious hitchhikers have occurred in our recent evolutionary past.     

Pseudogene evolution and natural selection for a compact genome

Pseudogenes are nonfunctional copies of protein-coding genes that are presumed to evolve without selective constraints on their coding function. They are of considerable utility in evolutionary genetics because, in the absence of selection, different types of mutations in pseudogenes should have equal probabilities of fixation. This theoretical inference justifies the estimation of patterns of spontaneous mutation from the analysis of patterns of substitutions in pseudogenes. Although it is possible to test whether pseudogene sequences evolve without constraints for their protein-coding function, it is much more difficult to ascertain whether pseudogenes may affect fitness in ways unrelated to their nucleotide sequence. Consider the possibility that a pseudogene affects fitness merely by increasing genome size. If a larger genome is deleterious--for example, because of increased energetic costs associated with genome replication and maintenance--then deletions, which decrease the length of a pseudogene, should be selectively advantageous relative to insertions or nucleotide substitutions. In this article we examine the implications of selection for genome size relative to small (1-400 bp) deletions, in light of empirical evidence pertaining to the size distribution of deletions observed in Drosophila and mammalian pseudogenes. There is a large difference in the deletion spectra between these organisms. We argue that this difference cannot easily be attributed to selection for overall genome size, since the magnitude of selection is unlikely to be strong enough to significantly affect the probability of fixation of small deletions in Drosophila.