Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution

An excess of nonsynonymous substitutions over synonymous ones is an important indicator of positive selection at the molecular level. A lineage that underwent Darwinian selection may have a nonsynonymous/synonymous rate ratio (dN/dS) that is different from those of other lineages or greater than one. In this paper, several codon-based likelihood models that allow for variable dN/dS ratios among lineages were developed. They were then used to construct likelihood ratio tests to examine whether the dN/dS ratio is variable among evolutionary lineages, whether the ratio for a few lineages of interest is different from the background ratio for other lineages in the phylogeny, and whether the dN/dS ratio for the lineages of interest is greater than one. The tests were applied to the lysozyme genes of 24 primate species. The dN/dS ratios were found to differ significantly among lineages, indicating that the evolution of primate lysozymes is episodic, which is incompatible with the neutral theory. Maximum- likelihood estimates of parameters suggested that about nine nonsynonymous and zero synonymous nucleotide substitutions occurred in the lineage leading to hominoids, and the dN/dS ratio for that lineage is significantly greater than one. The corresponding estimates for the lineage ancestral to colobine monkeys were nine and one, and the dN/dS ratio for the lineage is not significantly greater than one, although it is significantly higher than the background ratio. The likelihood analysis thus confirmed most, but not all, conclusions Messier and Stewart reached using reconstructed ancestral sequences to estimate synonymous and nonsynonymous rates for different lineages.      

Independent effects of alternative splicing and structural constraint on the evolution of mammalian coding exons

Alternative splicing (AS) is known to significantly affect exon-level protein evolutionary rates in mammals. Particularly, alternatively spliced exons (ASEs) have a higher nonsynonymous-to-synonymous substitution rate (dN/dS) ratio than constitutively spliced exons (CSEs), possibly because the former are required only occasionally for normal biological functions. Meanwhile, intrinsically disordered regions (IDRs), the protein regions lacking fixed 3D structures, are also reported to have an increased evolutionary rate due to lack of structural constraint. Interestingly, IDRs tend to be located in alternative protein regions. Yet which of these two factors is the major determinant of the increased dN/dS in mammalian ASEs remains unclear. By comparing human-macaque and human-mouse one-to-one orthologous genes, we demonstrate that AS and protein structural disorder have independent effects on mammalian exon evolution. We performed analyses of covariance to demonstrate that the slopes of the (dN/dS-percentage of IDR) regression lines differ significantly between CSEs and ASEs. In other words, the dN/dS ratios of both ASEs and CSEs increase with the proportion of IDR (PIDR), whereas ASEs have higher dN/dS ratios than CSEs when they have similar PIDRs. Since ASEs and IDRs may less frequently overlap with protein domains (which also affect dN/dS), we also examined the correlations between dN/dS ratio and exon type/PIDR by controlling for the density of protein domain. We found that the effects of exon type and PIDR on dN/dS are both independent of domain density. Our results imply that nature can select for different biological features with regard to ASEs and IDRs, even though the two biological features tend to be localized in the same protein regions.     

Mitochondrial DNA as a tool for reconstructing past life‐history traits in mammals

Reconstructing the ancestral characteristics of species is a major goal in evolutionary and comparative biology. Unfortunately, fossils are not always available and sufficiently informative, and phylogenetic methods based on models of character evolution can be unsatisfactory. Genomic data offer a new opportunity to estimate ancestral character states, through (i) the correlation between DNA evolutionary processes and species life‐history traits and (ii) available reliable methods for ancestral sequence inference. Here, we assess the relevance of mitochondrial DNA – the most popular molecular marker in animals – as a predictor of ancestral life‐history traits in mammals, using the order of Cetartiodactyla as a benchmark. Using the complete set of 13 mitochondrial protein‐coding genes, we show that the lineage‐specific nonsynonymous over synonymous substitution rate ratio (dN/dS) is closely correlated with the species body mass, longevity and age of sexual maturity in Cetartiodactyla and can be used as a marker of ancestral traits provided that the noise introduced by short branches is appropriately dealt with. Based on ancestral dN/dS estimates, we predict that the first cetartiodactyls were relatively small animals (around 20 kg). This finding is in accordance with Cope's rule and the fossil record but could not be recovered via continuous character evolution methods.     

Differential rates of evolution for the ZFY-related zinc finger genes,Zfy, Zfx,and Zfa in the mouse genus Mus

A comparative study of the last exon of the zinc finger genes Zfx, Zfy, and Zfa from species of mice in the genus Mus was conducted to assess the extent of gene-specific and chromosome-specific effects on the evolutionary patterns among related X-, Y-, and autosomal-linked genes. Phylogenetic analyses of 29 sequences from Zfx, Zfa, and Zfy from 10 taxa were performed to infer relatedness among the zinc finger loci, and codon-based maximum likelihood analyses were conducted to assess evolutionary pattern among genes. Five models of nucleotide sequence evolution were applied and compared using a likelihood ratio test. Estimates of nonsynonymous to synonymous changes (dN/dS) for these genes suggest that amino acid substitutions are occurring at a more rapid rate across the autosomal- and Y-specific lineages compared to the X-specific lineage, with the Y-specific lineage showing the highest rate under certain models. The data suggest the action of gene-specific effects on evolutionary pattern. In particular, Zfa and Zfy genes, both with presumed restricted expression, appear less functionally constrained relative to ubiquitously expressed Zfx. Slightly elevated dN/dS for Zfy genes in comparison to Zfa also suggest Y-specific effects.     

Lower prevalence but similar fitness in a parasitic fungus at higher radiation levels near Chernobyl

Nuclear disasters at Chernobyl and Fukushima provide examples of effects of acute ionizing radiation on mutations that can affect the fitness and distribution of species. Here, we investigated the prevalence of Microbotryum lychnidis‐dioicae, a pollinator‐transmitted fungal pathogen of plants causing anther‐smut disease in Chernobyl, its viability, fertility and karyotype variation, and the accumulation of nonsynonymous mutations in its genome. We collected diseased flowers of Silene latifolia from locations ranging by more than two orders of magnitude in background radiation, from 0.05 to 21.03 μGy/h. Disease prevalence decreased significantly with increasing radiation level, possibly due to lower pollinator abundance and altered pollinator behaviour. Viability and fertility, measured as the budding rate of haploid sporidia following meiosis from the diploid teliospores, did not vary with increasing radiation levels and neither did karyotype overall structure and level of chromosomal size heterozygosity. We sequenced the genomes of twelve samples from Chernobyl and of four samples collected from uncontaminated areas and analysed alignments of 6068 predicted genes, corresponding to 1.04 × 10⁷ base pairs. We found no dose‐dependent differences in substitution rates (neither dN, dS, nor dN/dS). Thus, we found no significant evidence of increased deleterious mutation rates at higher levels of background radiation in this plant pathogen. We even found lower levels of nonsynonymous substitution rates in contaminated areas compared to control regions, suggesting that purifying selection was stronger in contaminated than uncontaminated areas. We briefly discuss the possibilities for a mechanistic basis of radio resistance in this nonmelanized fungus.     

Cloning and characterization of wheat PDI (protein disulfide isomerase) homoeologous genes and promoter sequences

The genomic and cDNA sequences of three PDI homoeologous genes located on chromosomes 4A, 4B and 4D of bread wheat and their promoters were cloned and sequenced. The three sequences showed a very high conservation of the coding region and of the exon/intron structure, which consisted of ten exons. The comparison of wheat sequences with those of rice and Arabidopsis showed a significant conservation of the exon/intron structure across the three species. The expression of each gene was analysed by RT-PCR in different plant tissues (roots, coleoptiles, spikelets, leaves and developing caryopses). All the genes showed a higher expression in developing caryopses than in other analysed tissues, wherein some differences were detected. The promoter sequences of the three genes possessed some regulatory motifs typical of endosperm specific expression.     

Differentiation of Xylella fastidiosa strains via multilocus sequence analysis of environmentally mediated genes (MLSA-E)

Isolates of the plant pathogen Xylella fastidiosa are genetically very similar, but studies on their biological traits have indicated differences in virulence and infection symptomatology. Taxonomic analyses have identified several subspecies, and phylogenetic analyses of housekeeping genes have shown broad host-based genetic differences; however, results are still inconclusive for genetic differentiation of isolates within subspecies. This study employs multilocus sequence analysis of environmentally mediated genes (MLSA-E; genes influenced by environmental factors) to investigate X. fastidiosa relationships and differentiate isolates with low genetic variability. Potential environmentally mediated genes, including host colonization and survival genes related to infection establishment, were identified a priori. The ratio of the rate of nonsynonymous substitutions to the rate of synonymous substitutions (dN/dS) was calculated to select genes that may be under increased positive selection compared to previously studied housekeeping genes. Nine genes were sequenced from 54 X. fastidiosa isolates infecting different host plants across the United States. Results of maximum likelihood (ML) and Bayesian phylogenetic (BP) analyses are in agreement with known X. fastidiosa subspecies clades but show novel within-subspecies differentiation, including geographic differentiation, and provide additional information regarding host-based isolate variation and specificity. dN/dS ratios of environmentally mediated genes, though <1 due to high sequence similarity, are significantly greater than housekeeping gene dN/dS ratios and correlate with increased sequence variability. MLSA-E can more precisely resolve relationships between closely related bacterial strains with low genetic variability, such as X. fastidiosa isolates. Discovering the genetic relationships between X. fastidiosa isolates will provide new insights into the epidemiology of populations of X. fastidiosa, allowing improved disease management in economically important crops.     

Selective Pressures on <Emphasis Type="Italic">Drosophila</Emphasis> Chemosensory Receptor Genes

The evolution and patterns of selection of genes encoding 10 Drosophila odorant receptors (Or) and the sex pheromone receptor Gr68a were investigated by comparing orthologous sequences across five to eight ecologically diverse species of Drosophila. Using maximum likelihood estimates of dN/dS ratios we show that all 11 genes sampled are under purifying selection, indicating functional constraint. Four of these genes (Or33c, Or42a, Or85e, and Gr68a) may be under positive selection, and if so, there is good evidence that 12 specific amino acid sites may be under positive selection. All of these sites are predicted to be located either in loop regions or just inside membrane spanning regions, and interestingly one of the two sites in Gr68a is in a similar position to a previously described polymorphism in Gr5a that causes a shift in sensitivity to its ligand trehalose. For three Ors, possible evidence for positive selection was detected along a lineage. These include Or22a in the lineage leading to D. mauritiana and Or22b in the lineage leading to D. simulans. This is of interest in light of previous data showing a change in ligand response profile for these species in the sensory neuron (ab3A) which expresses both Or22a and Or22b in D. melanogaster. In summary, while the main chemosensory function and/or structural integrity of these 10 Or genes and Gr68a are evolutionarily preserved, positive selection appears to be acting on some of these genes, at specific sites and along certain lineages, and provides testable hypotheses for further functional experimentation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. [Reviewing Editor: Dr. David Pollock]     

Evolution of MHC class II SLA‐DRB1 locus in the Croatian wild boar (Sus scrofa) implies duplication and weak signals of positive selection

The wild boar is an ancestor of the domestic pig and an important game species with the widest geographical range of all ungulates. Although a large amount of data are available on major histocompatibility complex (MHC) variability in domestic pigs, only a few studies have been performed on wild boars. Due to their crucial role in appropriate immune responses and extreme polymorphism, MHC genes represent some of the best candidates for studying the processes of adaptive evolution. Here, we present the results on the variability and evolution of the entire MHC class II SLA‐DRB1 locus exon 2 in 133 wild boars from Croatia. Using direct sequencing and cloning methods, we identified 20 SLA‐DRB1 alleles, including eight new variants, with notable divergence. In some individuals, we documented functional locus duplication, and SLA‐DRB1*04:10 was identified as the allele involved in the duplication. The expression of a duplicated locus was confirmed by cloning and sequencing cDNA‐derived amplicons. Based on individual genotypes, we were able to assume that alleles SLA‐DRB1*04:10 and SLA‐DRB1*06:07 are linked as an allelic combination that co‐evolves as a two‐locus haplotype. Our investigation of evolutionary processes at the SLA‐DRB1 locus confirmed the role of intralocus recombination in generating allelic variability, whereas tests of positive selection based on the dN/dS (non‐synonymous/synonymous substitution rate ratio) test revealed atypically weak and ambiguous signals.     

Gene-specific signatures of elevated non-synonymous substitution rates correlate poorly across the Plasmodium genus

Background

Comparative genome analyses of parasites allow large scale investigation of selective pressures shaping their evolution. An acute limitation to such analysis of Plasmodium falciparum is that there is only very partial low-coverage genome sequence of the most closely related species, the chimpanzee parasite P. reichenowi. However, if orthologous genes have been under similar selective pressures throughout the Plasmodium genus then positive selection on the P. falciparum lineage might be predicted to some extent by analysis of other lineages.

Principal Findings

Here, three independent pairs of closely related species in different sub-generic clades (P. falciparum and P. reichenowi; P. vivax and P. knowlesi; P. yoelii and P. berghei) were compared for a set of 43 candidate ligand genes considered likely to be under positive directional selection and a set of 102 control genes for which there was no selective hypothesis. The ratios of non-synonymous to synonymous substitutions (dN/dS) were significantly elevated in the candidate ligand genes compared to control genes in each of the three clades. However, the rank order correlation of dN/dS ratios for individual candidate genes was very low, less than the correlation for the control genes.

Significance

The inability to predict positive selection on a gene in one lineage by identifying elevated dN/dS ratios in the orthologue within another lineage needs to be noted, as it reflects that adaptive mutations are generally rare events that lead to fixation in individual lineages. Thus it is essential to complete the genome sequences of particular species of phylogenetic importance, such as P. reichenowi.     

Selection in coastal Synechococcus (cyanobacteria) populations evaluated from environmental metagenomes

Environmental metagenomics provides snippets of genomic sequences from all organisms in an environmental sample and are an unprecedented resource of information for investigating microbial population genetics. Current analytical methods, however, are poorly equipped to handle metagenomic data, particularly of short, unlinked sequences. A custom analytical pipeline was developed to calculate dN/dS ratios, a common metric to evaluate the role of selection in the evolution of a gene, from environmental metagenomes sequenced using 454 technology of flow-sorted populations of marine Synechococcus, the dominant cyanobacteria in coastal environments. The large majority of genes (98%) have evolved under purifying selection (dN/dS<1). The metagenome sequence coverage of the reference genomes was not uniform and genes that were highly represented in the environment (i.e. high read coverage) tended to be more evolutionarily conserved. Of the genes that may have evolved under positive selection (dN/dS>1), 77 out of 83 (93%) were hypothetical. Notable among annotated genes, ribosomal protein L35 appears to be under positive selection in one Synechococcus population. Other annotated genes, in particular a possible porin, a large-conductance mechanosensitive channel, an ATP binding component of an ABC transporter, and a homologue of a pilus retraction protein had regions of the gene with elevated dN/dS. With the increasing use of next-generation sequencing in metagenomic investigations of microbial diversity and ecology, analytical methods need to accommodate the peculiarities of these data streams. By developing a means to analyze population diversity data from these environmental metagenomes, we have provided the first insight into the role of selection in the evolution of Synechococcus, a globally significant primary producer.     

Dynamic virulence‐related regions of the plant pathogenic fungus Verticillium dahliae display enhanced sequence conservation

Plant pathogens continuously evolve to evade host immune responses. During host colonization, many fungal pathogens secrete effectors to perturb such responses, but these in turn may become recognized by host immune receptors. To facilitate the evolution of effector repertoires, such as the elimination of recognized effectors, effector genes often reside in genomic regions that display increased plasticity, a phenomenon that is captured in the two‐speed genome hypothesis. The genome of the vascular wilt fungus Verticillium dahliae displays regions with extensive presence/absence polymorphisms, so‐called lineage‐specific regions, that are enriched in in planta‐induced putative effector genes. As expected, comparative genomics reveals differential degrees of sequence divergence between lineage‐specific regions and the core genome. Unanticipated, lineage‐specific regions display markedly higher sequence conservation in coding as well as noncoding regions than the core genome. We provide evidence that disqualifies horizontal transfer to explain the observed sequence conservation and conclude that sequence divergence occurs at a slower pace in lineage‐specific regions of the V. dahliae genome. We hypothesize that differences in chromatin organisation may explain lower nucleotide substitution rates in the plastic, lineage‐specific regions of V. dahliae.     

Preponderance of slightly deleterious polymorphism in mitochondrial DNA: nonsynonymous/synonymous rate ratio is much higher within species than between species

We estimated synonymous (dN) and nonsynonymous (dS) substitution rates for protein-coding genes of the mitochondrial genome from two individuals each of the species human, chimpanzee, and gorilla. The genes were analyzed both separately and in a combined data set. Pairwise sequence comparisons suggest that the dN/dS rate ratios are about 5-10 times higher in within-species comparisons than in between-species comparisons. This result is confirmed by a more rigorous likelihood ratio test, which rejected the null hypothesis that the dN/dS rate ratios are identical within and between species. The likelihood models account for the genetic code structure, transition/transversion rate ratio, and codon usage bias and are expected to produce more reliable results than the commonly used contingency test. Separate analyses of different genes show that the dN/dS rate ratios are higher within species than between species for all 13 mitochondrial genes, with the difference being statistically significant for all except three small or slowly evolving genes. Furthermore, in conserved genes, nonsynonymous rates within species tend to be higher than the between-species rates by a greater proportion than in fast-changing genes. Our findings confirm and extend earlier results obtained from smaller data sets and suggest the operation of slightly deleterious mutations throughout the mitochondrial genome in the hominoids. Implications of the results for evolutionary studies and, in particular, for studies of the origin of modern humans, are discussed.     

Evolution of MHC class IIB in the genome of wild and ornamental guppies, Poecilia reticulata

This is the first study to quantify genomic sequence variation of the major histocompatibility complex (MHC) in wild and ornamental guppies, Poecilia reticulata. We sequenced 196-219 bp of exon 2 MHC class IIB (DAB) in 56 wild Trinidadian guppies and 14 ornamental strain guppies. Each of two natural populations possessed high allelic richness (15-16 alleles), whereas only three or fewer DAB alleles were amplified from ornamental guppies. The disparity in allelic richness between wild and ornamental fish cannot be fully explained by fixation of alleles by inbreeding, nor by the presence of non-amplified sequences (ie null alleles). Rather, we suggest that the same allele is fixed at duplicated MHC DAB loci owing to gene conversion. Alternatively, the number of loci in the ornamental strains has contracted during >100 generations in captivity, a hypothesis consistent with the accordion model of MHC evolution. We furthermore analysed the substitution patterns by making pairwise comparisons of sequence variation at the putative peptide binding region (PBR). The rate of non-synonymous substitutions (dN) only marginally exceeded synonymous substitutions (dS) in PBR codons. Highly diverged sequences showed no evidence for diversifying selection, possibly because synonymous substitutions have accumulated since their divergence. Also, the substitution pattern of similar alleles did not show evidence for diversifying selection, plausibly because advantageous non-synonymous substitutions have not yet accumulated. Intermediately diverged sequences showed the highest relative rate of non-synonymous substitutions, with dN/dS>14 in some pairwise comparisons. Consequently, a curvilinear relationship was observed between the dN/dS ratio and the level of sequence divergence.