Molecular Biology and Evolution - current issue

Associate

Mon, 09 Nov 2009 18:35:11 PST

Eboard

Mon, 09 Nov 2009 18:35:11 PST

Subcription

Mon, 09 Nov 2009 18:35:11 PST

Cover

Mon, 09 Nov 2009 18:35:11 PST

Selection to Maintain Paralogous Amino Acid Differences Under the Pressure of Gene Conversion in the Heat-Shock Protein Genes in Yeast

Takuno, S., Innan, H. — Mon, 09 Nov 2009 18:35:11 PST

A genome scan for the signatures of selection for paralogous functional amino acid differences was performed with yeast genomes. This recently developed method makes it possible to localize the target sites of selection under the pressure of gene conversion. We found that two gene pairs have strong signatures of selection. The two pairs of duplicated genes happened to be heat shock genes (Ssa1/ Ssa2 and Ssb1/Ssb2), which have similar protein structures to each other, although the amino acid sequence identity between Ssa and Ssb is not high (~60%). Interestingly, the two gene pairs exhibit signature of selection at almost identical positions within the substrate-binding domain β. Because this domain specifies the substrate polypeptides, it is presumed that functional divergence may be advantageous in this domain. Evolutionary analysis demonstrated that the observed divergence in the two gene pairs has been maintained in many yeast species independently, suggesting long-term operation of strong selection.

Transcription, Translation, and the Evolution of Specialists and Generalists

Zhong, S., Miller, S. P., Dykhuizen, D. E., Dean, A. M. — Mon, 09 Nov 2009 18:35:11 PST

We used DNA microarrays to measure transcription and iTRAQ 2D liquid chromatography-mass spectrometry/mass spectrometry (a mass-tag labeling proteomic technique) to measure protein expression in 14 strains of Escherichia coli adapted for hundreds of generations to growth-limiting concentrations of either lactulose, methylgalactoside, or a 72:28 mixture of the two. The two ancestors, TD2 and TD10, differ only in their lac operons and have similar transcription and protein expression profiles. Changes in transcription and protein expression are observed at 30–250 genes depending on the evolved strain. Lactulose specialists carry duplications of the lac operon and show increased transcription and translation at lac. Methylgalactoside specialists are galS^– and so constitutively transcribe and translate mgl, which encodes a transporter of methylgalactoside. However, there are two strains that carry lac duplications, are galS^–, and show increased transcription and translation at both operons. One is a generalist, the other a lactulose specialist. The generalist fails to sweep to fixation because its lac⁺, galS⁺ competitor expresses the csg adhesin and sticks to the chemostat wall, thereby preventing complete washout. Transcription and translation are sometimes decoupled. Lactulose-adapted strains show increased protein expression at fru, a fructose transporter, without evidence of increased transcription. This suggests that fructose, produced by the action of β-galactosidase on lactulose, may leach from cells before being recouped. Reduced expression, at "late" flagella genes and the constitutive gat operon, is an adaptation to starvation. A comparison with two other long-term evolution experiments suggests that certain aspects of adaptation are predictable, some are characteristic of an experimental system, whereas others seem erratic.

Molecular Evolution of GYPC: Evidence for Recent Structural Innovation and Positive Selection in Humans

Wilder, J. A., Hewett, E. K., Gansner, M. E. — Mon, 09 Nov 2009 18:35:11 PST

GYPC encodes two erythrocyte surface sialoglycoproteins in humans, glycophorin C and glycophorin D (GPC and GPD), via initiation of translation at two start codons on a single transcript. The malaria-causing parasite Plasmodium falciparum uses GPC as a means of invasion into the human red blood cell. Here, we examine the molecular evolution of GYPC among the Hominoidea (Greater and Lesser Apes) and also the pattern of polymorphism at the locus in a global human sample. We find an excess of nonsynonymous divergence among species that appears to be caused solely by accelerated evolution of GYPC in the human lineage. Moreover, we find that the ability of GYPC to encode both GPC and GPD is a uniquely human trait, caused by the evolution of the GPC start codon in the human lineage. The pattern of polymorphism among humans is consistent with a hitchhiking event at the locus, suggesting that positive natural selection affected GYPC in the relatively recent past. Because GPC is exploited by P. falciparum for invasion of the red blood cell, we hypothesize that selection for evasion of P. falciparum has caused accelerated evolution of GYPC in humans (relative to other primates) and that this positive selection has continued to act in the recent evolution of our species. These data suggest that malaria has played a powerful role in shaping molecules on the surface of the human red blood cell. In addition, our examination of GYPC reveals a novel mechanism of protein evolution: co-option of untranslated region (UTR) sequence following the formation of a new start codon. In the case of human GYPC, the ancestral protein (GPD) continues to be produced through leaky translation. Because leaky translation is a widespread phenomenon among genes and organisms, we suggest that co-option of UTR sequence may be an important source of protein innovation.

General Heterotachy and Distance Method Adjustments

Wu, J., Susko, E. — Mon, 09 Nov 2009 18:35:11 PST

Heterotachy is a general term to describe positions in a sequence that evolve at different rates in different lineages. Kolaczkowski and Thornton (2004. Performance of maximum parsimony and likelihood phylogenetics when evolution is heterogeneous. Nature 431:980–984.) recently described an intriguing heterotachy model that leads to topological bias for likelihood-based methods and parsimony methods. In this article, we show that heterotachy can generally be viewed as multivariate rates-across-sites variation, which can be described as randomly drawing rates (or branch lengths) from a multivariate distribution for each branch at each site. Motivated by this idea, we propose a pairwise alpha heterotachy adjustment model, which gives us much improved topological estimation in the settings by Kolaczkowski and Thornton (2004).

Phylogeny of the "Forgotten" Cellular Slime Mold, Fonticula alba, Reveals a Key Evolutionary Branch within Opisthokonta

Brown, M. W., Spiegel, F. W., Silberman, J. D. — Mon, 09 Nov 2009 18:35:11 PST

The shared ancestry between Fungi and animals has been unequivocally demonstrated by abundant molecular and morphological data for well over a decade. Along with the animals and Fungi, multiple protists have been placed in the supergroup Opisthokonta making it exceptionally diverse. In an effort to place the cellular slime mold Fonticula alba, an amoeboid protist with aggregative, multicellular fruiting, we sequenced five nuclear encoded genes; small subunit ribosomal RNA, actin, beta-tubulin, elongation factor 1-alpha, and the cytosolic isoform of heat shock protein 70 for phylogenetic analyses. Molecular trees demonstrate that Fonticula is an opisthokont that branches sister to filose amoebae in the genus Nuclearia. Fonticula plus Nuclearia are sister to Fungi. We propose a new name for this well-supported clade, Nucletmycea, incorporating Nuclearia, Fonticula, and Fungi. Fonticula represents the first example of a cellular slime mold morphology within Opisthokonta. Thus, there are four types of multicellularity in the supergroup—animal, fungal, colonial, and now aggregative. Our data indicate that multicellularity in Fonticula evolved independent of that found in the fungal and animal radiations. With the rapidly expanding sequence and genomic data becoming available from many opisthokont lineages, Fonticula may be fundamental to understanding opisthokont evolution as well as any possible commonalities involved with the evolution of multicellularity.

Hemocyanin Suggests a Close Relationship of Remipedia and Hexapoda

Ertas, B., von Reumont, B. M., Wagele, J.-W., Misof, B., Burmester, T. — Mon, 09 Nov 2009 18:35:11 PST

The Remipedia are enigmatic crustaceans from anchialine cave systems, first described only 30 years ago, whose phylogenetic affinities are as yet unresolved. Here we report the sequence of hemocyanin from Speleonectes tulumensis Yager, 1987 (Remipedia, Speleonectidae). This is the first proof of the presence of this type of respiratory protein in a crustacean taxon other than Malacostraca. Speleonectes tulumensis hemocyanin consists of multiple distinct (at least three) subunits (StuHc1–3; Hc, hemocyanin). Surprisingly, the sequences are most similar to hexapod hemocyanins. Phylogenetic analyses showed that the S. tulumensis hemocyanin subunits StuHc1 and StuHc3 associate with the type 1 hexapod hemocyanin subunits, whereas StuHc2 associates with the type 2 subunits of hexapods. Together, remipede and hexapod hemocyanins are in the sister-group position to the hemocyanins of malacostracan crustaceans. Hemocyanins provide no indication of a close relationship of Myriapoda and Hexapoda but support Pancrustacea (Crustacea + Hexapoda). Our results also suggest that Crustacea are paraphyletic and that Hexapoda may have evolved from a Remipedia-like ancestor. Thus, Remipedia occupy a key position for the understanding of the evolution of hexapods, which are and have been one of the world's most speciose lineage of animals.

A Phylogenomic Approach to Resolve the Basal Pterygote Divergence

Simon, S., Strauss, S., von Haeseler, A., Hadrys, H. — Mon, 09 Nov 2009 18:35:11 PST

One of the most fascinating Bauplan transitions in the animal kingdom was the invention of insect wings, a change that also contributed to the success and enormous diversity of this animal group. However, the origin of insect flight and the relationships of basal winged insect orders are still controversial. Three hypotheses have been proposed to explain the phylogeny of winged insects: 1) the traditional Palaeoptera hypothesis (Ephemeroptera + Odonata, Neoptera), 2) the Metapterygota hypothesis (Ephemeroptera, Odonata + Neoptera), and 3) the Chiastomyaria hypothesis (Odonata, Ephemeroptera + Neoptera). Neither phylogenetic analyses of single genes nor even multiple marker systems (e.g., molecular markers + morphological characters) have yet been able to conclusively resolve basal pterygote divergences. A possible explanation for the lack of resolution is that the divergences took place in the mid-Devonian within a short period of time and attempts to solve this problem have been confounded by the major challenge of finding molecular markers to accurately track these short ancient internodes. Although phylogenomic data are available for Neoptera and some wingless (apterygote) orders, they are lacking for the crucial Odonata and Ephemeroptera orders. We adopt a multigene approach including data from two new expressed sequence tag projects—from the orders Ephemeroptera (Baetis sp.) and Odonata (Ischnura elegans)—to evaluate the potential of phylogenomic analyses in clarifying this unresolved issue. We analyzed two data sets that differed in represented taxa, genes, and overall sequence lengths: maxspe (15 taxa, 125 genes, and 31,643 amino acid positions) and maxgen (8 taxa, 150 genes, and 42,541 amino acid positions). Maximum likelihood and Bayesian inference analyses both place the Odonata at the base of the winged insects. Furthermore, statistical hypotheses testing rejected both the Palaeoptera and the Metapterygota hypotheses. The comprehensive molecular data set developed here provides conclusive support for odonates as the most basal winged insect order (Chiastomyaria hypothesis). Data quality assessment indicates that proteins involved in cellular processes and signaling harbor the most informative phylogenetic signal.

Benchmarking Next-Generation Transcriptome Sequencing for Functional and Evolutionary Genomics

Mon, 09 Nov 2009 18:35:11 PST

Next-generation sequencing has opened the door to genomic analysis of nonmodel organisms. Technologies generating long-sequence reads (200–400 bp) are increasingly used in evolutionary studies of nonmodel organisms, but the short-sequence reads (30–50 bp) that can be produced at lower cost are thought to be of limited utility for de novo sequencing applications. Here, we tested this assumption by short-read sequencing the transcriptomes of the tropical disease vectors Aedes aegypti and Anopheles gambiae, for which complete genome sequences are available. Comparison of our results to the reference genomes allowed us to accurately evaluate the quantity, quality, and functional and evolutionary information content of our "test" data. We produced more than 0.7 billion nucleotides of sequenced data per species that assembled into more than 21,000 test contigs larger than 100 bp per species and covered ~27% of the Aedes reference transcriptome. Remarkably, the substitution error rate in the test contigs was ~0.25% per site, with very few indels or assembly errors. Test contigs of both species were enriched for genes involved in energy production and protein synthesis and underrepresented in genes involved in transcription and differentiation. Ortholog prediction using the test contigs was accurate across hundreds of millions of years of evolution. Our results demonstrate the considerable utility of short-read transcriptome sequencing for genomic studies of nonmodel organisms and suggest an approach for assessing the information content of next-generation data for evolutionary studies.

Signal Conflicts in the Phylogeny of the Primary Photosynthetic Eukaryotes

Deschamps, P., Moreira, D. — Mon, 09 Nov 2009 18:35:11 PST

It is widely accepted that the first photosynthetic eukaryotes arose from a single primary endosymbiosis of a cyanobacterium in a phagotrophic eukaryotic host, which led to the emergence of three major lineages: Chloroplastida (green algae and land plants), Rhodophyta, and Glaucophyta. For a long time, Glaucophyta have been thought to represent the earliest branch among them. However, recent massive phylogenomic analyses of nuclear genes have challenged this view, because most of them suggested a basal position of Rhodophyta, though with moderate statistical support. We have addressed this question by phylogenomic analysis of a large data set of 124 proteins transferred from the chloroplast to the nuclear genome of the three Archaeplastida lineages. In contrast to previous analyses, we found strong support for the basal emergence of the Chloroplastida and the sister-group relationship of Glaucophyta and Rhodophyta. Moreover, the reanalysis of chloroplast gene sequences using methods more robust against compositional and evolutionary rate biases sustained the same result. Finally, we observed that the basal position of Rhodophyta found in the phylogenies based on nuclear genes depended on the sampling of sequences used as outgroup. When eukaryotes supposed to have never had plastids (animals and fungi) were used, the analysis strongly supported the early emergence of Glaucophyta instead of Rhodophyta. Therefore, there is a conflicting signal between genes of different evolutionary origins supporting either the basal branching of Glaucophyta or of Chloroplastida within the Archaeplastida. This second possibility would agree with the existence of the subkingdom Biliphyta, joining Glaucophyta and Rhodophyta.

Targets of Balancing Selection in the Human Genome

Mon, 09 Nov 2009 18:35:11 PST

Balancing selection is potentially an important biological force for maintaining advantageous genetic diversity in populations, including variation that is responsible for long-term adaptation to the environment. By serving as a means to maintain genetic variation, it may be particularly relevant to maintaining phenotypic variation in natural populations. Nevertheless, its prevalence and specific targets in the human genome remain largely unknown. We have analyzed the patterns of diversity and divergence of 13,400 genes in two human populations using an unbiased single-nucleotide polymorphism data set, a genome-wide approach, and a method that incorporates demography in neutrality tests. We identified an unbiased catalog of genes with signatures of long-term balancing selection, which includes immunity genes as well as genes encoding keratins and membrane channels; the catalog also shows enrichment in functional categories involved in cellular structure. Patterns are mostly concordant in the two populations, with a small fraction of genes showing population-specific signatures of selection. Power considerations indicate that our findings represent a subset of all targets in the genome, suggesting that although balancing selection may not have an obvious impact on a large proportion of human genes, it is a key force affecting the evolution of a number of genes in humans.

Tracing the History of Goat Pastoralism: New Clues from Mitochondrial and Y Chromosome DNA in North Africa

Mon, 09 Nov 2009 18:35:11 PST

Valuable insights into the history of human populations have been obtained by studying the genetic composition of their domesticated species. Here we address some of the long-standing questions about the origin and subsequent movements of goat pastoralism in Northern Africa. We present the first study combining results from mitochondrial DNA (mtDNA) and Y chromosome loci for the genetic characterization of a domestic goat population. Our analyses indicate a remarkably high diversity of maternal and paternal lineages in a sample of indigenous goats from the northwestern fringe of the African continent. Median-joining networks and a multidimensional scaling of ours and almost 2000 published mtDNA sequences revealed a considerable genetic affinity between goat populations from the Maghreb (Northwest Africa) and the Near East. It has been previously shown that goats have a weak phylogeographic structure compatible with high levels of gene flow, as demonstrated by the worldwide dispersal of the predominant mtDNA haplogroup A. In contrast, our results revealed a strong correlation between genetic and geographical distances in 20 populations from different regions of the world. The distribution of Y chromosome haplotypes in Maghrebi goats indicates a common origin for goat patrilines in both Mediterranean coastal regions. Taken together, these results suggest that the colonization and subsequent dispersal of domestic goats in Northern Africa was influenced by the maritime diffusion throughout the Mediterranean Sea and its coastal regions of pastoralist societies whose economy included goat herding. Finally, we also detected traces of gene flow between goat populations from the Maghreb and the Iberian Peninsula corroborating evidence of past cultural and commercial contacts across the Strait of Gibraltar.

A Comprehensive Classification and Evolutionary Analysis of Plant Homeobox Genes

Mukherjee, K., Brocchieri, L., Burglin, T. R. — Mon, 09 Nov 2009 18:35:11 PST

The full complement of homeobox transcription factor sequences, including genes and pseudogenes, was determined from the analysis of 10 complete genomes from flowering plants, moss, Selaginella, unicellular green algae, and red algae. Our exhaustive genome-wide searches resulted in the discovery in each class of a greater number of homeobox genes than previously reported. All homeobox genes can be unambiguously classified by sequence evolutionary analysis into 14 distinct classes also characterized by conserved intron–exon structure and by unique codomain architectures. We identified many new genes belonging to previously defined classes (HD-ZIP I to IV, BEL, KNOX, PLINC, WOX). Other newly identified genes allowed us to characterize PHD, DDT, NDX, and LD genes as members of four new evolutionary classes and to define two additional classes, which we named SAWADEE and PINTOX. Our comprehensive analysis allowed us to identify several newly characterized conserved motifs, including novel zinc finger motifs in SAWADEE and DDT. Members of the BEL and KNOX classes were found in Chlorobionta (green plants) and in Rhodophyta. We found representatives of the DDT, WOX, and PINTOX classes only in green plants, including unicellular green algae, moss, and vascular plants. All 14 homeobox gene classes were represented in flowering plants, Selaginella, and moss, suggesting that they had already differentiated in the last common ancestor of moss and vascular plants.

A Broadscale Phylogenetic Analysis of Group II Intron RNAs and Intron-Encoded Reverse Transcriptases

Simon, D. M., Kelchner, S. A., Zimmerly, S. — Mon, 09 Nov 2009 18:35:11 PST

Group II introns are self-splicing RNAs that are frequently assumed to be the ancestors of spliceosomal introns. They are widely distributed in bacteria and are also found in organelles of plants, fungi, and protists. In this study, we present a broadscale phylogenetic analysis of group II introns using sequence data from both the conserved RNA structure and the intron-encoded reverse transcriptase (RT). Two similar phylogenies are estimated for the RT open reading frame (ORF), based on either amino acid or nucleotide sequence, whereas one phylogeny is produced for the RNA. In making these estimates, we confronted nearly all the classic challenges to phylogenetic inference, including positional saturation, base composition heterogeneity, short internodes with low support, and sensitivity to taxon sampling. Although the major lineages are well-defined, robust resolution of topology is not possible between these lineages. The approximately unbiased (AU) and Shimodaira–Hasegawa topology tests indicated that the RT ORF and RNA ribozyme data sets are in significant conflict under a variety of models, revealing the possibility of imperfect coevolution between group II introns and their intron-encoded ORFs. The high level of sequence divergence, large timescale, and limited number of alignable characters in our study are representative of many RTs and group I introns, and our results suggest that phylogenetic analyses of any of these sequences could suffer from the same sources of error and instability identified in this study.

Global Microsatellite Content Distinguishes Humans, Primates, Animals, and Plants

Mon, 09 Nov 2009 18:35:11 PST

Microsatellites are highly mutable, repetitive sequences commonly used as genetic markers, but they have never been studied en masse. Using a custom microarray to measure hybridization intensities of every possible repetitive nucleotide motif from 1-mers to 6-mers, we examined 25 genomes. Here, we show that global microsatellite content varies predictably by species, as measured by array hybridization signal intensities, correlating with established taxonomic relationships, and particular motifs are characteristic of one species versus another. For instance, hominid-specific microsatellite motifs were identified despite alignment of the human reference, Celera, and Venter genomic sequences indicating substantial variation (30–50%) among individuals. Differential microsatellite motifs were mainly associated with genes involved in developmental processes, whereas those found in intergenic regions exhibited no discernible pattern. This is the first description of a method for evaluating microsatellite content to classify individual genomes.

Evolution of 7SK RNA and Its Protein Partners in Metazoa

Mon, 09 Nov 2009 18:35:11 PST

7SK RNA is a key player in the regulation of polymerase II transcription. 7SK RNA was considered as a highly conserved vertebrate innovation. The discovery of poorly conserved homologs in several insects and lophotrochozoans, however, implies a much earlier evolutionary origin. The mechanism of 7SK function requires interaction with the proteins HEXIM and La-related protein 7. Here, we present a comprehensive computational analysis of these two proteins in metazoa, and we extend the collection of 7SK RNAs by several additional candidates. In particular, we describe 7SK homologs in Caenorhabditis species. Furthermore, we derive an improved secondary structure model of 7SK RNA, which shows that the structure is quite well-conserved across animal phyla despite the extreme divergence at sequence level.

Genic Incompatibilities in Two Hybrid Bacteriophages

Rokyta, D. R., Wichman, H. A. — Mon, 09 Nov 2009 18:35:11 PST

Horizontal gene transfer and recombination play a major role in microbial evolution and have been detected in diverse groups, including many of medical relevance such as HIV and dengue virus. In the absence of mechanistic barriers, the evolutionary success of a particular recombination event is determined by whether the recombinant genotype suffers a fitness cost through the disruption of favorable epistatic interactions within the genome, and if so, the extent to which this fitness cost might be mitigated by subsequent compensatory evolution. To investigate the importance of epistatic interactions between genes and the evolutionary viability of a homologous recombination event between diverged ancestral genotypes, we constructed two recombinant microvirid bacteriophages by exchanging their alleles of the gene encoding the coat protein. The coding sequences for this gene differ by approximately 8% at the amino acid level and were interchanged between two ancestral phages related to X174 and well adapted to their culture conditions. Because the recombinant phages showed drastically reduced fitnesses, we further explored their evolutionary viability by subjecting replicate lines of each of them to selection. We found that all four lineages achieved fitnesses commensurate with ancestral fitnesses in as few as 60 generations, and on average, the first substitution accounted for more than half of the total fitness recovery. Fitness recovery required three to five substitutions in each lineage, and overall eight of the nine essential phage genes were involved, suggesting extensive epistatic interactions throughout the genome. Interestingly, the proteins with the most extensive and apparent physical interactions with the exchanged protein in the viral capsid did not appear to have much of a role in fitness recovery. This result appears to be a consequence of the conservation of the amino acid residues involved in the interactions. It suggests that strong epistatic interactions are less important than weaker, transient ones in producing genic incompatibilities because they preclude variability in the interacting regions of the proteins.

Novel Internal Regions of Fluorescent Proteins Undergo Divergent Evolutionary Patterns

Mon, 09 Nov 2009 18:35:11 PST

Over the past decade, fluorescent proteins (FPs) have become ubiquitous tools in biological research. Yet, little is known about the natural function or evolution of this superfamily of proteins that originate from marine organisms. Using molecular phylogenetic analyses of 102 naturally occurring cyan fluorescent proteins, green fluorescent proteins, red fluorescent proteins, as well as the nonfluorescent (purple-blue) protein sequences (including new FPs from Lizard Island, Australia) derived from organisms with known geographic origin, we show that FPs consist of two distinct and novel regions that have evolved under opposite and sharply divergent evolutionary pressures. A central region is highly conserved, and although it contains the residues that form the chromophore, its evolution does not track with fluorescent color and evolves independently from the rest of the protein. By contrast, the regions enclosing this central region are under strong positive selection pressure to vary its sequence and yet segregate well with fluorescence color emission. We did not find a significant correlation between geographic location of the organism from which the FP was isolated and molecular evolution of the protein. These results define for the first time two distinct regions based on evolution for this highly compact protein. The findings have implications for more sophisticated bioengineering of this molecule as well as studies directed toward understanding the natural function of FPs.

mtDNA Data Indicate a Single Origin for Dogs South of Yangtze River, Less Than 16,300 Years Ago, from Numerous Wolves

Mon, 09 Nov 2009 18:35:11 PST

There is no generally accepted picture of where, when, and how the domestic dog originated. Previous studies of mitochondrial DNA (mtDNA) have failed to establish the time and precise place of origin because of lack of phylogenetic resolution in the so far studied control region (CR), and inadequate sampling. We therefore analyzed entire mitochondrial genomes for 169 dogs to obtain maximal phylogenetic resolution and the CR for 1,543 dogs across the Old World for a comprehensive picture of geographical diversity. Hereby, a detailed picture of the origins of the dog can for the first time be suggested. We obtained evidence that the dog has a single origin in time and space and an estimation of the time of origin, number of founders, and approximate region, which also gives potential clues about the human culture involved. The analyses showed that dogs universally share a common homogenous gene pool containing 10 major haplogroups. However, the full range of genetic diversity, all 10 haplogroups, was found only in southeastern Asia south of Yangtze River, and diversity decreased following a gradient across Eurasia, through seven haplogroups in Central China and five in North China and Southwest (SW)Asia, down to only four haplogroups in Europe. The mean sequence distance to ancestral haplotypes indicates an origin 5,400–16,300 years ago (ya) from at least 51 female wolf founders. These results indicate that the domestic dog originated in southern China less than 16,300 ya, from several hundred wolves. The place and time coincide approximately with the origin of rice agriculture, suggesting that the dogs may have originated among sedentary hunter-gatherers or early farmers, and the numerous founders indicate that wolf taming was an important culture trait.

Evolutionary Trajectories of Primate Genes Involved in HIV Pathogenesis

Mon, 09 Nov 2009 18:35:11 PST

The current availability of five complete genomes of different primate species allows the analysis of genetic divergence over the last 40 million years of evolution. We hypothesized that the interspecies differences observed in susceptibility to HIV-1 would be influenced by the long-range selective pressures on host genes associated with HIV-1 pathogenesis. We established a list of human genes (n = 140) proposed to be involved in HIV-1 biology and pathogenesis and a control set of 100 random genes. We retrieved the orthologous genes from the genome of humans and of four nonhuman primates (Pan troglodytes, Pongo pygmaeus abeli, Macaca mulatta, and Callithrix jacchus) and analyzed the nucleotide substitution patterns of this data set using codon-based maximum likelihood procedures. In addition, we evaluated whether the candidate genes have been targets of recent positive selection in humans by analyzing HapMap Phase 2 single-nucleotide polymorphisms genotyped in a region centered on each candidate gene. A total of 1,064 sequences were used for the analyses. Similar median K_A/K_S values were estimated for the set of genes involved in HIV-1 pathogenesis and for control genes, 0.19 and 0.15, respectively. However, genes of the innate immunity had median values of 0.37 (P value = 0.0001, compared with control genes), and genes of intrinsic cellular defense had K_A/K_S values around or greater than 1.0 (P value = 0.0002). Detailed assessment allowed the identification of residues under positive selection in 13 proteins: AKT1, APOBEC3G, APOBEC3H, CD4, DEFB1, GML, IL4, IL8RA, L-SIGN/CLEC4M, PTPRC/CD45, Tetherin/BST2, TLR7, and TRIM5. A number of those residues are relevant for HIV-1 biology. The set of 140 genes involved in HIV-1 pathogenesis did not show a significant enrichment in signals of recent positive selection in humans (intraspecies selection). However, we identified within or near these genes 24 polymorphisms showing strong signatures of recent positive selection. Interestingly, the DEFB1 gene presented signatures of both interspecies positive selection in primates and intraspecies recent positive selection in humans. The systematic assessment of long-acting selective pressures on primate genomes is a useful tool to extend our understanding of genetic variation influencing contemporary susceptibility to HIV-1.

Lineage-Specific Adaptive Evolution of the Centromeric Protein CENH3 in Diploid and Allotetraploid Oryza Species

Hirsch, C. D., Wu, Y., Yan, H., Jiang, J. — Mon, 09 Nov 2009 18:35:11 PST

Centromeres in eukaryotic species are defined by the presence of a centromere-specific histone H3 variant, CENH3. CENH3 plays a key role in recruiting other centromeric proteins; thus, it is the central component in kinetochore formation and centromere function. The CENH3 proteins in several plant and animal species were found to be under positive selection, which was hypothesized to respond to the rapid changing of the repetitive DNA sequences associated with the centromeres. Here, we report the expression and evolution of the CenH3 genes in two allotetraploid rice species as well as their representative diploid progenitor species. Both copies of the CenH3 genes were transcribed in the two allotetraploid species and showed a nonpreferential expression pattern. Contrasting positive and stabilizing selection of the CenH3 genes was associated with different diploid Oryza species. This lineage-specific adaptive evolution of CENH3 was maintained in the two allotetraploid species. Thus, we demonstrate that the allopolyploidization events did not alter the expression or evolutionary patterns of the CenH3 genes in the Oryza species.