Evidence for Pervasive Adaptive Protein Evolution in Wild Mice

The relative contributions of neutral and adaptive substitutions to molecular evolution has been one of the most controversial issues in evolutionary biology for more than 40 years. The analysis of within-species nucleotide polymorphism and between-species divergence data supports a widespread role for adaptive protein evolution in certain taxa. For example, estimates of the proportion of adaptive amino acid substitutions (alpha) are 50% or more in enteric bacteria and Drosophila. In contrast, recent estimates of alpha for hominids have been at most 13%. Here, we estimate alpha for protein sequences of murid rodents based on nucleotide polymorphism data from multiple genes in a population of the house mouse subspecies Mus musculus castaneus, which inhabits the ancestral range of the Mus species complex and nucleotide divergence between M. m. castaneus and M. famulus or the rat. We estimate that 57% of amino acid substitutions in murids have been driven by positive selection. Hominids, therefore, are exceptional in having low apparent levels of adaptive protein evolution. The high frequency of adaptive amino acid substitutions in wild mice is consistent with their large effective population size, leading to effective natural selection at the molecular level. Effective natural selection also manifests itself as a paucity of effectively neutral nonsynonymous mutations in M. m. castaneus compared to humans

Molecular evolution of the membrane associated progesterone receptor in the Brachionus plicatilis (Rotifera, Monogononta) species complex

Author Posting. © Springer, 2010. This is the author's version of the work. It is posted here by permission of Springer for personal use, not for redistribution. The definitive version was published in Hydrobiologia 662 (2011): 99-106, doi:10.1007/s10750-010-0484-4.Many studies have investigated physiological roles of the membrane associated progesterone receptor (MAPR), but little is known of its evolution. Marked variations in response to exogenous progesterone have been reported for four brachionid rotifer species, suggesting differences in progesterone signaling and reception. Here we report sequence variation for the MAPR gene in the Brachionus plicatilis species complex. Phylogenetic analysis of this receptor is compared with relatedness based on cytochrome c oxidase subunit 1 sequences. Nonsynonymous to synonymous site substitution rate ratios, amino acid divergence, and variations in predicted phosphorylation sites are examined to assess evolution of the MAPR among brachionid clades.National Science Foundation grant BE/GenEn MCB-0412674E to TWS and DMW, and an NSF IGERT fellowship to HAS under DGE 0114400, supported this work

The Glyceraldehyde-3-Phosphate Dehydrogenase and the Small GTPase Rab 2 Are Crucial for Brucella Replication

The intracellular pathogen Brucella abortus survives and replicates inside host cells within an endoplasmic reticulum (ER)-derived replicative organelle named the “Brucella-containing vacuole” (BCV). Here, we developed a subcellular fractionation method to isolate BCVs and characterize for the first time the protein composition of its replicative niche. After identification of BCV membrane proteins by 2 dimensional (2D) gel electrophoresis and mass spectrometry, we focused on two eukaryotic proteins: the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and the small GTPase Rab 2 recruited to the vacuolar membrane of Brucella. These proteins were previously described to localize on vesicular and tubular clusters (VTC) and to regulate the VTC membrane traffic between the endoplasmic reticulum (ER) and the Golgi. Inhibition of either GAPDH or Rab 2 expression by small interfering RNA strongly inhibited B. abortus replication. Consistent with this result, inhibition of other partners of GAPDH and Rab 2, such as COPI and PKC ι, reduced B. abortus replication. Furthermore, blockage of Rab 2 GTPase in a GDP-locked form also inhibited B. abortus replication. Bacteria did not fuse with the ER and instead remained in lysosomal-associated membrane vacuoles. These results reveal an essential role for GAPDH and the small GTPase Rab 2 in B. abortus virulence within host cells

Genome Features of “Dark-Fly”, a Drosophila Line Reared Long-Term in a Dark Environment

Organisms are remarkably adapted to diverse environments by specialized metabolisms, morphology, or behaviors. To address the molecular mechanisms underlying environmental adaptation, we have utilized a Drosophila melanogaster line, termed “Dark-fly”, which has been maintained in constant dark conditions for 57 years (1400 generations). We found that Dark-fly exhibited higher fecundity in dark than in light conditions, indicating that Dark-fly possesses some traits advantageous in darkness. Using next-generation sequencing technology, we determined the whole genome sequence of Dark-fly and identified approximately 220,000 single nucleotide polymorphisms (SNPs) and 4,700 insertions or deletions (InDels) in the Dark-fly genome compared to the genome of the Oregon-R-S strain, a control strain. 1.8% of SNPs were classified as non-synonymous SNPs (nsSNPs: i.e., they alter the amino acid sequence of gene products). Among them, we detected 28 nonsense mutations (i.e., they produce a stop codon in the protein sequence) in the Dark-fly genome. These included genes encoding an olfactory receptor and a light receptor. We also searched runs of homozygosity (ROH) regions as putative regions selected during the population history, and found 21 ROH regions in the Dark-fly genome. We identified 241 genes carrying nsSNPs or InDels in the ROH regions. These include a cluster of alpha-esterase genes that are involved in detoxification processes. Furthermore, analysis of structural variants in the Dark-fly genome showed the deletion of a gene related to fatty acid metabolism. Our results revealed unique features of the Dark-fly genome and provided a list of potential candidate genes involved in environmental adaptation

Limits to the Rate of Adaptive Substitution in Sexual Populations

In large populations, many beneficial mutations may be simultaneously available and may compete with one another, slowing adaptation. By finding the probability of fixation of a favorable allele in a simple model of a haploid sexual population, we find limits to the rate of adaptive substitution, , that depend on simple parameter combinations. When variance in fitness is low and linkage is loose, the baseline rate of substitution is , where is the population size, is the rate of beneficial mutations per genome, and is their mean selective advantage. Heritable variance in log fitness due to unlinked loci reduces by under polygamy and under monogamy. With a linear genetic map of length Morgans, interference is yet stronger. We use a scaling argument to show that the density of adaptive substitutions depends on , , , and only through the baseline density: . Under the approximation that the interference due to different sweeps adds up, we show that , implying that interference prevents the rate of adaptive substitution from exceeding one per centimorgan per 200 generations. Simulations and numerical calculations confirm the scaling argument and confirm the additive approximation for ; for higher , the rate of adaptation grows above , but only very slowly. We also consider the effect of sweeps on neutral diversity and show that, while even occasional sweeps can greatly reduce neutral diversity, this effect saturates as sweeps become more common—diversity can be maintained even in populations experiencing very strong interference. Our results indicate that for some organisms the rate of adaptive substitution may be primarily recombination-limited, depending only weakly on the mutation supply and the strength of selection