Neutral Model Of Evolution Research Articles

Exact enumeration of cherries and pitchforks in ranked trees under the coalescent model

We consider exact enumerations and probabilistic properties of ranked trees when generated under the random coalescent process. Using a new approach, based on generating functions, we derive several statistics such as the exact probability of finding k cherries in a ranked tree of fixed size n. We then extend our method to consider also the number of pitchforks. We find a recursive formula to calculate the joint and conditional probabilities of cherries and pitchforks when the size of the tree is fixed. These results provide insights into structural properties of coalescent trees under the model of neutral evolution.

Adaptive morphological shifts to novel habitats in marine sculpin fishes

Sculpin fishes of the North American Pacific Coast provide an ideal opportunity to examine whether adaptive morphological character shifts have facilitated occupation of novel habitat types because of their well-described phylogeny and ecology. In this group, the basal-rooted species primarily occupy the subtidal habitat, whereas the species in the most distal clades are found in the intertidal. We tested multiple evolutionary models to determine whether changes in body size and changes in number of scales are adaptive for habitat use in sculpins. Based on a statistically robust, highly resolved molecular phylogeny of 26 species of sculpins, in combination with morphometric and habitat affinity data, our analyses show that an adaptive model based on habitat use best explains changes in body size and number of scales. The habitat model was statistically supported over models of neutral evolution, stabilizing selection across all habitats, and three clade-based models. We suggest that loss of scales and reduction of body size in the intertidal may facilitate cutaneous breathing in air when tidepools become hypoxic during low tides. This study demonstrates how the combined use of phylogenetic, ecological and statistical approaches helps to identify traits that are likely adaptive to novel habitats.

Nucleotide diversity and linkage disequilibrium of adaptive significant genes in Larix (Pinaceae)

Nucleotide polymorphism in genes potentially responsible for the adaptation to the latitudinal gradient of climate was investigated in the Siberian larch (Larix sibirica). The genes were selected among those previously studied in Douglas fir (Pseudotsuga menziesii), the variability of which indicates the influence of selection or shows the association with phenotypic traits that are important for adaptation to low temperatures and to geographically heterogeneous environment. Nucleotide sequences of eight genes, including coding and noncoding regions, were amplified and sequenced using newly developed primers specific to Larix. The sample size was 123 megagametophytes per locus in L. sibirica and 16.2 in L. occidentalis, which was taken as the out-group. The length of the sequenced fragments was 246-1700 bp. Nucleotide polymorphism π averaged 0.00536 (0.002 - 0.008), and haplotype diversity was H(d): 0.822 (0.625-0.948). Tajima's D was negative in all fragments and significant in three, while statistics D* and F* was significant in three and four segments, respectively, and F(s) was significant in three. This may indicate the presence of purifying selection on these genes or population growth. The HKA test revealed no significant deviations from the neutral model of evolution in all genes. The recombination parameter ρ/θ = 0.28 was close to the value obtained from P. menziesii. To investigate the association of polymorphic sites (factor) in these eight genes to the latitude of investigated individuals (trait), the generalized linear model (GLM) was used taking into account the population structure. After the False Discovery Rate correction no significant associations were found. The age of the split of American and Eurasian Larix lineages based on the nucleotide differences in the eight genes between L. sibirica and L. occidentalis is estimated to be 12 million years, which is much younger than the age of the most ancient Larix discoveries.

Peculiarities of Phosphoglycerate Kinase-1 Pseudogene Evolution in Schrenck Salamander (Salamandrella schrenckiiStrauch, 1870)

Processed copies of genes generally evolve in neutral mode as pseudogenes, however, some of them might be important sources of new functional genes. The psiPGK1 pseudogene has been discovered in Schrenck salamander (Salamandrella schrenckii, Amphibia, Caudata, Hynobiidae) via polymerase chain reaction used to amplify the phosphoglycerate kinase 1 gene (PGK1). This pseudogene is an intronless copy of PGK1 gene absent of exon 6. Analysis of psiPGK1 pseudogene polymorphism has demonstrated that it lacks mutations, which results in shifts in the stop codons and reading frames, as well as that the interspecies variation of this pseudogene was inconsistent with the neutral model of evolution. In addition, the pattern of phylogeographic differentiation of the psiPGK1 variants mainly coincides with that observed in mitochondrial DNA. These observations allow it to be suggested that the psiPGK1 pseudogene is a new functional gene in the Schrenck salamander.

Unravelling the Phellinus pini s.l. complex in North America: a multilocus phylogeny and differentiation analysis of Porodaedalea

SummaryPhellinus sensu lato (s.l.) is a complex of segregate genera that act as aggressive pathogens of woody plants. Nearly all of the genera in this complex have unresolved taxonomies, including Porodaedalea, which is one of the most important trunk rot pathogens of coniferous trees throughout the northern hemisphere. In an attempt to elucidate the species within Porodaedalea, a multilocus phylogenetic analysis was performed with partial sequences from four loci (internal transcribed spacer, nuclear large subunit, tef1 and rpb2) using 41 isolates that originated from North America and Europe. For reference, we analysed the neotype isolates of Porodaedalea pini and P. chrysoloma. Our results confirmed that Porodaedalea pini s.s. and P. chrysoloma s.s. are unique phylogenetic species that do not occur in North America. We detected two discrete clades of Porodaedalea originating from the southwestern and southeastern United States. Isolates from these regions grouped with significant statistical support and represent undescribed taxa. With the exception of P. cancriformans, our analyses revealed monophyly among 28 isolates originating from the northern United States, Canada and Fennoscandia, a group we have labelled the ‘Holarctic group’. Holarctic group isolates were collected from Larix, Picea, Pinus, Pseudotsuga and Tsuga and were presumed to represent at least four morphological species (P. gilbertsonii, P. laricis, P. pini s.l. and P. piceina). Tests of gene flow and genetic differentiation detected significant differences among Holarctic group isolates by region of origin, and three subgroups were designated: (i) Atlantic‐Boreal; (ii) Interior; and (iii) Pacific. Neutrality tests using the Holarctic group demonstrated significant departures from the standard neutral model of evolution and could indicate that a diversifying selection has maintained rare phenotypes in the population, which has fostered taxonomic confusion in Porodaedalea.

Nucleotide sequence diversity of floral pigment genes in Mexican populations of Ipomoea purpurea (morning glory) accord with a neutral model of evolution

The common morning glory (Ipomoea purpurea) is an annual vine native to Central and Southern Mexico. The genetics of flower color polymorphisms and interactions with the biotic environment have been extensively studied in I. purpurea and in its sister species I. nil. In this study, we examine nucleotide sequence polymorphism in 11 loci, 9 of which are known to participate in a pathway that produces floral pigments. A sample of 30 I. purpurea accessions from the native range of Central and Southern Mexico comprise the data, along with one accession from each of the two sister species I. alba and I. nil. We observe moderate levels of nucleotide sequence polymorphism of ~1%. The ratio of recombination to mutation parameter estimates (ρ/θ) of ~2.5 appears consistent with a mixed-mating system. Ipomoea resequencing data from these genic regions are noteworthy in providing a good fit to the standard neutral model of molecular evolution. The derived silent site frequency spectrum is very close to that predicted by coalescent simulations of a drift-mutation process, and Tajima's D values are not significantly different from expectations under neutrality.

Nucleotide Sequence Diversity of Floral Pigment Genes in Mexican Populations of Ipomoea purpurea (Morning Glory) Accord with a Neutral Model of Evolution

The common morning glory (Ipomoea purpurea) is an annual vine native to Central and Southern Mexico. The genetics of flower color polymorphisms and interactions with the biotic environment have been extensively studied in I. purpurea and in its sister species I. nil. In this study, we examine nucleotide sequence polymorphism in 11 loci, 9 of which are known to participate in a pathway that produces floral pigments. A sample of 30 I. purpurea accessions from the native range of Central and Southern Mexico comprise the data, along with one accession from each of the two sister species I. alba and I. nil. We observe moderate levels of nucleotide sequence polymorphism of ~1%. The ratio of recombination to mutation parameter estimates (ρ/θ) of ~2.5 appears consistent with a mixed-mating system. Ipomoea resequencing data from these genic regions are noteworthy in providing a good fit to the standard neutral model of molecular evolution. The derived silent site frequency spectrum is very close to that predicted by coalescent simulations of a drift-mutation process, and Tajima's D values are not significantly different from expectations under neutrality.

Reduced Selective Constraint in Endosymbionts: Elevation in Radical Amino Acid Replacements Occurs Genome-Wide

As predicted by the nearly neutral model of evolution, numerous studies have shown that reduced Ne accelerates the accumulation of slightly deleterious changes under genetic drift. While such studies have mostly focused on eukaryotes, bacteria also offer excellent models to explore the effects of Ne. Most notably, the genomes of host-dependent bacteria with small Ne show signatures of genetic drift, including elevated Ka/Ks. Here, I explore the utility of an alternative measure of selective constraint: the per-site rate of radical and conservative amino acid substitutions (Dr/Dc). I test the hypothesis that purifying selection against radical amino acid changes is less effective in two insect endosymbiont groups (Blochmannia of ants and Buchnera of aphids), compared to related gamma-Proteobacteria. Genome comparisons demonstrate a significant elevation in Dr/Dc in endosymbionts that affects the majority (66–79%) of shared orthologs examined. The elevation of Dr/Dc in endosymbionts affects all functional categories examined. Simulations indicate that Dr/Dc estimates are sensitive to codon frequencies and mutational parameters; however, estimation biases occur in the opposite direction as the patterns observed in genome comparisons, thereby making the inference of elevated Dr/Dc more conservative. Increased Dr/Dc and other signatures of genome degradation in endosymbionts are consistent with strong effects of genetic drift in their small populations, as well as linkage to selected sites in these asexual bacteria. While relaxed selection against radical substitutions may contribute, genome-wide processes such as genetic drift and linkage best explain the pervasive elevation in Dr/Dc across diverse functional categories that include basic cellular processes. Although the current study focuses on a few bacterial lineages, it suggests Dr/Dc is a useful gauge of selective constraint and may provide a valuable alternative to Ka/Ks when high sequence divergences preclude estimates of Ks. Broader application of Dr/Dc will benefit from approaches less prone to estimation biases.

Molecular Population Genetics of Elicitor-Induced Resistance Genes in European Aspen (Populus tremula L., Salicaceae)

Owing to their long life span and ecological dominance in many communities, forest trees are subject to attack from a diverse array of herbivores throughout their range, and have therefore developed a large number of both constitutive and inducible defenses. We used molecular population genetics methods to examine the evolution of eight genes in European aspen, Populus tremula, that are all associated with defensive responses against pests and/or pathogens, and have earlier been shown to become strongly up-regulated in poplars as a response to wounding and insect herbivory. Our results show that the majority of these defense genes show patterns of intraspecific polymorphism and site-frequency spectra that are consistent with a neutral model of evolution. However, two of the genes, both belonging to a small gene family of polyphenol oxidases, show multiple deviations from the neutral model. The gene PPO1 has a 600 bp region with a highly elevated KA/KS ratio and reduced synonymous diversity. PPO1 also shows a skew toward intermediate frequency variants in the SFS, and a pronounced fixation of non-synonymous mutations, all pointing to the fact that PPO1 has been subjected to recurrent selective sweeps. The gene PPO2 shows a marked excess of high frequency, derived variants and shows many of the same trends as PPO1 does, even though the pattern is less pronounced, suggesting that PPO2 might have been the target of a recent selective sweep. Our results supports data from both Populus and other species which have found that the the majority of defense-associated genes show few signs of selection but that a number of genes involved in mediating defense against herbivores show signs of adaptive evolution.

The characteristics of allelic polymorphism in killer-immunoglobulin-like receptor framework genes in African Americans

The frequencies of alleles of killer cell immunoglobulin-like receptor genes, KIR3DL3 and KIR3DL2, and the carrier frequency of KIR2DL4 alleles have been determined from a population of African Americans (n = 100) by DNA sequencing of the coding regions. Fifty alleles of KIR3DL3 were observed with the most frequent, KIR3DL3*00901 (13%). KIR3DL2 was also diverse; 32 alleles with KIR3DL2*00103 the most frequent (17%). For KIR2DL4, of the 18 alleles observed, one allele, KIR2DL4*00103, was found in 64 of the 100 individuals. Thirty-six novel alleles encoding a total of 28 unique receptors are described. Pairwise comparisons among all of the alleles at each locus suggest a predominance of synonymous substitutions. The variation at all three framework loci fits a neutral model of evolution.

Evolutionary patterns and processes in the radiation of phyllostomid bats

BackgroundThe phyllostomid bats present the most extensive ecological and phenotypic radiation known among mammal families. This group is an important model system for studies of cranial ecomorphology and functional optimisation because of the constraints imposed by the requirements of flight. A number of studies supporting phyllostomid adaptation have focused on qualitative descriptions or correlating functional variables and diet, but explicit tests of possible evolutionary mechanisms and scenarios for phenotypic diversification have not been performed. We used a combination of morphometric and comparative methods to test hypotheses regarding the evolutionary processes behind the diversification of phenotype (mandible shape and size) and diet during the phyllostomid radiation.ResultsThe different phyllostomid lineages radiate in mandible shape space, with each feeding specialisation evolving towards different axes. Size and shape evolve quite independently, as the main directions of shape variation are associated with mandible elongation (nectarivores) or the relative size of tooth rows and mandibular processes (sanguivores and frugivores), which are not associated with size changes in the mandible. The early period of phyllostomid diversification is marked by a burst of shape, size, and diet disparity (before 20 Mya), larger than expected by neutral evolution models, settling later to a period of relative phenotypic and ecological stasis. The best fitting evolutionary model for both mandible shape and size divergence was an Ornstein-Uhlenbeck process with five adaptive peaks (insectivory, carnivory, sanguivory, nectarivory and frugivory).ConclusionsThe radiation of phyllostomid bats presented adaptive and non-adaptive components nested together through the time frame of the family's evolution. The first 10 My of the radiation were marked by strong phenotypic and ecological divergence among ancestors of modern lineages, whereas the remaining 20 My were marked by stasis around a number of probable adaptive peaks. A considerable amount of cladogenesis and speciation in this period is likely to be the result of non-adaptive allopatric divergence or adaptations to peaks within major dietary categories.

Low genetic polymorphism of merozoite surface proteins 7 and 10 in Colombian Plasmodium vivax isolates

The merozoite surface protein (MSP) family is involved in the initial interaction between merozoites and erythrocytes in Plasmodium species, its members are therefore becoming major targets for vaccine development. Considering that antigens included in a subunit malaria vaccine should be both accessible to the immune system and lack genetic diversity or have very limited polymorphism, we have analyzed the genetic diversity of three msp genes (msp-7A, msp-7K and msp-10) in different geographical regions of Colombia. The results showed that these genes follow the neutral model of evolution and also display low genetic diversity. The strong conservation found for msp-7 haplotypes in isolates from geographically different regions further suggests that these proteins could be good components of a vaccine against Plasmodium vivax malaria, thereby avoiding strain-specific immune responses.

Extended Phylogeny of the Equine Arteritis Virus Sequences Including South American Sequences

Objective: To perform genetic analysis of the ORF5 of equine arteritis virus (EAV) may provide new insights into the genetic evolution and origin of the Argentinean EAV sequences. Methods: A total of 76 sequences were analyzed by neighbor joining (NJ), maximum parsimony and maximum likelihood algorithms. The analysis of the selective pressures was performed using the Tajima’s test. Results: The trees showed similar topologies. Two clades were identified: the first clade was formed by strains isolated mainly from a donkey, whereas the second clade presented four large groups from different geographic regions exclusively from Equus caballus. In this clade, we identified a group formed by South African and another one by South American and European sequences. In the latter, the monophyletic group was formed by seven Argentinean sequences. In the NJ tree, we identified a group formed by six Argentinean sequences. The Tajima’s test showed a D value of 1.73663, indicating that the sequences analyzed follow a neutral evolution model. Conclusion: We concluded that the Argentinean sequences have a paraphyletic origin and that the fixation of point mutation might follow the neutral model evolution; however, we identified purifying pressures that may be involved in the differentiation of the EAV sequences.

Effects of selection and mutation on mitochondrial variation and inferences of historical population expansion in a Caribbean reef fish

Mitochondrial sequence data has contributed to the understanding of historical demography through the application of neutrality tests and coalescent estimators of population growth rates. Characteristics of the mitochondrial genome, such as high mutation rate and lack of recombination, render it particularly suitable for these types of studies. However, selection can also affect patterns of mitochondrial variation. Furthermore, conclusions based on single mitochondrial loci can be sensitive to differences among mitochondrial genes in mutation rate and pattern and levels of homoplasy. We tested the contributions of these factors to patterns of mitochondrial variation in the Caribbean reef fish Halichoeres bivittatus using a multilocus sequence-based approach. Mitochondrial protein-coding loci deviated strongly from a neutral model of evolution and indicated high rates of estimated population growth. In contrast, the mitochondrial control region and a nuclear intron showed little or no deviation from neutrality and low estimated growth rates. The level of variation among loci is inconsistent with a demographic explanation and likely stems from the influence of mutation and selection on the mitochondrial genome. In H. bivittatus, a finding of high rates of population growth is likely an artifact of selection on mitochondrial proteins. This result suggests caution in the interpretation of variation at single mitochondrial loci, and highlights the importance of the use of unlinked nuclear loci to test demographic inferences made from mitochondrial DNA.

Computing and visually analyzing mutual information in molecular co-evolution.

BackgroundSelective pressure in molecular evolution leads to uneven distributions of amino acids and nucleotides. In fact one observes correlations among such constituents due to a large number of biophysical mechanisms (folding properties, electrostatics, ...). To quantify these correlations the mutual information -after proper normalization - has proven most effective. The challenge is to navigate the large amount of data, which in a study for a typical protein cannot simply be plotted.ResultsTo visually analyze mutual information we developed a matrix visualization tool that allows different views on the mutual information matrix: filtering, sorting, and weighting are among them. The user can interactively navigate a huge matrix in real-time and search e.g., for patterns and unusual high or low values. A computation of the mutual information matrix for a sequence alignment in FASTA-format is possible. The respective stand-alone program computes in addition proper normalizations for a null model of neutral evolution and maps the mutual information to Z-scores with respect to the null model.ConclusionsThe new tool allows to compute and visually analyze sequence data for possible co-evolutionary signals. The tool has already been successfully employed in evolutionary studies on HIV1 protease and acetylcholinesterase. The functionality of the tool was defined by users using the tool in real-world research. The software can also be used for visual analysis of other matrix-like data, such as information obtained by DNA microarray experiments. The package is platform-independently implemented in Java and free for academic use under a GPL license.

Adaptive evolution in subterranean termite antifungal peptides

We identified and analysed mRNA sequences of two immune proteins from the subterranean termites Reticulitermes flavipes and Reticulitermes virginicus. These proteins correspond to two immune proteins described in the distantly related termite genus Nasutitermes; termicin, which is a small antifungal peptide, and GNBP2, which functions both as a broad pattern recognition receptor and a direct antifungal effector. A population genetic analysis of nucleotide intraspecific polymorphism and interspecific divergence indicates that a selective sweep has reduced polymorphism in the termicins. Moreover, this selective sweep appears to have been driven by the positive selection of beneficial amino acid changes in the antifungal peptide. In contrast, the pattern of polymorphism and divergence in GNBP2 corresponds to the standard neutral model of evolution.

Mathematical Modeling of Ultradeep Sequencing Data Reveals that Acute CD8 + T-Lymphocyte Responses Exert Strong Selective Pressure in Simian Immunodeficiency Virus-Infected Macaques but Still Fail To Clear Founder Epitope Sequences

The prominent role of antiviral cytotoxic CD8(+) T-lymphocytes (CD8-TL) in containing the acute viremia of human and simian immunodeficiency viruses (HIV-1 and SIV) has rationalized the development of T-cell-based vaccines. However, the presence of escape mutations in the acute stage of infection has raised a concern that accelerated escape from vaccine-induced CD8-TL responses might undermine vaccine efficacy. We reanalyzed previously published data of 101,822 viral genomes of three CD8-TL epitopes, Nef(103-111)RM9 (RM9), Tat(28-35)SL8 (SL8), and Gag(181-189)CM9 (CM9), sampled by ultradeep pyrosequencing from eight macaques. Multiple epitope variants appeared during the resolution of acute viremia, followed by the predominance of a single mutant epitope. By fitting a mathematical model, we estimated the first acute escape rate as 0.36 day(-1) within escape-prone epitopes, RM9 and SL8, and the chronic escape rate as 0.014 day(-1) within the CM9 epitope. Our estimate of SIV acute escape rates was found to be comparable to very early HIV-1 escape rates. The timing of the first escape was more highly correlated with the timing of the peak CD8-TL response than with the magnitude of the CD8-TL response. The transmitted epitope decayed more than 400 times faster during the acute viral decline stage than predicted by a neutral evolution model. However, the founder epitope persisted as a minor population even at the viral set point; in contrast, the majority of acute escape epitopes were completely cleared. Our results suggest that a reservoir of SIV infection is preferentially formed by virus with the transmitted epitope.

Polymorphisms in the CPB2 Gene Are Maintained by Balancing Selection and Result in Haplotype-Preferential Splicing of Exon 7

The CPB2 gene encodes thrombin-activatable fibrinolysis inhibitor (TAFI), a hepatically secreted zymogen acting as a molecular link among coagulation, fibrinolysis, and inflammation. Variants in CPB2 have been associated with several human conditions. We resequenced and analyzed the two regions carrying previously known nonsynonimous single-nucleotide polymorphisms (Ala147Thr and Ile325Thr) and variants affecting transcript stability. Our data indicate that whereas the gene portion extending from exon 9 to the 3' untranslated region fits a model of neutral evolution, variants in the region encompassing exons 6-7 have been maintained by balancing selection. Indeed, we verified that the region displays high nucleotide diversity, many intermediate frequency variants, and an excess of polymorphism compared with interspecific divergence. Consistently, haplotype analysis indicated the presence of two major haplotype clades separated by deep branches. Transcript analysis revealed that in both HepG2 cells and human liver samples, CPB2 exon 7 undergoes haplotype-preferential skipping. Therefore, we indicate that balancing selection has been maintaining functional variants that promote alternative exon 7 splicing. Although transcripts lacking exon 7 represent a minority of total CPB2 products, the effect on antifibrinolytic activity might be much greater as the intrinsic instability of TAFI is a major determinant of its antifibrinolytic potential. These data highlight the contribution of population genetics approaches to the analysis of functional genetic variation and may orient further biochemical and genetics studies on the pathophysiologic role of CPB2 gene products.

Selection at the MHC class IIB locus across guppy (Poecilia reticulata) populations

The highly diverse genes of the major histocompatibility complex (MHC) are important in the adaptive immune system and are expected to be under selection from pathogens. Thus, the MHC genes provide an exceptional opportunity to investigate patterns of selection within and across populations. In this study, we analyzed genetic variation at the MHC class IIB gene and six microsatellite loci across 10 populations of guppies (Poecilia reticulata) in the northern range of Trinidad. We found a high level of diversity at the MHC, with a total of 43 alleles in 142 individuals. At the population level, we found that neutral evolution could not fully account for the variability found at the MHC. Instead, we found that MHC F(ST) statistics were lower than F(ST) derived from the microsatellite loci; 33 of 45 population pairwise estimates for the MHC were significantly lower than those for the microsatellite loci, and MHC F(ST) estimates were consistently lower than those predicted by a coalescent model of neutral evolution. These results suggest a similar selection acting across populations, and we discuss the potential roles of directional and balancing selection. At the sequence level, we found evidence for both positive and purifying selection. Furthermore, positive selection was detected within and adjacent to the putative peptide-binding region (PBR) of the MHC. Surprisingly, we also found a purifying selection at two sites within the putative PBR. Overall, our data provide evidence for selection for functional diversity at the MHC class IIB gene at both the population and nucleotide levels of guppy populations.

Modeling language change: An evaluation of Trudgill's theory of the emergence of New Zealand English

AbstractTrudgill (2004) proposed that the emergence of New Zealand English, and of isolated new dialects generally, is purely deterministic. It can be explained solely in terms of the frequency of occurrence of particular variants and the frequency of interactions between different speakers in the society. Trudgill's theory is closely related to usage-based models of language, in which frequency plays a role in the representation of linguistic knowledge and in language change. Trudgill's theory also corresponds to a neutral evolution model of language change. We use a mathematical model based on Croft's usage-based evolutionary framework for language change (Baxter, Blythe, Croft, & McKane, 2006), and investigate whether Trudgill's theory is a plausible model of the emergence of new dialects. The results of our modeling indicate that determinism cannot be a sufficient mechanism for the emergence of a new dialect. Our approach illustrates the utility of mathematical modeling of theories and of empirical data for the study of language change.