SEARCH · PubMed Health

Results for “Purifying selection”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36Linked to original sources

PubMed2006

Fast sequence evolution of Hox and Hox-derived genes in the genus Drosophila.

BACKGROUND: It is expected that genes that are expressed early in development and have a complex expression pattern are under strong purifying selection and thus evolve slowly. Hox genes fulfill these criteria and thus, should have a low evolutionary rate. However, some observations point to a completely different scenario. Hox genes are usually highly conserved inside the homeobox, but very variable outside it. RESULTS: We have measured the rates of nucleotide divergence and indel fixation of three Hox genes, labial (lab), proboscipedia (pb) and abdominal-A (abd-A), and compared them with those of three genes derived by duplication from Hox3, bicoid (bcd), zerknüllt (zen) and zerknüllt-related (zen2), and 15 non-Hox genes in sets of orthologous sequences of three species of the genus Drosophila. These rates were compared to test the hypothesis that Hox genes evolve slowly. Our results show that the evolutionary rate of Hox genes is higher than that of non-Hox genes when both amino acid differences and indels are taken into account: 43.39% of the amino acid sequence is altered in Hox genes, versus 30.97% in non-Hox genes and 64.73% in Hox-derived genes. Microsatellites scattered along the coding sequence of Hox genes explain partially, but not fully, their fast sequence evolution. CONCLUSION: These results show that Hox genes have a higher evolutionary dynamics than other developmental genes, and emphasize the need to take into account indels in addition to nucleotide substitutions in order to accurately estimate evolutionary rates.

Amino Acid Sequence↗

PubMed2005

Investigating hookworm genomes by comparative analysis of two Ancylostoma species.

BACKGROUND: Hookworms, infecting over one billion people, are the mostly closely related major human parasites to the model nematode Caenorhabditis elegans. Applying genomics techniques to these species, we analyzed 3,840 and 3,149 genes from Ancylostoma caninum and A. ceylanicum. RESULTS: Transcripts originated from libraries representing infective L3 larva, stimulated L3, arrested L3, and adults. Most genes are represented in single stages including abundant transcripts like hsp-20 in infective L3 and vit-3 in adults. Over 80% of the genes have homologs in C. elegans, and nearly 30% of these were with observable RNA interference phenotypes. Homologies were identified to nematode-specific and clade V specific gene families. To study the evolution of hookworm genes, 574 A. caninum/A. ceylanicum orthologs were identified, all of which were found to be under purifying selection with distribution ratios of nonsynonymous to synonymous amino acid substitutions similar to that reported for C. elegans/C. briggsae orthologs. The phylogenetic distance between A. caninum and A. ceylanicum is almost identical to that for C. elegans/C. briggsae. CONCLUSION: The genes discovered should substantially accelerate research toward better understanding of the parasites' basic biology as well as new therapies including vaccines and novel anthelmintics.

Ancylostoma↗

PubMed2005

Pigs in sequence space: a 0.66X coverage pig genome survey based on shotgun sequencing.

BACKGROUND: Comparative whole genome analysis of Mammalia can benefit from the addition of more species. The pig is an obvious choice due to its economic and medical importance as well as its evolutionary position in the artiodactyls. RESULTS: We have generated approximately 3.84 million shotgun sequences (0.66X coverage) from the pig genome. The data are hereby released (NCBI Trace repository with center name "SDJVP", and project name "Sino-Danish Pig Genome Project") together with an initial evolutionary analysis. The non-repetitive fraction of the sequences was aligned to the UCSC human-mouse alignment and the resulting three-species alignments were annotated using the human genome annotation. Ultra-conserved elements and miRNAs were identified. The results show that for each of these types of orthologous data, pig is much closer to human than mouse is. Purifying selection has been more efficient in pig compared to human, but not as efficient as in mouse, and pig seems to have an isochore structure most similar to the structure in human. CONCLUSION: The addition of the pig to the set of species sequenced at low coverage adds to the understanding of selective pressures that have acted on the human genome by bisecting the evolutionary branch between human and mouse with the mouse branch being approximately 3 times as long as the human branch. Additionally, the joint alignment of the shot-gun sequences to the human-mouse alignment offers the investigator a rapid way to defining specific regions for analysis and resequencing.

Animals↗

PubMed2006

Molecular phylogeny of the antiangiogenic and neurotrophic serpin, pigment epithelium derived factor in vertebrates.

BACKGROUND: Pigment epithelium derived factor (PEDF), a member of the serpin family, regulates cell proliferation, promotes survival of neurons, and blocks growth of new blood vessels in mammals. Defining the molecular phylogeny of PEDF by bioinformatic analysis is one approach to understanding the link between its gene structure and its function in these biological processes. RESULTS: From a comprehensive search of available DNA databases we identified a single PEDF gene in all vertebrate species examined. These included four mammalian and six non-mammalian vertebrate species in which PEDF had not previously been described. A five gene cluster around PEDF was found in an approximate 100 kb region in mammals, birds, and amphibians. In ray-finned fish these genes are scattered over three chromosomes although only one PEDF gene was consistently found. The PEDF gene is absent in invertebrates including Drosophila melanogaster (D. melanogaster), Caenorhabditis elegans (C. elegans), and sea squirt (C. intestinalis). The PEDF gene is transcribed in all vertebrate phyla, suggesting it is biologically active throughout vertebrate evolution. The multiple actions of PEDF are likely conserved in evolution since it has the same gene structure across phyla, although the size of the gene ranges from 48.3 kb in X. tropicalis to 2.9 kb in fugu, with human PEDF at a size of 15.6 kb. A strong similarity in the proximal 200 bp of the PEDF promoter in mammals suggests the existence of a possible regulatory region across phyla. Using a non-synonymous/synonymous substitution rate ratio we show that mammalian and fish PEDFs have similar ratios of <0.13, reflecting a strong purifying selection of PEDF gene. A large number of repetitive transposable elements of the SINE and LINE class were found with random distribution in both the promoter and introns of mammalian PEDF. CONCLUSION: The PEDF gene first appears in vertebrates and our studies suggest that the regulation and biological actions of this gene are preserved across vertebrates. This comprehensive analysis of the PEDF gene across phyla provides new information that will aid further characterization of common functional motifs of this serpin in biological processes.

Animals↗

PubMed2006

Evolutionary anatomies of positions and types of disease-associated and neutral amino acid mutations in the human genome.

BACKGROUND: Amino acid mutations in a large number of human proteins are known to be associated with heritable genetic disease. These disease-associated mutations (DAMs) are known to occur predominantly in positions essential to the structure and function of the proteins. Here, we examine how the relative perpetuation and conservation of amino acid positions modulate the genome-wide patterns of 8,627 human disease-associated mutations (DAMs) reported in 541 genes. We compare these patterns with 5,308 non-synonymous Single Nucleotide Polymorphisms (nSNPs) in 2,592 genes from primary SNP resources. RESULTS: The abundance of DAMs shows a negative relationship with the evolutionary rate of the amino acid positions harboring them. An opposite trend describes the distribution of nSNPs. DAMs are also preferentially found in the amino acid positions that are retained (or present) in multiple vertebrate species, whereas the nSNPs are over-abundant in the positions that have been lost (or absent) in the non-human vertebrates. These observations are consistent with the effect of purifying selection on natural variation, which also explains the existence of lower minor nSNP allele frequencies at highly-conserved amino acid positions. The biochemical severity of the inter-specific amino acid changes is also modulated by natural selection, with the fast-evolving positions containing more radical amino acid differences among species. Similarly, DAMs associated with early-onset diseases are more radical than those associated with the late-onset diseases. A small fraction of DAMs (10%) overlap with the amino acid differences between species within the same position, but are biochemically the most conservative group of amino acid differences in our datasets. Overlapping DAMs are found disproportionately in fast-evolving amino acid positions, which, along with the conservative nature of the amino acid changes, may have allowed some of them to escape natural selection until compensatory changes occur. CONCLUSION: The consistency and predictability of genome-wide patterns of disease- associated and neutral amino acid variants reported here underscores the importance of the consideration of evolutionary rates of amino acid positions in clinical and population genetic analyses aimed at understanding the nature and fate of disease-associated and neutral population variation. Establishing such general patterns is an early step in efforts to diagnose the pathogenic potentials of novel amino acid mutations.

Amino Acid Substitution↗

PubMed2006

Fast rate of evolution in alternatively spliced coding regions of mammalian genes.

BACKGROUND: At least half of mammalian genes are alternatively spliced. Alternative isoforms are often genome-specific and it has been suggested that alternative splicing is one of the major mechanisms for generating protein diversity in the course of evolution. Another way of looking at alternative splicing is to consider sequence evolution of constitutive and alternative regions of protein-coding genes. Indeed, it turns out that constitutive and alternative regions evolve in different ways. RESULTS: A set of 3029 orthologous pairs of human and mouse alternatively spliced genes was considered. The rate of nonsynonymous substitutions (dN), the rate of synonymous substitutions (dS), and their ratio (omega = dN/dS) appear to be significantly higher in alternatively spliced coding regions compared to constitutive regions. When N-terminal, internal and C-terminal alternatives are analysed separately, C-terminal alternatives appear to make the main contribution to the observed difference. The effects become even more pronounced in a subset of fast evolving genes. CONCLUSION: These results provide evidence of weaker purifying selection and/or stronger positive selection in alternative regions and thus one more confirmation of accelerated evolution in alternative regions. This study corroborates the theory that alternative splicing serves as a testing ground for molecular evolution.

Alternative Splicing↗

PubMed2006

Phylogeographic support for horizontal gene transfer involving sympatric bruchid species.

BACKGROUND: We report on the probable horizontal transfer of a mitochondrial gene, cytb, between species of Neotropical bruchid beetles, in a zone where these species are sympatric. The bruchid beetles Acanthoscelides obtectus, A. obvelatus, A. argillaceus and Zabrotes subfasciatus develop on various bean species in Mexico. Whereas A. obtectus and A. obvelatus develop on Phaseolus vulgaris in the Mexican Altiplano, A. argillaceus feeds on P. lunatus in the Pacific coast. The generalist Z. subfasciatus feeds on both bean species, and is sympatric with A. obtectus and A. obvelatus in the Mexican Altiplano, and with A. argillaceus in the Pacific coast. In order to assess the phylogenetic position of these four species, we amplified and sequenced one nuclear (28S rRNA) and two mitochondrial (cytb, COI) genes. RESULTS: Whereas species were well segregated in topologies obtained for COI and 28S rRNA, an unexpected pattern was obtained in the cytb phylogenetic tree. In this tree, individuals from A. obtectus and A. obvelatus, as well as Z. subfasciatus individuals from the Mexican Altiplano, clustered together in a unique little variable monophyletic unit. In contrast, A. argillaceus and Z. subfasciatus individuals from the Pacific coast clustered in two separated clades, identically to the pattern obtained for COI and 28S rRNA. An additional analysis showed that Z. subfasciatus individuals from the Mexican Altiplano also possessed the cytb gene present in individuals of this species from the Pacific coast. Zabrotes subfasciatus individuals from the Mexican Altiplano thus demonstrated two cytb genes, an "original" one and an "infectious" one, showing 25% of nucleotide divergence. The "infectious" cytb gene seems to be under purifying selection and to be expressed in mitochondria. CONCLUSION: The high degree of incongruence of the cytb tree with patterns for other genes is discussed in the light of three hypotheses: experimental contamination, hybridization, and pseudogenisation. However, none of these seem able to explain the patterns observed. A fourth hypothesis, involving recent horizontal gene transfer (HGT) between A. obtectus and A. obvelatus, and from one of these species to Z. subfasciatus in the Mexican Altiplano, seems the only plausible explanation. The HGT between our study species seems to have occurred recently, and only in a zone where the three beetles are sympatric and share common host plants. This suggests that transfer could have been effected by some external vector such as a eukaryotic or viral parasite, which might still host the transferred fragment. REVIEWERS: This article was reviewed by Eric Bapteste, Adam Eyre-Walker and Alexey Kondrashov.

Journal Article↗

PubMed2026

Genome-wide characterization of the bZIP gene family in Rattus norvegicus and expression profiling analysis during brain development.

BACKGROUND: The brown rat (Rattus norvegicus) serves as a cornerstone model organism in biomedical research, particularly for understanding physiological homeostasis and stress responses. The basic leucine zipper (bZIP) transcription factor family is a pivotal regulatory network involved in growth, organogenesis, and neurodevelopment. Despite its importance, a systematic characterization of the bZIP gene family in rats has remained elusive. RESULTS: In this study, we performed a genome-wide identification of 61 RnbZIP genes, which were categorized into 10 distinct subfamilies based on phylogenetic relationships and chromosomal localization. Structural analysis revealed conserved motif arrangements within subfamilies, while collinearity analysis identified significant gene duplication events-predominantly tandem and segmental duplications-that have driven the evolutionary expansion of the RnbZIP family. Quantitative analysis showed that members within the same subfamily shared 45%-92% sequence similarity (calculated using the BLOSUM62 scoring matrix), and all duplicated gene pairs underwent strong purifying selection (Ka/Ks < 1). Comparative genomics across seven rodent species further underscored the evolutionary conservation and divergence of these factors. Expression profiling across diverse organs and brain developmental stages indicated that RnbZIP genes exhibit high tissue specificity. Notably, 10 candidate genes, including RnbZIP01, RnbZIP02, and RnbZIP08, demonstrated dynamic expression patterns during brain maturation, suggesting their essential roles in neurodevelopmental processes. CONCLUSIONS: Our findings provide a comprehensive structural and evolutionary framework for the RnbZIP gene family, highlighting their potential regulatory functions in rat organogenesis and brain development. This study establishes a valuable resource for further functional characterization of specific bZIP members in mammalian neurological systems.

Animals↗

PubMed2026

Identification of the R2R3-MYB gene family in wild jujube (Ziziphus jujuba var. spinosa) and analysis of its expression under drought stress.

BACKGROUND: R2R3-MYB gene family serves as a pivotal regulatory factor in plant growth, development, and responses to environmental stresses. To investigate its function in the drought stress response of wild jujube (Ziziphus jujuba Mill. var. spinosa), a typical eco-economic forest species, this study performed genome-wide identification and relevant analyses of R2R3-MYB genes. RESULTS: A total of 91 R2R3-MYB genes (designated as ZjMYB1 to ZjMYB91) were identified, which were unevenly distributed across 12 chromosomes. These genes mainly encode hydrophilic and unstable proteins, 97.8% of which are localized in the nucleus. Phylogenetic analysis classified these genes into 25 clades, showing evolutionary conservation and species-specific divergence with the R2R3-MYB protein family. The expansion of the ZjMYB family is mainly characterized by segmental duplication, and all duplicated gene pairs have undergone purifying selection. ZjMYBs are widely involved in plant growth and development as well as abiotic stress responses, with the highest expression level particularly in leaf tissues; a total of 13 genes were specifically annotated as water deficit response-related genes in drought stress and abscisic acid (ABA) signaling pathways. Integrating the above analyses together with transcriptome data and qRT-PCR validation results revealed that ZjMYB5, ZjMYB53, ZjMYB57 and ZjMYB85 function as core drought-responsive genes, which display both tissue-specific and time-dependent expression patterns under drought stress. CONCLUSIONS: This study systematically elucidated the functional characteristics and regulatory network of the R2R3-MYB gene family in wild jujube, providing critical genetic resources and a theoretical basis for dissecting the molecular mechanisms underlying drought tolerance in wild jujube and breeding drought-resistant cultivars.

Ziziphus↗

PubMed2005

Pleiotropic functions of a conserved insect-specific Hox peptide motif.

The proteins that regulate developmental processes in animals have generally been well conserved during evolution. A few cases are known where protein activities have functionally evolved. These rare examples raise the issue of how highly conserved regulatory proteins with many roles evolve new functions while maintaining old functions. We have investigated this by analyzing the function of the ;QA' peptide motif of the Hox protein Ultrabithorax (Ubx), a motif that has been conserved throughout insect evolution since its establishment early in the lineage. We precisely deleted the QA motif at the endogenous locus via allelic replacement in Drosophila melanogaster. Although the QA motif was originally characterized as involved in the repression of limb formation, we have found that it is highly pleiotropic. Curiously, deleting the QA motif had strong effects in some tissues while barely affecting others, suggesting that QA function is preferentially required for a subset of Ubx target genes. QA deletion homozygotes had a normal complement of limbs, but, at reduced doses of Ubx and the abdominal-A (abd-A) Hox gene, ectopic limb primordia and adult abdominal limbs formed when the QA motif was absent. These results show that redundancy and the additive contributions of activity-regulating peptide motifs play important roles in moderating the phenotypic consequences of Hox protein evolution, and that pleiotropic peptide motifs that contribute quantitatively to several functions are subject to intense purifying selection.

Amino Acid Motifs↗

PubMed2000

Biochemical characterization of 60S acidic ribosomal P proteins from porcine liver and the inhibition of their immunocomplex formation with sera from systemic lupus erythematosus (SLE) patients by glycyrrhizin in vitro.

The three casein kinase II (CK-II) phosphate acceptors (p35, p17 and p15) in the Superdex CK-II fraction prepared from a 1.5 M NaCl extract of porcine liver were selectively purified by glycyrrhizin (GL)-affinity column chromatography (HPLC) as a heterocomplex associated with CK-II. Determination of the N-terminal amino acid sequences and immunological tests confirmed that these three CK-II phosphate acceptors belong to the family of 60S acidic ribosomal proteins (P0, P1 and P2). Three polyphenol-containing anti-oxidant compounds [catechin, epigallocatechin gallate (EGCG) and quercetin] inhibited CK-II activity (phosphorylation of these ribosomal P proteins) in a dose-dependent manner in vitro. Quercetin (ID50 = approx. 50 nM) was found to be an effective CK-II inhibitor. In contrast, CK-II activity was significantly stimulated by lower doses (0.3-3 microl) of GL, but was inhibited at high doses above 30 microM. As expected, GL at high doses above 200 microM inhibited the immunocomplex formation of 60S acidic ribosomal P proteins with their specific antibodies in the sera from patients with systemic lupus erythematosus (SLE). These results suggest that (i) a GL-affinity column is useful for effective purification of 60S acidic ribosomal P proteins from various mammalian cells as a heterocomplex associated with CK-II; and (ii) a relative high dose of GL may prevent the immunocomplex formation of 60S acidic ribosomal P proteins with their specific antibodies in the sera of SLE patients.

Amino Acid Sequence↗

PubMed2001

Characterization of secretory type IIA phospholipase A2 (sPLA2-IIA) as a glycyrrhizin (GL)-binding protein and the GL-induced inhibition of the CK-II-mediated stimulation of sPLA2-IIA activity in vitro.

By means of heparin-affinity and glycyrrhizin (GL)-affinity column chromatographies (HPLC), a GL-binding phospholipase A2 (gbPLA2) was selectively purified from the synovial fluids of patients with rheumatoid arthritis. This purified gbPLA2 was identified as a secretory type IIA PLA2 (sPLA2-IIA) since it was crossreacted with anti-sPLA2-IIA serum. The activity of purified sPLA2-IIA was inhibited by glycyrrhetinic acid (GA) and a GA derivative (oGA) in a dose-dependent manner, but it was more sensitive to GA than GL. Furthermore, it was found that (i) purified sPLA2-IIA is phosphorylated by casein kinase II (CK-II) in vitro; (ii) this phosphorylation induces in a significant stimulation of PLA2 activity; and (iii) oGA at one-tenth the concentration of GL inhibits the CK-II-mediated stimulation of sPLA2-IIA activity. These results show that (i) sPLA2-IIA is a GL-binding protein; and (ii) CK-II mediates stimulation of its PLA2 activity in vitro.

Anti-Inflammatory Agents, Non-Steroidal↗

PubMed2001

Molecular cloning of horse Hsp90 cDNA and its comparative analysis with other vertebrate Hsp90 sequences.

Heat shock protein 90 (Hsp90), a molecular chaperone, is ubiquitous and involved in numerous cellular processes. To contribute to the relatively small collection of vertebrate Hsp90 sequences in the gene data bank, we cloned and sequenced horse (Equus caballus) Hsp90 alpha and beta cDNAs. This enabled identification of horse-specific primers for development of a convenient PCR-based method that could monitor horse stress tolerance. We analyzed the sequence data comparatively and phylogenetically with other Hsp90 cDNA sequences, and identified vertebrate-specific and isoform-specific conserved regions to facilitate future molecular investigations of Hsp90 functions. We found 4 highly conserved regions to vertebrate Hsp90 exclusively and 27 amino acids conserved among but differing between Hsp90 alpha and Hsp90 beta sequences. Protein-based phylogenetic trees revealed high conservation between mammal species within Hsp90 alpha and beta clusters. Comparison of nucleotide and amino acid substitution levels suggests that horse Hsp90 beta has undergone strong purifying selection, while rat Hsp90 beta and hamster Hsp90 alpha have been positively selected. Surprisingly, fish Hsp90 alpha genes clearly clustered with Hsp90 beta genes, and no distinct placement of fish Hsp90 alpha protein was found. The Hsp90 alpha isoform is apparently the result of beta gene duplication. Our results highlight the importance of organism- and isoform-specific Hsp90 functional analyses in describing the role of Hsp90 in cells.

Amino Acid Sequence↗

PubMed2005

Low variability of the POLG (CAG)n repeat in north Eurasian populations.

We investigated the frequency of different repeat-length alleles of the trinucleotide CAG microsatellite repeat in the coding sequence of the nuclear gene for the catalytic subunit of mitochondrial DNA polymerase gamma (POLG) in 12 ethnic groups from northern Eurasia. The population sample consisted of 1,330 individuals from 3 large geographic areas: Europe, Southwest Asia, and Siberia/East Asia. We found that the 10-repeat allele of the POLG gene is the most frequent in all analyzed populations, with a frequency of 88-96%. The heterozygosity level ranges from 22% in Europe to 13.6% in Southwest Asia with the lowest value of 7.4% in Siberia/East Asia. The present study provides evidence of clinal distribution of POLG gene heterozygosity in North Eurasian populations. In general, we found an extremely low variability of the trinucleotide CAG microsatellite repeat, suggesting that purifying selection acts against deleterious alleles, although low mutability of the repeated region cannot be ruled out.

Alleles↗

PubMed2005

Emergence of young human genes after a burst of retroposition in primates.

The origin of new genes through gene duplication is fundamental to the evolution of lineage- or species-specific phenotypic traits. In this report, we estimate the number of functional retrogenes on the lineage leading to humans generated by the high rate of retroposition (retroduplication) in primates. Extensive comparative sequencing and expression studies coupled with evolutionary analyses and simulations suggest that a significant proportion of recent retrocopies represent bona fide human genes. We estimate that at least one new retrogene per million years emerged on the human lineage during the past approximately 63 million years of primate evolution. Detailed analysis of a subset of the data shows that the majority of retrogenes are specifically expressed in testis, whereas their parental genes show broad expression patterns. Consistently, most retrogenes evolved functional roles in spermatogenesis. Proteins encoded by X chromosome-derived retrogenes were strongly preserved by purifying selection following the duplication event, supporting the view that they may act as functional autosomal substitutes during X-inactivation of late spermatogenesis genes. Also, some retrogenes acquired a new or more adapted function driven by positive selection. We conclude that retroduplication significantly contributed to the formation of recent human genes and that most new retrogenes were progressively recruited during primate evolution by natural and/or sexual selection to enhance male germline function.

Animals↗

PubMed2006

Assessing the accuracy of ancestral protein reconstruction methods.

The phylogenetic inference of ancestral protein sequences is a powerful technique for the study of molecular evolution, but any conclusions drawn from such studies are only as good as the accuracy of the reconstruction method. Every inference method leads to errors in the ancestral protein sequence, resulting in potentially misleading estimates of the ancestral protein's properties. To assess the accuracy of ancestral protein reconstruction methods, we performed computational population evolution simulations featuring near-neutral evolution under purifying selection, speciation, and divergence using an off-lattice protein model where fitness depends on the ability to be stable in a specified target structure. We were thus able to compare the thermodynamic properties of the true ancestral sequences with the properties of "ancestral sequences" inferred by maximum parsimony, maximum likelihood, and Bayesian methods. Surprisingly, we found that methods such as maximum parsimony and maximum likelihood that reconstruct a "best guess" amino acid at each position overestimate thermostability, while a Bayesian method that sometimes chooses less-probable residues from the posterior probability distribution does not. Maximum likelihood and maximum parsimony apparently tend to eliminate variants at a position that are slightly detrimental to structural stability simply because such detrimental variants are less frequent. Other properties of ancestral proteins might be similarly overestimated. This suggests that ancestral reconstruction studies require greater care to come to credible conclusions regarding functional evolution. Inferred functional patterns that mimic reconstruction bias should be reevaluated.

Algorithms↗

PubMed2006

The evolution of word composition in metazoan promoter sequence.

The field of molecular evolution provides many examples of the principle that molecular differences between species contain information about evolutionary history. One surprising case can be found in the frequency of short words in DNA: more closely related species have more similar word compositions. Interest in this has often focused on its utility in deducing phylogenetic relationships. However, it is also of interest because of the opportunity it provides for studying the evolution of genome function. Word-frequency differences between species change too slowly to be purely the result of random mutational drift. Rather, their slow pattern of change reflects the direct or indirect action of purifying selection and the presence of functional constraints. Many such constraints are likely to exist, and an important challenge is to distinguish them. Here we develop a method to do so by isolating the effects acting at different word sizes. We apply our method to 2-, 4-, and 8-base-pair (bp) words across several classes of noncoding sequence. Our major result is that similarities in 8-bp word frequencies scale with evolutionary time for regions immediately upstream of genes. This association is present although weaker in intronic sequence, but cannot be detected in intergenic sequence using our method. In contrast, 2-bp and 4-bp word frequencies scale with time in all classes of noncoding sequence. These results suggest that different genomic processes are involved at different word sizes. The pattern in 2-bp and 4-bp words may be due to evolutionary changes in processes such as DNA replication and repair, as has been suggested before. The pattern in 8-bp words may reflect evolutionary changes in gene-regulatory machinery, such as changes in the frequencies of transcription-factor binding sites, or in the affinity of transcription factors for particular sequences.

Amino Acids↗

PubMed2026

An ancient alpharetrovirus lineage in bats: Evolutionary insights and possible roles in reproduction.

Alpharetroviruses are an important group of pathogens known to cause leukemias and tumors, and were historically considered to be restricted to avian hosts. The identification of alpharetrovirus-like envelopes in bat genomes has hinted at a potentially wider host range, although their relations to modern alpharetroviruses and distribution remains unclear. Through a paleovirological screening of 818 vertebrate genomes we identified CHIRalphaEnv, a lineage that belongs firmly within alpharetroviruses, and emerged from a cross-class transmission from saurian hosts. We determine that CHIRalphaEnv envelope genes have been co-opted across bats on eight separate occasions between 43.8 and 18.9 million years ago and are preserved in most bat genomes screened. CHIRalphaEnv elements encode full-length envelope proteins and have been maintained under purifying selection, demonstrating multiple instances of exaptation by their bat hosts and a likely ubiquitous function. We observe high expression levels of CHIRalphaEnv envelopes in endometrium tissue from Carollia perspicillata, suggesting an involvement in reproductive function. We find CHIRalphaEnv sequence relatives in multiple mammalian clades (Afrotherians, rodents and bats), expanding the host range and extending origins of alpharetroviruses beyond 43 million years. We also propose the first mammalian co-opted Endogenous retrovirus (ERV) derived from an Alpharetrovirus envelope and explore the convergent functional recruitment of CHIRalphaEnv in hemochorial placentation in bats, elephant shrews and spiny mice. These findings highlight alpharetroviruses as a previously underappreciated source of functional exaptation in mammals.

Animals↗