Positive and purifying selection in mitochondrial genomes of a bird with mitonuclear discordance

Diversifying selection on metabolic pathways can reduce intraspecific gene flow and promote population divergence. An opportunity to explore this arises from mitonuclear discordance observed in an Australian bird Eopsaltria australis. Across >1500 km, nuclear differentiation is low and latitudinally structured by isolation by distance, whereas two highly divergent, parapatric mitochondrial lineages (>6.6% in ND2) show a discordant longitudinal geographic pattern and experience different climates. Vicariance, incomplete lineage sorting and sex‐biased dispersal were shown earlier to be unlikely drivers of the mitonuclear discordance; instead, natural selection on a female‐linked trait was the preferred hypothesis. Accordingly, here we tested for signals of positive, divergent selection on mitochondrial genes in E. australis. We used codon models and physicochemical profiles of amino acid replacements to analyse complete mitochondrial genomes of the two mitochondrial lineages in E. australis, its sister species Eopsaltria griseogularis, and outgroups. We found evidence of positive selection on at least five amino acids, encoded by genes of two oxidative phosphorylation pathway complexes NADH dehydrogenase (ND4 and ND4L) and cytochrome bc₁ (cyt‐b) against a background of widespread purifying selection on all mitochondrial genes. Three of these amino acid replacements were fixed in ND4 of the geographically most widespread E. australis lineage. The other two replacements were fixed in ND4L and cyt‐b of the geographically more restricted E. australis lineage. We discuss whether this selection may reflect local environmental adaptation, a by‐product of other selective processes, or genetic incompatibilities, and propose how these hypotheses can be tested in future.     

The evolution of mitochondrial genomes in subterranean caviomorph rodents: adaptation against a background of purifying selection

South American tuco-tucos (Ctenomys) and the related coruro (Spalacopus) are two rodent lineages that have independently colonised the subterranean niche. The energetically demanding lifestyles of these species, coupled with the hypoxic atmospheres characteristic of subterranean environments, may have altered the selective regimes on genes encoding proteins related to cellular respiration. Here, we examined the molecular evolution of 13 protein-coding genes in the mitochondrial genome of seven caviomorph rodents, including these two subterranean genera and their above-ground relatives. Using maximum-likelihood and Bayesian approaches, we estimated rates of synonymous (dS) and nonsynonymous (dN) substitutions. We found a significantly higher ω ratio (dN/dS) in subterranean groups as compared to their non-subterranean counterparts in 11 of 13 genes, although no ω ratio was larger than 1. Additionally, we applied a method based on quantitative physicochemical properties to test for positive selection. Amino acid changes implicated in radical structural or functional shifts in the protein property were found to be ubiquitous across the phylogeny, but concentrated in the subterranean lineages. Convergent changes were also found between the subterranean genera used in this study and other mammals adapted to hypoxia. The results of this study suggest a link between niche shifts and weak directional (or episodic) selection at the molecular level against a background of purifying selection.     

Ultrastructural organization of recombinant Marburg virus nucleoprotein: comparison with Marburg virus inclusions

HeLa cells expressing the recombinant Marburg virus (MBGV) nucleoprotein (NP) have been studied by immunoelectron microscopy. It was found that MBGV NPs assembled into large aggregates which were in close association with membranes of the rough endoplasmic reticulum. Further analysis of these aggregates revealed that NPs formed tubule-like structures which were arranged in a hexagonal pattern. A similar pattern of preformed nucleocapsids was detected in intracellular inclusions induced by MBGV infection. Our data indicated that MBGV NP is able to form nucleocapsid-like structures in the absence of the authentic viral genome and other nucleocapsid-associated proteins.     

Genealogies and weak purifying selection

The assumption that selection alters the genealogical tree of a sample of alleles from a population relative to the neutral expectation underlies several "tests of neutrality." Two recent papers have studied the effect of purifying selection; their suggestive but incomplete results indicate that, in the single site case, the shape of a gene genealogy for a locus may differ only from the neutral expectation. We verify this finding for weak selection using the "ancestral selection graph." We consider a wider range of models, including both a four-allele single-site model and an infinite-sites model. Our results confirm the previous claim for the symmetric-mutation single site model. We emphasize, however, that a neutral-seeming genealogy is consistent with detectable effects of selection on the distribution of allele frequences within the sample. With selection operating, the information about a sample cannot be reduced to the genealogy. As a result, a distinction needs to be made between the selected sites themselves, for which the genealogy offers insufficient information, and linked neutral variation. This distinction seems to have been overlooked in previous papers, yet it has significant implications for the interpretation of data on DNA sequence variation. In particular, it predicts that under purifying selection, the frequency spectrum of neutral mutations will not reflect the skew toward rare polymorphisms at replacement sites even if there is no recombination between them. We caution, however, that the effect of weak selection on the genealogy is specific to the model; a (more realistic) model of multiple linked sites could lead to a more distorted genealogy than is observed for a single site.     

Patterns of intra- and interhost nonsynonymous variation reveal strong purifying selection in dengue virus

Considerable uncertainty surrounds the evolutionary rates of and selection pressures acting on arthropod-borne RNA viruses (arboviruses). In particular, it is unclear why arboviruses such as dengue virus show substantial genetic variation within individual humans and mosquitoes yet low long-term rates of amino acid substitution. To address this question, I compared patterns of nonsynonymous variation in populations of dengue virus sampled at different levels of evolutionary divergence. Although nonsynonymous variation was abundant in viral populations within individual humans, there was a marked reduction in the frequency of nonsynonymous mutations in interhost comparisons. Moreover, intrahost genetic variation corresponded to a random pattern of mutation, and most of the sites that exhibited nonsynonymous variation within hosts were invariant at deeper phylogenetic levels. This loss of long-term nonsynonymous variation is the signature of extensive purifying selection such that more than 90% of all nonsynonymous mutations are deleterious. Consequently, although arboviruses are able to successfully adapt to diverse cell types, they are characterized by a high rate of deleterious mutation.     

Evolution of the human immunodeficiency virus envelope gene is dominated by purifying selection

The evolution of the human immunodeficiency virus (HIV-1) during chronic infection involves the rapid, continuous turnover of genetic diversity. However, the role of natural selection, relative to random genetic drift, in governing this process is unclear. We tested a stochastic model of genetic drift using partial envelope sequences sampled longitudinally in 28 infected children. In each case the Bayesian posterior (empirical) distribution of coalescent genealogies was estimated using Markov chain Monte Carlo methods. Posterior predictive simulation was then used to generate a null distribution of genealogies assuming neutrality, with the null and empirical distributions compared using four genealogy-based summary statistics sensitive to nonneutral evolution. Because both null and empirical distributions were generated within a coalescent framework, we were able to explicitly account for the confounding influence of demography. From the distribution of corrected P-values across patients, we conclude that empirical genealogies are more asymmetric than expected if evolution is driven by mutation and genetic drift only, with an excess of low-frequency polymorphisms in the population. This indicates that although drift may still play an important role, natural selection has a strong influence on the evolution of HIV-1 envelope. A negative relationship between effective population size and substitution rate indicates that as the efficacy of selection increases, a smaller proportion of mutations approach fixation in the population. This suggests the presence of deleterious mutations. We therefore conclude that intrahost HIV-1 evolution in envelope is dominated by purifying selection against low-frequency deleterious mutations that do not reach fixation.     

Outbreake of Marburg virus disease in Johannesburg.

The first recognised outbreak of Marburg virus disease in Africa, and the first since the original epidemic in West Germany and Yugoslavia in 1967, occurred in South Africa in February 1975. The primary case was in a young Australian man , who was admitted to the Johannesburg Hospital after having toured Rhodesia. Two secondary cases occurred, one being in the first patient''s travelling companion, and the other in a nurse. Features of the illness included high fever, myalgia, vomiting and diarrhoea, hepatitis, a characteristic maculopapular rash, leucopenia, thrombocytopenia, and a bleeding tendency. The first patient died on the seventh day from haemorrhage resulting from a combination of disseminated intravascular coagulation and hepatic failure. The other two patients were given vigorous supportive treatment and prophylactic heparin and recovered after an acute phase lasting about seven days. During this period on developed pancreatitis, the serum amylase remaining raised until the 32nd day after the onset of the illness. The other developed unilateral uveitis after having been asymptomatic for two months. This persisted for several weeks and Marburg virus was cultured from the anterior chamber of the eye.     

Detecting amino acid sites under positive selection and purifying selection

An excess of nonsynonymous over synonymous substitution at individual amino acid sites is an important indicator that positive selection has affected the evolution of a protein between the extant sequences under study and their most recent common ancestor. Several methods exist to detect the presence, and sometimes location, of positively selected sites in alignments of protein-coding sequences. This article describes the "sitewise likelihood-ratio" (SLR) method for detecting nonneutral evolution, a statistical test that can identify sites that are unusually conserved as well as those that are unusually variable. We show that the SLR method can be more powerful than currently published methods for detecting the location of positive selection, especially in difficult cases where the strength of selection is low. The increase in power is achieved while relaxing assumptions about how the strength of selection varies over sites and without elevated rates of false-positive results that have been reported with some other methods. We also show that the SLR method performs well even under circumstances where the results from some previous methods can be misleading.     

Conditional gene genealogies under strong purifying selection

The ancestral selection graph, conditioned on the allelic types in the sample, is used to obtain a limiting gene genealogical process under strong selection. In an equilibrium, two-allele system with strong selection, neutral gene genealogies are predicted for random samples and for samples containing at most one unfavorable allele. Samples containing more than one unfavorable allele have gene genealogies that differ greatly from neutral predictions. However, they are related to neutral gene genealogies via the well-known Ewens sampling formula. Simulations show rapid convergence to limiting analytical predictions as the strength of selection increases. These results extend the idea of a soft selective sweep to deleterious alleles and have implications for the interpretation of polymorphism among disease-causing alleles in humans.     

Strong purifying selection in transmission of mammalian mitochondrial DNA

There is an intense debate concerning whether selection or demographics has been most important in shaping the sequence variation observed in modern human mitochondrial DNA (mtDNA). Purifying selection is thought to be important in shaping mtDNA sequence evolution, but the strength of this selection has been debated, mainly due to the threshold effect of pathogenic mtDNA mutations and an observed excess of new mtDNA mutations in human population data. We experimentally addressed this issue by studying the maternal transmission of random mtDNA mutations in mtDNA mutator mice expressing a proofreading-deficient mitochondrial DNA polymerase. We report a rapid and strong elimination of nonsynonymous changes in protein-coding genes; the hallmark of purifying selection. There are striking similarities between the mutational patterns in our experimental mouse system and human mtDNA polymorphisms. These data show strong purifying selection against mutations within mtDNA protein-coding genes. To our knowledge, our study presents the first direct experimental observations of the fate of random mtDNA mutations in the mammalian germ line and demonstrates the importance of purifying selection in shaping mitochondrial sequence diversity.     

Carbohydrate structure of Marburg virus glycoprotein.

Marburg virus was propagated in E6 cells, a cloned cell line of Vero cells, in the presence of [6-3H]glucosamine. Radiolabelled viral glycoprotein was digested with trypsin, and oligosaccharides were liberated by sequential treatment with endo-beta-N-acetylglucosaminidase H, peptide-N4-(N-acetyl-beta-glucosaminyl)asparagine amidase F and O-glycosidase, by beta-elimination, and by alkaline hydrolysis. After fractionation by HPLC and gel filtration, glycans were characterized chromatographically, by digestion with exoglycosidases and, in part, by methylation analysis and liquid secondary ion mass spectrometry. The oligosaccharide structures thus established include oligomannosidic and hybrid-type N-glycans, as well as neutral fucosylated bi-, tri- and tetraantennary species, most of which carry an additional bisecting N-acetylglucosamine. In addition, high amounts of neutral mucin-type O-glycans with type-1 and type-2 core structures were detected. None of the glycans present in this viral glycoprotein carried sialic acid residues.     

Annotation of selection strengths in viral genomes

MOTIVATION: Viral genomes tend to code in overlapping reading frames to maximize informational content. This may result in atypical codon bias and particular evolutionary constraints. Due to the fast mutation rate of viruses, there is additional strong evidence for varying selection between intra- and intergenomic regions. The presence of multiple coding regions complicates the concept of K(a)/K(s) ratio, and thus begs for an alternative approach when investigating selection strengths. Building on the paper by McCauley and Hein, we develop a method for annotating a viral genome coding in overlapping reading frames. We introduce an evolutionary model capable of accounting for varying levels of selection along the genome, and incorporate it into our prior single sequence HMM methodology, extending it now to a phylogenetic HMM. Given an alignment of several homologous viruses to a reference sequence, we may thus achieve an annotation both of coding regions as well as selection strengths, allowing us to investigate different selection patterns and hypotheses. RESULTS: We illustrate our method by applying it to a multiple alignment of four HIV2 sequences, as well as of three Hepatitis B sequences. We obtain an annotation of the coding regions, as well as a posterior probability for each site of the strength of selection acting on it. From this we may deduce the average posterior selection acting on the different genes. Whilst we are encouraged to see in HIV2, that the known to be conserved genes gag and pol are indeed annotated as such, we also discover several sites of less stringent negative selection within the env gene. To the best of our knowledge, we are the first to subsequently provide a full selection annotation of the Hepatitis B genome by explicitly modelling the evolution within overlapping reading frames, and not relying on simple K(a)/K(s) ratios.     

De novo mutational signature discovery in tumor genomes using SparseSignatures

Cancer is the result of mutagenic processes that can be inferred from tumor genomes by analyzing rate spectra of point mutations, or “mutational signatures”. Here we present SparseSignatures, a novel framework to extract signatures from somatic point mutation data. Our approach incorporates a user-specified background signature, employs regularization to reduce noise in non-background signatures, uses cross-validation to identify the number of signatures, and is scalable to large datasets. We show that SparseSignatures outperforms current state-of-the-art methods on simulated data using a variety of standard metrics. We then apply SparseSignatures to whole genome sequences of pancreatic and breast tumors, discovering well-differentiated signatures that are linked to known mutagenic mechanisms and are strongly associated with patient clinical features.