Molecular Biology and Evolution - current issue

Convergent Evolution of Clustering of Iroquois Homeobox Genes across Metazoans

Irimia, M., Maeso, I., Garcia-Fernandez, J. — 2008-07-15

Vertebrate and Drosophila Iroquois genes are organized in clusters of 3 genes sharing blocks of conserved regulatory sequences. Here, we report a 3-gene cluster in the basal, preduplicative chordate amphioxus. Surprisingly, however, the origin of the amphioxus cluster is independent of those in vertebrates and drosophilids. Investigation of genomic organization of Iroquois genes in other 17 metazoan genomes revealed a fourth independent 3-gene cluster organization in polychaetes, as well as additional 2- and 4-gene clusters in other clades, in one of the most striking examples of convergence in genomic organization described so far. The recurrent independent evolution of Iroquois clusters suggests a functional importance of this organization for these genes, perhaps related to the sharing of regulatory elements. Consistent with this, comparative analysis of genomic regions flanking the 3 amphioxus Irx genes revealed several blocks of sequences, conserved for at least 100 Myr. Finally, we discuss the possible causes and implications of the convergent evolution of this genomic and regulatory organization throughout metazoans.

The Guanine Nucleotide Exchange Factors Sec2 and PRONE: Candidate Synapomorphies for the Opisthokonta and the Archaeplastida

Elias, M. — 2008-07-15

Although recent multigene phylogenetic analyses support close relationship of Metazoa and Fungi (the eukaryotic supergroup Opisthokonta) and monophyly of eukaryotes with the primary plastid, that is, Chloroplastida, Rhodophyta, and Glaucophyta (the supergroup Archaeplastida or Plantae), some authors still challenge this scheme. I found that 2 particular types of guanine nucleotide exchange factors (GEFs, i.e., cofactors of GTPases) might provide a new piece of evidence to resolve this controversy. An exhaustive analysis of available sequence data revealed that Sec2-related proteins, known to serve as GEF for exocytic GTPases of the Rab8/Sec4 subfamily, are restricted to opisthokonts, whereas proteins with the PRONE domain, recently described as novel plant-specific GEFs for RHO family GTPases, occur only in Chloroplastida and Rhodophyta. The results thus point to possible evolutionary innovations in the exocytic apparatus of the ancestral opisthokonts and reveal the probably first plastid-independent trait (i.e., a unique mode of RHO GTPase regulation) exclusive for Chloroplastida + Rhodophyta, further supporting monophyly of these 2 groups.

Natural Selection against Protein Aggregation on Self-Interacting and Essential Proteins in Yeast, Fly, and Worm

Chen, Y., Dokholyan, N. V. — 2008-07-15

Protein aggregation is the phenomenon of protein self-association potentially leading to detrimental effects on physiology, which is closely related to numerous human diseases such as Alzheimer's and Parkinson's disease. Despite progress in understanding the mechanism of protein aggregation, how natural selection against protein aggregation acts on subunits of protein complexes and on proteins with different contributions to organism fitness remains largely unknown. Here, we perform a proteome-wide analysis by using an experimentally validated algorithm TANGO and utilizing sequence, interactomic and phenotype-based functional genomic data from yeast, fly, and nematode. We find that proteins that are capable of forming homooligomeric complex have lower aggregation propensity compared with proteins that do not function as homooligomer. Further, proteins that are essential to the fitness of an organism have lower aggregation propensity compared with nonessential ones. Our finding suggests that the selection force against protein aggregation acts across different hierarchies of biological system.

Genetic Exchange versus Genetic Differentiation in a Medium-Sized Inversion of Drosophila: The A2/Ast Arrangements of Drosophila subobscura

Nobrega, C., Khadem, M., Aguade, M., Segarra, C. — 2008-07-15

Chromosomal inversion polymorphism affects nucleotide variation at loci associated with inversions. In Drosophila subobscura, a species with a rich chromosomal inversion polymorphism and the largest recombinational map so far reported in the Drosophila genus, extensive genetic structure of nucleotide variation was detected in the segment affected by the O₃ inversion, a moderately sized inversion at Muller's element E. Indeed, a strong genetic differentiation all over O₃ and no evidence of a higher genetic exchange in the center of the inversion than at breakpoints were detected. In order to ascertain, whether other polymorphic and differently sized inversions of D. subobscura also exhibited a strong genetic structure, nucleotide variation in 5 gene regions (P236, P275, P150, Sxl, and P125) located along the A₂ inversion was analyzed in A_st and A₂ chromosomes of D. subobscura. A₂ is a medium-sized inversion at Muller's element A and forms a single inversion loop in heterokaryotypes. The lower level of variation in A₂ relative to A_st and the significant excess of low-frequency variants at polymorphic sites indicate that nucleotide variation at A₂ is not at mutation–drift equilibrium. The closest region to an inversion breakpoint, P236, exhibits the highest level of genetic differentiation (F_ST) and of linkage disequilibrium (LD) between arrangements and variants at nucleotide polymorphic sites. The remaining 4 regions show a higher level of genetic exchange between A₂ and A_st chromosomes than P236, as revealed by F_ST and LD estimates. However, significant genetic differentiation between the A_st and A₂ arrangements was detected not only at P236 but also in the other 4 regions separated from the nearest breakpoint by 1.2–2.9 Mb. Therefore, the extent of genetic exchange between arrangements has not been high enough to homogenize nucleotide variation in the center of the A₂ inversion. A₂ can be considered a typical successful inversion of D. subobscura according to its relative length. Chromosomal inversion polymorphism of D. subobscura might thus cause the genome of this species to be highly structured and to harbor different gene pools that might contribute to maintain adaptations to particular environments.

Influence of the Mutant Spectrum in Viral Evolution: Focused Selection of Antigenic Variants in a Reconstructed Viral Quasispecies

Martin, V., Domingo, E. — 2008-07-15

RNA viruses replicate as complex mutant distributions termed viral quasispecies. Despite this, studies on virus populations subjected to positive selection have generally been performed and analyzed as if the viral population consisted of a defined genomic nucleotide sequence; such a simplification may not reflect accurately the molecular events underlying the selection process. In the present study, we have reconstructed a foot-and-mouth disease virus quasispecies with multiple, low-frequency, genetically distinguishable mutants that can escape neutralization by a monoclonal antibody. Some of the mutants included an amino acid substitution that affected an integrin recognition motif that overlaps with the antibody-binding site, whereas other mutants included an amino acid substitution that affected antibody binding but not integrin recognition. We have monitored consensus and clonal nucleotide sequences of populations passaged either in the absence or the presence of the neutralizing antibody. In both cases, the populations focused toward a specific mutant that was surrounded by a cloud of mutants with different antigenic and cell recognition specificities. In the absence of antibody selection, an antigenic variant that maintained integrin recognition became dominant, but the mutant cloud included as one of its minority components a variant with altered integrin recognition. Conversely, in the presence of antibody selection, a variant with altered integrin recognition motif became dominant, but it was surrounded by a cloud of antigenic variants that maintained integrin recognition. The results have documented that a mutant spectrum can exert an influence on a viral population subjected to a sustained positive selection pressure and have unveiled a mechanism of antigenic flexibility in viral populations, consisting in the presence in the selected quasispecies of mutants with different antigenic and cell recognition specificities.

Coexpression of Linked Genes in Mammalian Genomes Is Generally Disadvantageous

Liao, B.-Y., Zhang, J. — 2008-07-15

Similarity in gene expression pattern between closely linked genes is known in several eukaryotes. Two models have been proposed to explain the presence of such coexpression patterns. The adaptive model assumes that coexpression is advantageous and is established by relocation of initially unlinked but coexpressed genes, whereas the neutral model asserts that coexpression is a type of leaky expression due to similar expressional environments of linked genes, but is neither advantageous nor detrimental. However, these models are incompatible with several empirical observations. Here, we propose that coexpression of linked genes is a form of transcriptional interference that is disadvantageous to the organism. We show that even distantly linked genes that are tens of megabases away exhibit significant coexpression in the human genome. However, the linkage is more likely to be broken during evolution between genes of high coexpression than those of low coexpression and the breakage of linkage reduces gene coexpression. These results support our hypothesis that coexpression of linked genes in mammalian genomes is generally disadvantageous, implying that many mammalian genes may never reach their optimal expression pattern due to the interference of their genomic environment and that such transcriptional interference may be a force promoting recurrent relocation of genes in the genome.

Substitution of the Gene for Chloroplast RPS16 Was Assisted by Generation of a Dual Targeting Signal

2008-07-15

Organelle (mitochondria and chloroplasts in plants) genomes lost a large number of genes after endosymbiosis occurred. Even after this major gene loss, organelle genomes still lose their own genes, even those that are essential, via gene transfer to the nucleus and gene substitution of either different organelle origin or de novo genes. Gene transfer and substitution events are important processes in the evolution of the eukaryotic cell. Gene loss is an ongoing process in the mitochondria and chloroplasts of higher plants. The gene for ribosomal protein S16 (rps16) is encoded in the chloroplast genome of most higher plants but not in Medicago truncatula and Populus alba. Here, we show that these 2 species have compensated for loss of the rps16 from the chloroplast genome by having a mitochondrial rps16 that can target the chloroplasts as well as mitochondria. Furthermore, in Arabidopsis thaliana, Lycopersicon esculentum, and Oryza sativa, whose chloroplast genomes encode the rps16, we show that the product of the mitochondrial rps16 has dual targeting ability. These results suggest that the dual targeting of RPS16 to the mitochondria and chloroplasts emerged before the divergence of monocots and dicots (140–150 MYA). The gene substitution of the chloroplast rps16 by the nuclear-encoded rps16 in higher plants is the first report about ongoing gene substitution by dual targeting and provides evidence for an intermediate stage in the formation of this heterogeneous organelle.

How Well Does the HoT Score Reflect Sequence Alignment Accuracy?

Hall, B. G. — 2008-07-15

Multiple sequence alignment is an essential tool in many areas of biological research, and the accuracy of an alignment can strongly affect the accuracy of a downstream application such as phylogenetic analysis, identification of functional motifs, or polymerase chain reaction primer design. The heads or tails (HoT) method (Landan G, Graur D. 2007. Heads or tails: a simple reliability check for multiple sequence alignments. Mol Biol Evol. 24:1380–1383.) assesses the consistency of an alignment by comparing the alignment of a set of sequences with the alignment of the same set of sequences written in reverse order. This study shows that HoT scores and the alignment accuracies are positively correlated, so alignments with higher HoT scores are preferable. However, HoT scores are overestimates of alignment accuracy in general, with the extent of overestimation depending on the method used for multiple sequence alignment.

Duplication and Functional Diversification of HAP3 Genes Leading to the Origin of the Seed-Developmental Regulatory Gene, LEAFY COTYLEDON1 (LEC1), in Nonseed Plant Genomes

Xie, Z., Li, X., Glover, B. J., Bai, S., Rao, G.-Y., Luo, J., Yang, J. — 2008-07-15

The HAP3 gene encodes a subunit of the CCAAT-box-binding factor (CBF), a highly conserved trimeric activator that recognizes and binds the ubiquitous CCAAT promoter element with high affinity. Two types of HAP3 gene have been identified in plant genomes. The LEAFY COTYLEDON1 (LEC1)–type HAP3 genes encode a functionally specialized subunit of CBF, which is expressed specifically in developing seeds. In contrast, most non–LEC1-type HAP3 genes are expressed in various tissues. It has been proposed that the LEC1-type HAP3 genes originated from the duplication and functional divergence of non–LEC1-type HAP3 genes. However, it is not yet known when this duplication event took place or whether the LEC1-type HAP3 genes appeared at the same time as the origin of seed plants. Here we describe a comprehensive comparison of the duplication patterns of HAP3 genes in different plant genomes. We recognize a major expansion of the HAP3 gene family accompanying the origin and early diversification of land plants and postulate that retrotransposition and other mechanisms of gene duplication have been involved in the expansion of the plant HAP3 gene family. We provide evidence that the LEC1-type HAP3 genes originated in nonseed vascular plant genomes and demonstrate that they are inductively expressed under drought stress in nonseed plants. These genes, however, were recruited to a novel regulatory network in the early stages of seed plant evolution and steadily expressed during seed development and maturation.

Distinct Evolutionary Patterns between Chemoreceptors of 2 Vertebrate Olfactory Systems and the Differential Tuning Hypothesis

Grus, W. E., Zhang, J. — 2008-07-15

Most tetrapod vertebrates have 2 olfactory systems, the main olfactory system (MOS) and the vomeronasal system (VNS). According to the dual olfactory hypothesis, the MOS detects environmental odorants, whereas the VNS recognizes intraspecific pheromonal cues. However, this strict functional distinction has been blurred by recent reports that both systems can perceive both types of signals. Studies of a limited number of receptors suggest that MOS receptors are broadly tuned generalists, whereas VNS receptors are narrowly tuned specialists. However, whether this distinction applies to all MOS and VNS receptors remains unknown. The differential tuning hypothesis predicts that generalist MOS receptors detect an overlapping set of ligands and thus are more likely to be conserved over evolutionary time than specialist VNS receptors, which would evolve in a more lineage-specific manner. Here we test this prediction for all olfactory chemoreceptors by examining the evolutionary patterns of MOS-expressed odorant receptors (ORs) and trace amine–associated receptors (TAARs) and VNS-expressed vomeronasal type 1 receptors (V1Rs) and vomeronasal type 2 receptors (V2Rs) in 7 tetrapods (mouse, rat, dog, opossum, platypus, chicken, and frog). The phylogenies of V1Rs and V2Rs show abundant lineage-specific gene gains/losses and virtually no one-to-one orthologs between species. Opposite patterns are found for ORs and TAARs. Analysis of functional data and ligand-binding sites of ORs confirms that paralogous chemoreceptors are more likely than orthologs to have different ligands and that functional divergence between paralogous chemoreceptors is established relatively quickly following gene duplication. Together, these results strongly suggest that the functional profile of the VNS chemoreceptor repertoire evolves much faster than that of the MOS chemoreceptor repertoire and that the differential tuning hypothesis applies to the majority, if not all, of MOS and VNS receptors.

DNA Methylation and Structural and Functional Bimodality of Vertebrate Promoters

Elango, N., Yi, S. V. — 2008-07-15

Human promoters divide into 2 classes, the low CpG (LCG) and the high CpG (HCG), based on their CpG dinucleotide content. The LCG class of promoters is hypermethylated and is associated with tissue-specific genes, whereas the HCG class is hypomethylated and associated with broadly expressed genes. By analyzing several chordate genomes separated for hundreds of millions of years, here we show that the divide between low CpG and high CpG promoters is conserved in several distantly related vertebrate taxa (including human, chicken, frog, lizard, and fish) but not in close invertebrate outgroups (sea squirts). Furthermore, LCG and HCG promoters are distinctively associated with tissue-specific and broadly expressed genes in these distantly related vertebrate taxa. Our results indicate that the function of DNA methylation on gene expression is conserved across these vertebrate taxa and suggest that the 2 classes of promoters have evolved early in vertebrate evolution, as a consequence of the advent of global DNA methylation.

Long-Term Balancing Selection at the West Nile Virus Resistance Gene, Oas1b, Maintains Transspecific Polymorphisms in the House Mouse

Ferguson, W., Dvora, S., Gallo, J., Orth, A., Boissinot, S. — 2008-07-15

Oligoadenylate synthetases (OASs) are interferon-inducible enzymes that participate in the first line of defense against a wide range of viral infection. Recent studies have determined that Oas1b, a member of the OAS gene family in the house mouse (Mus musculus), provides specific protection against flavivirus infection (e.g., West Nile virus, dengue fever virus, and yellow fever virus). We characterized the nucleotide sequence variation in coding and noncoding regions of the Oas1b gene for a large number of wild-derived strains of M. musculus and related species. Our sequence analyses determined that this gene is one of the most polymorphic genes ever described in any mammal. The level of variation in noncoding regions of Oas1b is an order of magnitude higher than the level reported for other regions of the mouse genome and is significantly different from the level of intraspecific variation expected under neutrality. Furthermore, a phylogenetic analysis of intronic sequences demonstrated that Oas1b alleles are ancient and that their divergence predates several speciation events, resulting in transspecific polymorphisms. The amino acid sequence of Oas1b is also extremely variable, with 1 out of 7 amino acid positions being polymorphic within M. musculus. Oas1b alleles are comparatively more divergent at synonymous positions than most autosomal genes and the ratio of nonsynonymous to synonymous substitution is remarkably high, suggesting that positive selection has been acting on Oas1b. The ancestry of Oas1b polymorphisms and the high level of amino acid polymorphisms strongly suggest that the allelic variation at Oas1b has been maintained in mouse populations by long-term balancing selection.

The Deep Archaeal Roots of Eukaryotes

Yutin, N., Makarova, K. S., Mekhedov, S. L., Wolf, Y. I., Koonin, E. V. — 2008-07-15

The set of conserved eukaryotic protein-coding genes includes distinct subsets one of which appears to be most closely related to and, by inference, derived from archaea, whereas another one appears to be of bacterial, possibly, endosymbiotic origin. The "archaeal" genes of eukaryotes, primarily, encode components of information-processing systems, whereas the "bacterial" genes are predominantly operational. The precise nature of the archaeo–eukaryotic relationship remains uncertain, and it has been variously argued that eukaryotic informational genes evolved from the homologous genes of Euryarchaeota or Crenarchaeota (the major branches of extant archaea) or that the origin of eukaryotes lies outside the known diversity of archaea. We describe a comprehensive set of 355 eukaryotic genes of apparent archaeal origin identified through ortholog detection and phylogenetic analysis. Phylogenetic hypothesis testing using constrained trees, combined with a systematic search for shared derived characters in the form of homologous inserts in conserved proteins, indicate that, for the majority of these genes, the preferred tree topology is one with the eukaryotic branch placed outside the extant diversity of archaea although small subsets of genes show crenarchaeal and euryarchaeal affinities. Thus, the archaeal genes in eukaryotes appear to descend from a distinct, ancient, and otherwise uncharacterized archaeal lineage that acquired some euryarchaeal and crenarchaeal genes via early horizontal gene transfer.

Overdispersion of the Molecular Clock: Temporal Variation of Gene-Specific Substitution Rates in Drosophila

Bedford, T., Hartl, D. L. — 2008-07-15

Simple models of molecular evolution assume that sequences evolve by a Poisson process in which nucleotide or amino acid substitutions occur as rare independent events. In these models, the expected ratio of the variance to the mean of substitution counts equals 1, and substitution processes with a ratio greater than 1 are called overdispersed. Comparing the genomes of 10 closely related species of Drosophila, we extend earlier evidence for overdispersion in amino acid replacements as well as in four-fold synonymous substitutions. The observed deviation from the Poisson expectation can be described as a linear function of the rate at which substitutions occur on a phylogeny, which implies that deviations from the Poisson expectation arise from gene-specific temporal variation in substitution rates. Amino acid sequences show greater temporal variation in substitution rates than do four-fold synonymous sequences. Our findings provide a general phenomenological framework for understanding overdispersion in the molecular clock. Also, the presence of substantial variation in gene-specific substitution rates has broad implications for work in phylogeny reconstruction and evolutionary rate estimation.

Effects of X-Linkage and Sex-Biased Gene Expression on the Rate of Adaptive Protein Evolution in Drosophila

Baines, J. F., Sawyer, S. A., Hartl, D. L., Parsch, J. — 2008-07-15

Patterns of polymorphism and divergence in Drosophila protein-coding genes suggest that a considerable fraction of amino acid differences between species can be attributed to positive selection and that genes with sex-biased expression, that is, those expressed predominantly in one sex, have especially high rates of adaptive evolution. Previous studies, however, have been restricted to autosomal sex-biased genes and, thus, do not provide a complete picture of the evolutionary forces acting on sex-biased genes across the genome. To determine the effects of X-linkage on sex-biased gene evolution, we surveyed DNA sequence polymorphism and divergence in 45 X-linked genes, including 17 with male-biased expression, 13 with female-biased expression, and 15 with equal expression in the 2 sexes. Using both single- and multilocus tests for selection, we found evidence for adaptive evolution in both groups of sex-biased genes. The signal of adaptive evolution was particularly strong for X-linked male-biased genes. A comparison with data from 91 autosomal genes revealed a "fast-X" effect, in which the rate of adaptive evolution was greater for X-linked than for autosomal genes. This effect was strongest for male-biased genes but could be seen in the other groups as well. A genome-wide analysis of coding sequence divergence that accounted for sex-biased expression also uncovered a fast-X effect for male-biased and unbiased genes, suggesting that recessive beneficial mutations play an important role in adaptation.

Mitochondrial DNA Phylogeny in Eastern and Western Slavs

2008-07-15

To resolve the phylogeny of certain mitochondrial DNA (mtDNA) haplogroups in eastern Europe and estimate their evolutionary age, a total of 73 samples representing mitochondrial haplogroups U4, HV*, and R1 were selected for complete mitochondrial genome sequencing from a collection of about 2,000 control region sequences sampled in eastern (Russians, Belorussians, and Ukrainians) and western (Poles, Czechs, and Slovaks) Slavs. On the basis of whole-genome resolution, we fully characterized a number of haplogroups (HV3, HV4, U4a1, U4a2, U4a3, U4b, U4c, U4d, and R1a) that were previously described only partially. Our findings demonstrate that haplogroups HV3, HV4, and U4a1 could be traced back to the pre-Neolithic times (~12,000–19,000 years before present [YBP]) in eastern Europe. In addition, an ancient connection between the Caucasus/Europe and India has been revealed by analysis of haplogroup R1 diversity, with a split between the Indian and Caucasus/European R1a lineages occurring about 16,500 years ago. Meanwhile, some mtDNA subgroups detected in Slavs (such as U4a2a, U4a2*, HV3a, and R1a1) are definitely younger being dated between 6,400 and 8,200 YBP. However, robust age estimations appear to be problematic due to the high ratios of nonsynonymous to synonymous substitutions found in young mtDNA subclusters.

Signature Genes as a Phylogenomic Tool

Dutilh, B. E., Snel, B., Ettema, T. J. G., Huynen, M. A. — 2008-07-15

Gene content has been shown to contain a strong phylogenetic signal, yet its usage for phylogenetic questions is hampered by horizontal gene transfer and parallel gene loss and until now required completely sequenced genomes. Here, we introduce an approach that allows the phylogenetic signal in gene content to be applied to any set of sequences, using signature genes for phylogenetic classification. The hundreds of publicly available genomes allow us to identify signature genes at various taxonomic depths, and we show how the presence of signature genes in an unspecified sample can be used to characterize its taxonomic composition. We identify 8,362 signature genes specific for 112 prokaryotic taxa. We show that these signature genes can be used to address phylogenetic questions on the basis of gene content in cases where classic gene content or sequence analyses provide an ambiguous answer, such as for Nanoarchaeum equitans, and even in cases where complete genomes are not available, such as for metagenomics data. Cross-validation experiments leaving out up to 30% of the species show that ~92% of the signature genes correctly place the species in a related clade. Analyses of metagenomics data sets with the signature gene approach are in good agreement with the previously reported species distributions based on phylogenetic analysis of marker genes. Summarizing, signature genes can complement traditional sequence-based methods in addressing taxonomic questions.

Widespread Ultraconservation Divergence in Primates

Ovcharenko, I. — 2008-07-15

The distribution and evolution of ultraconserved elements (UCEs, DNA stretches that are perfectly identical in primates and rodents) were examined in genomes of 3 primate species (human, chimpanzee, and rhesus macaque). It was found that the number of UCEs has decreased throughout primate evolution. At least 26% of ancestral UCEs have diverged in hominoids, whereas an additional 17% have accumulated one or more single nucleotide polymorphisms in the human genome. Sequence polymorphism analyses indicate that mutation fixation within an UCE can trigger a relaxation in the selective constraint on that element. Homogeneous mutation accumulations in UCEs served as a template by which purifying selection acted more effectively on protein-coding UCEs. Gene ontology annotation suggests that UCE sequence variation, primarily occurring in noncoding regions, might be linked to the reprogramming of the expression pattern of transcription factors and developmentally important genes. Many of these genes are expressed in the central nervous system. Finally, UCE sequence variability within human populations has been identified, including population-specific nonsynonymous changes in protein-coding regions.

Reconstructing Evolutionary Graphs: 3D Parsimony

Lake, J. A. — 2008-07-15

The increasing recognition that symbioses have greatly altered evolution through genome fusions is creating a need for algorithms that can reliably detect and reconstruct fusions. Here, we generalize the bootstrappers gambit algorithm (a quartet method) in order to permit it to analyze both bifurcations and fusions under a single mathematical model, and thereby detect past genomic branchings and endosymbioses. This new method, 3-dimensional parsimony, can be applied to aligned sequences, such as gene, indel, or other genomic presence/absence sequences. It also provides a statistical measure of support for each possible graph. The usefulness of this method is demonstrated by applying it to the ring of life.

Spatial and Temporal Heterogeneity in Nucleotide Sequence Evolution

Whelan, S. — 2008-07-15

Models of nucleotide substitution make many simplifying assumptions about the evolutionary process, including that the same process acts on all sites in an alignment and on all branches on the phylogenetic tree. Many studies have shown that in reality the substitution process is heterogeneous and that this variability can introduce systematic errors into many forms of phylogenetic analyses. I propose a new rigorous approach for describing heterogeneity called a temporal hidden Markov model (THMM), which can distinguish between among site (spatial) heterogeneity and among lineage (temporal) heterogeneity. Several versions of the THMM are applied to 16 sets of aligned sequences to quantitatively assess the different forms of heterogeneity acting within them. The most general THMM provides the best fit in all the data sets examined, providing strong evidence of pervasive heterogeneity during evolution. Investigating individual forms of heterogeneity provides further insights. In agreement with previous studies, spatial rate heterogeneity (rates across sites [RAS]) is inferred to be the single most prevalent form of heterogeneity. Interestingly, RAS appears so dominant that failure to independently include it in the THMM masks other forms of heterogeneity, particularly temporal heterogeneity. Incorporating RAS into the THMM reveals substantial temporal and spatial heterogeneity in nucleotide composition and bias toward transition substitution in all alignments examined, although the relative importance of different forms of heterogeneity varies between data sets. Furthermore, the improvements in model fit observed by adding complexity to the model suggest that the THMMs used in this study do not capture all the evolutionary heterogeneity occurring in the data. These observations all indicate that current tests may consistently underestimate the degree of temporal heterogeneity occurring in data. Finally, there is a weak link between the amount of heterogeneity detected and the level of divergence between the sequences, suggesting that variability in the evolutionary process will be a particular problem for deep phylogeny.

Nontoxic Strains of Cyanobacteria Are the Result of Major Gene Deletion Events Induced by a Transposable Element

Christiansen, G., Molitor, C., Philmus, B., Kurmayer, R. — 2008-07-15

Blooms that are formed by cyanobacteria consist of toxic and nontoxic strains. The mechanisms that result in the occurrence of nontoxic strains are enigmatic. All the nontoxic strains of the filamentous cyanobacterium Planktothrix that were isolated from 9 European countries were found to have lost 90% of a large microcystin synthetase (mcy) gene cluster that encoded the synthesis of the toxic peptide microcystin (MC). Those strains still contain the flanking regions of the mcy gene cluster along with remnants of the transposable elements that are found in between. The majority of the strains still contain a gene coding for a distinct thioesterase type II (mcyT), which is putatively involved in MC synthesis. The insertional inactivation of mcyT in an MC-producing strain resulted in the reduction of MC synthesis by 94 ± 2% (1 standard deviation). Nontoxic strains that occur in shallow lakes throughout Europe form a monophyletic lineage. A second lineage consists of strains that contain the mcy gene cluster but differ in their photosynthetic pigment composition, which is due to the occurrence of strains that contain phycocyanin or large amounts of phycoerythrin in addition to phycocyanin. Strains containing phycoerythrin typically occur in deep-stratified lakes. The rare occurrence of gene cluster deletion, paired with the evolutionary diversification of the lineages of strains that lost or still contain the mcy gene cluster, needs to be invoked in order to explain the absence or dominance of toxic cyanobacteria in various habitats.

The Sperm Proteins from Amphioxus Mirror Its Basal Position among Chordates and Redefine the Origin of Vertebrate Protamines

Eirin-Lopez, J. M., Frehlick, L. J., Chiva, M., Saperas, N., Ausio, J. — 2008-07-15

The sperm nuclear basic proteins (SNBPs) that participate in chromatin condensation in spermatozoa belong to 3 groups: histone (H), protamine-like (PL), and protamine (P) type. They share a common origin with histone H1 resulting from the segregation of PL components, corresponding to different regions of an H1 precursor molecule (N-terminal, winged-helix, C-terminal domains), becoming independent and following a subsequent process of parallel vertical evolution (H PL P). In the present work, we describe the sequence and primary structure of the main SNBP component in the sperm of the cephalochordate Branchiostoma floridae (amphioxus), revealing that it represents the deuterostome counterpart of the PL-III SNBP component from molluscs corresponding to the H1 N-terminal region. Until now, this has been a missing piece needed to complete the evolutionary history of SNBPs in metazoan genomes. The discovery of this PL lineage in deuterostomes definitively validates the parallel vertical evolution of SNBPs across metazoans, giving further support to the "basal" position of amphioxus among chordates, with respect to tunicates. Sequence analyses suggest that later on in evolution, the appearance of positively selected arginine-rich protamines, derived from the H1 C-terminal region, led to the extinction of this PL lineage in the genomes of early protostomes and deuterostomes. Given that tunicates are now viewed as a sister group of vertebrates, the lysine to arginine transition responsible for the origin of vertebrate protamines must be set a step back from tunicates.

Inferring Natural Selection on Fine-Scale Chromatin Organization in Yeast

Babbitt, G. A., Kim, Y. — 2008-07-15

Despite its potential role in the evolution of complex phenotypes, the detection of negative (purifying) and positive selection on noncoding regulatory sequence has been elusive because of the inherent difficulty in predicting the functional consequences of mutations on noncoding sequence. Because the functioning of regulatory sequence depends upon both chromatin configuration and cis-regulatory factor binding, we investigate the idea that the functional conservation of regulatory regions should be associated with the conservation of sequence-dependent bending properties of DNA that determine its affinity for the nucleosome. Recent advances in the computational prediction of sequence-dependent affinity to nucleosomes provide an opportunity to distinguish between neutral and nonneutral evolution of fine-scale chromatin organization. Here, a statistical test is presented for detecting evolutionary conservation and/or adaptive evolution of nucleosome affinity from interspecies comparisons of DNA sequences. Local nucleosome affinities of homologous sequences were calculated using 2 recently published methods. A randomization test was applied to sites of mutation to evaluate the similarity of DNA–nucleosome affinity between several closely related species of Saccharomyces yeast. For most of the genes we analyzed, the conservation of local nucleosome affinity was detected at a few distinct locations in the upstream noncoding region. Our results also demonstrate that different patterns of chromatin evolution have shaped DNA–nucleosome interaction at the core promoters of TATA-containing and TATA-less genes and that elevated purifying selection has maintained low affinity for nucleosome in the core promoters of the latter group. Across the entire yeast genome, DNA–nucleosome interaction was also discovered to be significantly more conserved in TATA-less genes compared with TATA-containing genes.

The Impact of Founder Events on Chromosomal Variability in Multiply Mating Species

Pool, J. E., Nielsen, R. — 2008-07-15

In species with heterogametic males, the relative levels of X chromosome versus autosome diversity hold key information about the evolutionary forces at work in a population. It has been shown that population size changes alter the ratio of X linked to autosomal (X/A) variation, with population size reductions and recent bottlenecks leading to decreased X/A diversity ratios. Here we use theory and simulation to investigate a separate demographic effect—that of founder events involving multiply mated females—and find that it leads to much stronger reductions in X/A diversity ratios than are produced by simple population size changes. Investigating the potential of this process to account for sharply reduced X-linked diversity in European Drosophila melanogaster, we find that this model yields predictions that are compatible with the empirical data.

Programmed Genetic Instability: A Tumor-Permissive Mechanism for Maintaining the Evolvability of Higher Species through Methylation-Dependent Mutation of DNA Repair Genes in the Male Germ Line

Zhao, Y., Epstein, R. J. — 2008-07-15

Tumor suppressor genes are classified by their somatic behavior either as caretakers (CTs) that maintain DNA integrity or as gatekeepers (GKs) that regulate cell survival, but the germ line role of these disease-related gene subgroups may differ. To test this hypothesis, we have used genomic data mining to compare the features of human CTs (n = 38), GKs (n = 36), DNA repair genes (n = 165), apoptosis genes (n = 622), and their orthologs. This analysis reveals that repair genes are numerically less common than apoptosis genes in the genomes of multicellular organisms (P < 0.01), whereas CT orthologs are commoner than GK orthologs in unicellular organisms (P < 0.05). Gene targeting data show that CTs are less essential than GKs for survival of multicellular organisms (P < 0.0005) and that CT knockouts often permit offspring viability at the cost of male sterility. Patterns of human familial oncogenic mutations confirm that isolated CT loss is commoner than is isolated GK loss (P < 0.00001). In sexually reproducing species, CTs appear subject to less efficient purifying selection (i.e., higher Ka/Ks) than GKs (P = 0.000003); the faster evolution of CTs seems likely to be mediated by gene methylation and reduced transcription-coupled repair, based on differences in dinucleotide patterns (P = 0.001). These data suggest that germ line CT/repair gene function is relatively dispensable for survival, and imply that milder (e.g., epimutational) male prezygotic repair defects could enhance sperm variation—and hence environmental adaptation and speciation—while sparing fertility. We submit that CTs and repair genes are general targets for epigenetically initiated adaptive evolution, and propose a model in which human cancers arise in part as an evolutionarily programmed side effect of age- and damage-inducible genetic instability affecting both somatic and germ line lineages.

Gene Flow and Natural Selection in Oceanic Human Populations Inferred from Genome-Wide SNP Typing

2008-07-15

It is suggested that the major prehistoric human colonizations of Oceania occurred twice, namely, about 50,000 and 4,000 years ago. The first settlers are considered as ancestors of indigenous people in New Guinea and Australia. The second settlers are Austronesian-speaking people who dispersed by voyaging in the Pacific Ocean. In this study, we performed genome-wide single-nucleotide polymorphism (SNP) typing on an indigenous Melanesian (Papuan) population, Gidra, and a Polynesian population, Tongans, by using the Affymetrix 500K assay. The SNP data were analyzed together with the data of the HapMap samples provided by Affymetrix. In agreement with previous studies, our phylogenetic analysis indicated that indigenous Melanesians are genetically closer to Asians than to Africans and European Americans. Population structure analyses revealed that the Tongan population is genetically originated from Asians at 70% and indigenous Melanesians at 30%, which thus supports the so-called Slow train model. We also applied the SNP data to genome-wide scans for positive selection by examining haplotypic variation and identified many candidates of locally selected genes. Providing a clue to understand human adaptation to environments, our approach based on evolutionary genetics must contribute to revealing unknown gene functions as well as functional differences between alleles. Conversely, this approach can also shed some light onto the invisible phenotypic differences between populations.

Frequent, Phylogenetically Local Horizontal Transfer of the cox1 Group I Intron in Flowering Plant Mitochondria

Sanchez-Puerta, M. V., Cho, Y., Mower, J. P., Alverson, A. J., Palmer, J. D. — 2008-07-15

Horizontal gene transfer is surprisingly common among plant mitochondrial genomes. The first well-established case involves a homing group I intron in the mitochondrial cox1 gene shown to have been frequently acquired via horizontal transfer in angiosperms. Here, we report extensive additional sampling of angiosperms, including 85 newly sequenced introns from 30 families. Analysis of all available data leads us to conclude that, among the 640 angiosperms (from 212 families) whose cox1 intron status has been characterized thus far, the intron has been acquired via roughly 70 separate horizontal transfer events. We propose that the intron was originally seeded into angiosperms by a single transfer from fungi, with all subsequent inferred transfers occurring from one angiosperm to another. The pattern of angiosperm-to-angiosperm transfer is biased toward exchanges between plants belonging to the same family. Illegitimate pollination is proposed as one potential factor responsible for this pattern, given that aberrant, cross-species pollination is more likely between close relatives. Other potential factors include shared vectoring agents or common geographic locations. We report the first apparent cases of loss of the cox1 intron; losses are accompanied by retention of the exonic coconversion tract, which is located immediately downstream of the intron and which is a product of the intron's self-insertion mechanism. We discuss the many reasons why the cox1 intron is so frequently and detectably transferred, and rarely lost, and conclude that it should be regarded as the "canary in the coal mine" with respect to horizontal transfer in angiosperm mitochondria.

Compensatory Evolution in RNA Secondary Structures Increases Substitution Rate Variation among Sites

2008-07-15

There is growing evidence that interactions between biological molecules (e.g., RNA–RNA, protein–protein, RNA–protein) place limits on the rate and trajectory of molecular evolution. Here, by extending Kimura's model of compensatory evolution at interacting sites, we show that the ratio of transition to transversion substitutions () at interacting sites should be equal to the square of the ratio at independent sites. Because transition mutations generally occur at a higher rate than transversions, the model predicts that should be higher at interacting sites than at independent sites. We tested this prediction in 10 RNA secondary structures by comparing phylogenetically derived estimates of in paired sites within stems (_p) and unpaired sites within loops (_u). Eight of the 10 structures showed an excellent match to the quantitative predictions of the model, and 9 of the 10 structures matched the qualitative prediction _p > _u. Only the Rev response element from the human immunovirus (HIV) genome showed the reverse pattern, with _p < _u. Although a variety of evolutionary forces could produce quantitative deviations from the model predictions, the reversal in magnitude of _p and _u could be achieved only by violating the model assumption that the underlying transition (or transversion) mutation rates were identical in paired and unpaired regions of the molecule. We explore the ability of the APOBEC3 enzymes, host defense mechanisms against retroviruses, which induce transition mutations preferentially in single-stranded regions of the HIV genome, to explain this exception to the rule. Taken as a whole, our findings suggest that may have utility as a simple diagnostic to evaluate proposed secondary structures.

Investigation of Loss and Gain of Introns in the Compact Genomes of Pufferfishes (fugu and Tetraodon)

Loh, Y.-H., Brenner, S., Venkatesh, B. — 2008-07-15

Non-Mendelian Inheritance of Paralogs of 2 Cytoskeletal Genes in the Ciliate Chilodonella uncinata

Robinson, T., Katz, L. A. — 2008-07-15