Neutral Model Of Molecular Evolution Research Articles

Interpreting the pervasive observation of U-shaped Site Frequency Spectra.

The standard neutral model of molecular evolution has traditionally been used as the null model for population genomics. We gathered a collection of 45 genome-wide site frequency spectra from a diverse set of species, most of which display an excess of low and high frequency variants compared to the expectation of the standard neutral model, resulting in U-shaped spectra. We show that multiple merger coalescent models often provide a better fit to these observations than the standard Kingman coalescent. Hence, in many circumstances these under-utilized models may serve as the more appropriate reference for genomic analyses. We further discuss the underlying evolutionary processes that may result in the widespread U-shape of frequency spectra.

Molecular Phylogenetics and Population Genetics of the Dengue Vector Aedes aegypti From the Arabian Peninsula.

Aedes aegypti mosquito is the principal dengue vector in the Kingdom of Saudi Arabia (KSA); however, no study has addressed its ecology and population structure yet. Therefore, we report on Ae. aegypti phylo- and population genetics using three DNA markers: COI, ND4, and rDNA-ITS2. Sampling the immature stages of Ae. aegypti revealed that water storage tanks (34.3% of habitats) were the most productive and contained 33% of immatures stages. Other important habitats included containers for wastewater drainage (including air-conditioning and water cooler trays) and containers associated with ornamentation. Shallow water leakage spots (2.7% of habitats, 8% of immatures) can be considered rare-but-epidemiologically-important containers. Neighbor-joining (NJ) phylogenetic analysis of Ae. aegypti identified 8, 14, and 9 haplotypes of COI, ND4, and ITS2, respectively, and revealed high levels of genetic variation in Ae. aegypti populations of KSA. Global distribution of haplotypes also indicated multiple gene introductions into these populations, with high levels of intra-population genetic variation and continuous gene exchange. The neutrality values indicated a deficiency of alleles and suggested that the KSA Ae. aegypti loci tested did not follow a neutral model of molecular evolution. Fst values and AMOVA indicated that most of the genetic variation in the KSA Ae. aegypti populations is due to intra- rather than inter-population differences. This is the first comprehensive report on the phylo- and population genetics of Ae. aegypti from the Arabian Peninsula. This information expands our understanding of the ecology and population dynamics of this important arboviral vector for informed control efforts.

How to Make a Rodent Giant: Genomic Basis and Tradeoffs of Gigantism in the Capybara, the World's Largest Rodent.

Gigantism results when one lineage within a clade evolves extremely large body size relative to its small-bodied ancestors, a common phenomenon in animals. Theory predicts that the evolution of giants should be constrained by two tradeoffs. First, because body size is negatively correlated with population size, purifying selection is expected to be less efficient in species of large body size, leading to increased mutational load. Second, gigantism is achieved through generating a higher number of cells along with higher rates of cell proliferation, thus increasing the likelihood of cancer. To explore the genetic basis of gigantism in rodents and uncover genomic signatures of gigantism-related tradeoffs, we assembled a draft genome of the capybara (Hydrochoerus hydrochaeris), the world’s largest living rodent. We found that the genome-wide ratio of nonsynonymous to synonymous mutations (ω) is elevated in the capybara relative to other rodents, likely caused by a generation-time effect and consistent with a nearly neutral model of molecular evolution. A genome-wide scan for adaptive protein evolution in the capybara highlighted several genes controlling postnatal bone growth regulation and musculoskeletal development, which are relevant to anatomical and developmental modifications for an increase in overall body size. Capybara-specific gene-family expansions included a putative novel anticancer adaptation that involves T-cell-mediated tumor suppression, offering a potential resolution to the increased cancer risk in this lineage. Our comparative genomic results uncovered the signature of an intragenomic conflict where the evolution of gigantism in the capybara involved selection on genes and pathways that are directly linked to cancer.

Molecular evolution of the VacA p55 binding domain of Helicobacter pylori in mestizos from a high gastric cancer region of Colombia.

The stomach bacterium Helicobacter pylori is one of the most prevalent pathogens in humans, closely linked with serious diseases such as gastric cancer. The microbe has been associated with its host for more than 100,000 years and escorted modern humans out of Africa. H. pylori is predominantly transmitted within families and dispersed globally, resulting in distinct phylogeographic patterns, which can be utilized to investigate migrations and bioturbation events in human history. Latin America was affected by several human migratory waves due to the Spanish colonisation that drastically changed the genetic load and composition of the bacteria and its host. Genetic evidence indicates that independent evolutionary lines of H. pylori have evolved in mestizos from Colombia and other countries in the region during more than 500 years since colonisation. The vacuolating cytotoxin VacA represents a major virulence factor of the pathogen comprising two domains, p33 and p55, the latter of which is essential for binding to the host epithelial cell. The evolution of the VacA toxin in Colombia has been strongly biased due to the effects of Spanish colonization. However, the variation patterns and microevolution of the p55 domain have not yet been described for this population. In the present study, we determined the genetic polymorphisms and deviations in the neutral model of molecular evolution in the p55 domain of 101 clinical H. pylori isolates collected in Bogotá, a city located in Andean mountains characterized by its high gastric cancer risk and its dominant mestizo population. The microevolutionary patterns of the p55 domain were shaped by recombination, purifying and episodic diversifying positive selection. Furthermore, amino acid positions 261 and 321 in the p55 domain of VacA show a high variability among mestizos clinical subsets, suggesting that natural selection in H. pylori may operate differentially in patients with different gastric diseases.

Rapid evolution of the Helicobacter pylori AlpA adhesin in a high gastric cancer risk region from Colombia.

To be able to survive, Helicobacter pylori must adhere to the gastric epithelial cells of its human host. For this purpose, the bacterium employs an array of adhesins, for example, AlpA. The adhesin AlpA has been proposed as a major adhesin because of its critical role in human stomach colonization. Therefore, understanding how AlpA evolved could be important for the development of new diagnostic strategies. However, the genetic variation and microevolutionary patterns of alpA have not been described in Colombia. The study aim was to describe the variation patterns and microevolutionary process of alpA in Colombian clinical isolates of H. pylori. The existing polymorphisms, which are deviations from the neutral model of molecular evolution, and the genetic differentiation of the alpA gene from Colombian clinical isolates of H. pylori were determined. The analysis shows that gene conversion and purifying selection have shaped the evolution of three different variants of alpA in Colombia.

Fluctuating Selection in the Moran.

Contrary to classical population genetics theory, experiments demonstrate that fluctuating selection can protect a haploid polymorphism in the absence of frequency dependent effects on fitness. Using forward simulations with the Moran model, we confirm our analytical results showing that a fluctuating selection regime, with a mean selection coefficient of zero, promotes polymorphism. We find that increases in heterozygosity over neutral expectations are especially pronounced when fluctuations are rapid, mutation is weak, the population size is large, and the variance in selection is big. Lowering the frequency of fluctuations makes selection more directional, and so heterozygosity declines. We also show that fluctuating selection raises dn/ds ratios for polymorphism, not only by sweeping selected alleles into the population, but also by purging the neutral variants of selected alleles as they undergo repeated bottlenecks. Our analysis shows that randomly fluctuating selection increases the rate of evolution by increasing the probability of fixation. The impact is especially noticeable when the selection is strong and mutation is weak. Simulations show the increase in the rate of evolution declines as the rate of new mutations entering the population increases, an effect attributable to clonal interference. Intriguingly, fluctuating selection increases the dn/ds ratios for divergence more than for polymorphism, a pattern commonly seen in comparative genomics. Our model, which extends the classical neutral model of molecular evolution by incorporating random fluctuations in selection, accommodates a wide variety of observations, both neutral and selected, with economy.

Disparate Modes of Evolution Shaped Modern Prion (PRNP) and Prion-Related Doppel (PRND) Variation in Domestic Cattle.

Previous investigations aimed at determining whether the mammalian prion protein actually facilitates tangible molecular aspects of either a discrete or pleiotropic functional niche have been debated, especially given the apparent absence of overt behavioral or physiological phenotypes associated with several mammalian prion gene (PRNP) knockout experiments. Moreover, a previous evaluation of PRNP knockout cattle concluded that they were normal, suggesting that the bovine prion protein is physiologically dispensable. Herein, we examined the frequency and distribution of nucleotide sequence variation within the coding regions of bovine PRNP and the adjacent Doppel (PRND) gene, a proximal paralogue to PRNP on BTA13. Evaluation of PRND variation demonstrated that the gene does not depart from a strictly neutral model of molecular evolution, and would therefore not be expected to influence tests of selection within PRNP. Collectively, our analyses confirm that intense purifying selection is indeed occurring directly on bovine PRNP, which is indicative of a protein with an important role. These results suggest that the lack of observed fitness effects may not manifest in the controlled environmental conditions used to care for and raise PRNP knockout animals.

Genetic diversity and molecular evolution of the internal transcribed spacer (ITSs) of nuclear ribosomal DNA in the Tunisian fig cultivars (Ficus carica L.; Moracea)

Abstract Nuclear ribosomal DNAs were explored to establish genetic relationships among Ficus carica cultivars and elucidate its molecular evolution. Results suggest the occurrence of haplotypic and nucleotide diversity. The neighbour-joining dendrograms show a continuous diversity that characterize local resources. Furthermore, our results demonstrated that the ITS2 spacer is seating to a larger number of substitutions than the ITS1 spacer. Sequence analysis demonstrates that the ITS2 spacer is evolving 1.12 times faster than the ITS1 one. The ratio of transition/transversion of 0.278 suggests that the 5.8S gene is evolving 2.84 and 3.20 times less rapidly than the spacers ITS1 and ITS2, respectively. Molecular evolution analysis confirmed an explicit rejection of the null hypothesis in F. carica . ITS1, ITS2 spacers and the 5.8S gene evolved under a strictly neutral model of molecular evolution. A scenario of positive selection and recent expansion of F. carica genotypes across Tunisia seems to be retained.

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.

Experimental illumination of a fitness landscape

The genes of all organisms have been shaped by selective pressures. The relationship between gene sequence and fitness has tremendous implications for understanding both evolutionary processes and functional constraints on the encoded proteins. Here, we have exploited deep sequencing technology to experimentally determine the fitness of all possible individual point mutants under controlled conditions for a nine-amino acid region of Hsp90. Over the past five decades, limited glimpses into the relationship between gene sequence and function have sparked a long debate regarding the distribution, relative proportion, and evolutionary significance of deleterious, neutral, and advantageous mutations. Our systematic experimental measurement of fitness effects of Hsp90 mutants in yeast, evaluated in the light of existing population genetic theory, are remarkably consistent with a nearly neutral model of molecular evolution.

Sequence analysis of coding DNA fragments of pfcrt and pfmdr-1 genes in Plasmodium falciparum isolates from Odisha, India

The global emergence and spread of malaria parasites resistant to antimalarial drugs is the major problem in malaria control. The genetic basis of the parasite's resistance to the antimalarial drug chloroquine (CQ) is well-documented, allowing for the analysis of field isolates of malaria parasites to address evolutionary questions concerning the origin and spread of CQ-resistance. Here, we present DNA sequence analyses of both the second exon of the Plasmodium falciparum CQ-resistance transporter (pfcrt) gene and the 5' end of the P. falciparum multidrug-resistance 1 (pfmdr-1) gene in 40 P. falciparum field isolates collected from eight different localities of Odisha, India. First, we genotyped the samples for the pfcrt K76T and pfmdr-1 N86Y mutations in these two genes, which are the mutations primarily implicated in CQ-resistance. We further analyzed amino acid changes in codons 72-76 of the pfcrt haplotypes. Interestingly, both the K76T and N86Y mutations were found to co-exist in 32 out of the total 40 isolates, which were of either the CVIET or SVMNT haplotype, while the remaining eight isolates were of the CVMNK haplotype. In total, eight nonsynonymous single nucleotide polymorphisms (SNPs) were observed, six in the pfcrt gene and two in the pfmdr-1 gene. One poorly studied SNP in the pfcrt gene (A97T) was found at a high frequency in many P. falciparum samples. Using population genetics to analyze these two gene fragments, we revealed comparatively higher nucleotide diversity in the pfcrt gene than in the pfmdr-1 gene. Furthermore, linkage disequilibrium was found to be tight between closely spaced SNPs of the pfcrt gene. Finally, both the pfcrt and the pfmdr-1 genes were found to evolve under the standard neutral model of molecular evolution.

Multilocus nuclear DNA markers and genetic parameters in an Indian Anopheles minimus population

Estimation of population genetic parameters is highly dependent on the choice of genetic markers. Furthermore, inferences based on single genes could lead to erroneous conclusions and population genetic outcomes, thus usage of multiple loci is suggested. Considering malaria is a highly fatal vector-borne infectious disease, inference on population genetic structure and demography could be of help in the long run for malaria vector management and control. Using the published genome sequence information of Anopheles gambiae we designed EPIC primers to amplify DNA fragments in An. minimus nuclear genome. Eight such DNA fragments could be successfully amplified and sequenced and homology to corresponding genes of An. gambiae was established. All the eight DNA fragments were found to be polymorphic for single nucleotide polymorphisms (SNPs) in a population sample of An. minimus from India. Several tests of neutrality confirmed that all the eight fragments evolve under a standard neutral model of molecular evolution. Furthermore, multilocus linkage disequilibrium studies revealed that the DNA fragments were not genetically linked to each other and thus are independently evolving. Tests of past population demographic events clearly revealed that this Indian population of An. minimus follows demographic equilibrium model, without any significant recent population bottleneck or expansion. The eight multilocus nuclear DNA fragments thus could be considered as ‘putatively neutral’ and be used to infer population structure and demographic history of An. minimus, a major malaria vector in the Southeast Asia and India. Moreover, the estimations of population demography using these putatively neutral markers can provide a baseline against which, test for the role of natural selection in functionally relevant genes of An. minimus would be possible.

Molecular evolution of chloroplast DNA in fig ( Ficus carica L.): Footprints of sweep selection and recent expansion

Here, we report the nucleotide variation in two non-coding regions of the chloroplast DNA (cpDNA) to construct a possible evolutionary scenario in Ficus carica L. Our results suggest the occurrence of haplotypic and nucleotide diversity with a large variation level of chloroplast non-coding regions. Furthermore, our results demonstrated an explicit rejection of the null hypothesis that within F. carica the intron trnL and the spacer trnL- trnF evolved under a strictly neutral model of molecular evolution. Although, recent population expansion could serve as one alternative explanation for the detected excess of singleton, our results imply a positive selection and the genetic hitchhiking effect is unlikely. Parameters performed supported scenario of sweep selection and recent expansion of F. carica across Tunisia. Our results indicate that both positive selection and demographic histories have jointly contributed to the observed patterns of nucleotide diversity and haplotypes structure. Based on the results, we characterize the fig resources and provide several suggestions for effective conservation and improvement programs.

Natural selection on cis and trans regulation in yeasts

Gene expression is regulated both by cis elements, which are DNA segments closely linked to the genes they regulate, and by trans factors, which are usually proteins capable of diffusing to unlinked genes. Understanding the patterns and sources of regulatory variation is crucial for understanding phenotypic and genome evolution. Here, we measure genome-wide allele-specific expression by deep sequencing to investigate the patterns of cis and trans expression variation between two strains of Saccharomyces cerevisiae. We propose a statistical modeling framework based on the binomial distribution that simultaneously addresses normalization of read counts derived from different parents and estimating the cis and trans expression variation parameters. We find that expression polymorphism in yeast is common for both cis and trans, though trans variation is more common. Constraint in expression evolution is correlated with other hallmarks of constraint, including gene essentiality, number of protein interaction partners, and constraint in amino acid substitution, indicating that both cis and trans polymorphism are clearly under purifying selection, though trans variation appears to be more sensitive to selective constraint. Comparing interspecific expression divergence between S. cerevisiae and S. paradoxus to our intraspecific variation suggests a significant departure from a neutral model of molecular evolution. A further examination of correlation between polymorphism and divergence within each category suggests that cis divergence is more frequently mediated by positive Darwinian selection than is trans divergence.

Rapid evolution of yeast centromeres in the absence of drive.

To find the most rapidly evolving regions in the yeast genome we compared most of chromosome III from three closely related lineages of the wild yeast Saccharomyces paradoxus. Unexpectedly, the centromere appears to be the fastest-evolving part of the chromosome, evolving even faster than DNA sequences unlikely to be under selective constraint (i.e., synonymous sites after correcting for codon usage bias and remnant transposable elements). Centromeres on other chromosomes also show an elevated rate of nucleotide substitution. Rapid centromere evolution has also been reported for some plants and animals and has been attributed to selection for inclusion in the egg or the ovule at female meiosis. But Saccharomyces yeasts have symmetrical meioses with all four products surviving, thus providing no opportunity for meiotic drive. In addition, yeast centromeres show the high levels of polymorphism expected under a neutral model of molecular evolution. We suggest that yeast centromeres suffer an elevated rate of mutation relative to other chromosomal regions and they change through a process of "centromere drift," not drive.

Mammalian Nonsynonymous Sites Are Not Overdispersed: Comparative Genomic Analysis of Index of Dispersion of Mammalian Proteins

It is often stated that patterns of nonsynonymous rate variation among mammalian lineages are more irregular than expected or overdispersed under the neutral model, whereas synonymous sites conform to the neutral model. Here we reexamined genome-wide patterns of the variance to mean ratio, or index of dispersion (R), of substitutions in proteins from human, mouse, and dog. Contrary to the prevailing notion, we found that the mean index of dispersion for nonsynonymous sites of mammalian proteins is not significantly different from 1. We propose that earlier analyses were biased because the data included disproportionately more protein hormones, which tend to be more dispersed than genes in other functional categories. Synonymous sites exhibit greater degree of dispersion than nonsynonymous sites, although similar to earlier estimates and potentially due to errors associated with correction for multiple hits. Overall, our analysis identifies strong genome-wide generation-time effect and natural selection as important determinants of among-lineage variation of protein evolutionary rates. Furthermore, patterns of lineage-specific selective constraint are consistent with the nearly neutral model of molecular evolution.

DNA Sequence Polymorphism and Divergence at the erect wing and suppressor of sable Loci of Drosophila melanogaster and D. simulans

Several evolutionary models of linked selection (e.g., genetic hitchhiking, background selection, and random environment) predict a reduction in polymorphism relative to divergence in genomic regions where the rate of crossing over per physical distance is restricted. We tested this prediction near the telomere of the Drosophila melanogaster and D. simulans X chromosome at two loci, erect wing (ewg) and suppressor of sable [su(s)]. Consistent with this prediction, polymorphism is reduced at both loci, while divergence is normal. The reduction is greater at ewg, the more distal of the two regions. Two models can be discriminated by comparing the observed site frequency spectra with those predicted by the models. The hitchhiking model predicts a skew toward rare variants in a sample, while the spectra under the background-selection model are similar to those of the neutral model of molecular evolution. Statistical tests of the fit to the predictions of these models require many sampled alleles and segregating sites. Thus we used SSCP and stratified DNA sequencing to cover a large number of randomly sampled alleles (approximately 50) from each of three populations. The result is a clear trend toward negative values of Tajima's D, indicating an excess of rare variants at ewg, the more distal of the two loci. One fixed difference among the populations and high FST values indicate strong population subdivision among the three populations at ewg. These results indicate genetic hitchhiking at ewg, in particular, geographically localized hitchhiking events within Africa. The reduction of polymorphism at su(s) combined with the excess of high-frequency variants in D. simulans is inconsistent with the hitchhiking and background-selection models.

Selection for highly biased amino acid frequency in the TolA cell envelope protein of Proteobacteria.

The bacterial cell envelope protein TolA functions to maintain the integrity of the cell membrane. This protein contains high levels of alanine and lysine that are used in the formation of alpha helices, which are required for normal protein function. The neutral model of molecular evolution predicts that amino acid composition and nucleotide composition are driven by the underlying GC content, as a result of mutation bias. However, this study shows that selection has acted to maintain high levels of alanine and lysine in the TolA protein of Proteobacteria, which in turn has biased nucleotide composition in the corresponding tolA gene.

Polymorphism patterns in two tightly linked developmental genes, Idgf1 and Idgf3, of Drosophila melanogaster.

A new developmental gene family, recently identified in D. melanogaster, has been called imaginal disc growth factors (IDGF) because the proteins promote growth of cell lineages derived from imaginal discs. These are the first genes reported that encode polypeptide factors with mitotic activity in invertebrates. Characteristics such as similar arrangement of introns and exons, small size, and different cytological localization make this family an excellent candidate for evolutionary studies. We focus on the loci Idgf1 and Idgf3, two genes that possess the most distinctive features. We examine the pattern of intra- and interspecific nucleotide variation in the sequences from 20 isogenic lines of D. melanogaster and sequences from D. simulans and D. yakuba. While MK, HKA, and Tajima's tests of neutrality fail to reject a neutral model of molecular evolution, Fu and Li's test with outgroup and McDonald's test suggest that balancing selection is modulating the evolution of the Idgf1 locus. The rate of recombination between the two loci is high enough to uncouple any linkage disequilibrium arising between Idgf1 and Idgf3, despite their close physical proximity.