Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Evolutionary analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

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

At least 541 records · Page 30Linked to original sources

Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila.

Evolutionary analysis of codon bias in Drosophila indicates that synonymous mutations are not neutral, but rather are subject to weak selection at the translation level. Here we show that the effectiveness of natural selection on synonymous sites is strongly correlated with the rate of recombination, in accord with the nearly neutral hypothesis. This correlation, however, is apparent only in genes encoding short proteins. Long coding regions have both a lower codon bias and higher synonymous substitution rates, suggesting that they are affected less efficiently by selection. Therefore, both the length of the coding region and the recombination rate modulate codon bias. In addition, the data indicate that selection coefficients for synonymous mutations must vary by a minimum of one or two orders of magnitude. Two hypotheses are proposed to explain the relationship among the coding region length, the codon bias, and the synonymous divergence and polymorphism levels across the range of recombination rates in Drosophila. The first hypothesis is that selection coefficients on synonymous mutations are inversely related to the total length of the coding region. The second hypothesis proposes that interference among synonymous mutations reduces the efficacy of selection on these mutations. We investigated this second hypothesis by carrying out forward simulations of weakly selected mutations in model populations. These simulations show that even with realistic recombination rates, this interference, which we call the "small-scale" Hill-Robertson effect, can have a moderately strong influence on codon bias.

Animals↗

The correlation between intron length and recombination in drosophila. Dynamic equilibrium between mutational and selective forces.

Intron length is negatively correlated with recombination in both Drosophila melanogaster and humans. This correlation is not likely to be the result of mutational processes alone: evolutionary analysis of intron length polymorphism in D. melanogaster reveals equivalent ratios of deletion to insertion in regions of high and low recombination. The polymorphism data do reveal, however, an excess of deletions relative to insertions (i.e., a deletion bias), with an overall deletion-to-insertion events ratio of 1.35. We propose two types of selection favoring longer intron lengths. First, the natural mutational bias toward deletion must be opposed by strong selection in very short introns to maintain the minimum intron length needed for the intron splicing reaction. Second, selection will favor insertions in introns that increase recombination between mutations under the influence of selection in adjacent exons. Mutations that increase recombination, even slightly, will be selectively favored because they reduce interference among selected mutations. Interference selection acting on intron length mutations must be very weak, as indicated by frequency spectrum analysis of Drosophila intron length polymorphism, making the equilibrium for intron length sensitive to changes in the recombinational environment and population size. One consequence of this sensitivity is that the advantage of longer introns is expected to decrease inversely with the rate of recombination, thus leading to a negative correlation between intron length and recombination rate. Also in accord with this model, intron length differs between closely related Drosophila species, with the longest variant present more often in D. melanogaster than in D. simulans. We suggest that the study of the proposed dynamic model, taking into account interference among selected sites, might shed light on many aspects of the comparative biology of genome sizes including the C value paradox.

Animals↗

Evolution of the human ASPM gene, a major determinant of brain size.

The size of human brain tripled over a period of approximately 2 million years (MY) that ended 0.2-0.4 MY ago. This evolutionary expansion is believed to be important to the emergence of human language and other high-order cognitive functions, yet its genetic basis remains unknown. An evolutionary analysis of genes controlling brain development may shed light on it. ASPM (abnormal spindle-like microcephaly associated) is one of such genes, as nonsense mutations lead to primary microcephaly, a human disease characterized by a 70% reduction in brain size. Here I provide evidence suggesting that human ASPM went through an episode of accelerated sequence evolution by positive Darwinian selection after the split of humans and chimpanzees but before the separation of modern non-Africans from Africans. Because positive selection acts on a gene only when the gene function is altered and the organismal fitness is increased, my results suggest that adaptive functional modifications occurred in human ASPM and that it may be a major genetic component underlying the evolution of the human brain.

Anthropometry↗

A mutation in GJB3 is associated with recessive erythrokeratodermia variabilis (EKV) and leads to defective trafficking of the connexin 31 protein.

Erythrokeratodermia variabilis (EKV) is a skin disorder characterized by variable (transient) erythemas and fixed keratosis. The disorder maps to chromosome 1p34-35, a location that contains the GJB3 gene encoding the gap junction protein connexin 31. Until now, only heterozygote mutations in the form of dominant inheritance have been described in this gene associated with EKV. We report here a homozygote mutation in the connexin 31 gene, found in a family that shows recessive inheritance of the disorder, thus providing the first molecular support for a recessive variant of EKV. The entire GJB3 coding sequence was scanned for mutations by sequencing. We detected a T-->C transition at position 101 of the coding sequence, which replaces a leucine with a proline at residue 34 of the protein (L34P). Evolutionary analysis shows that this mutation is located at a highly conserved region of connexin in the first putative transmembrane helix (TMH). In transfected keratinocytes, L34P connexin 31 had a cytoplasmic distribution, suggesting that the mutant form of this protein will not form normal gap junctions between adjacent cells. The change of leucine to proline is likely to alter the structure of the first TMH of connexin by inducing a kink, thus influencing connexon structure and function.

Amino Acid Sequence↗

The IL-9 receptor gene, located in the Xq/Yq pseudoautosomal region, has an autosomal origin, escapes X inactivation and is expressed from the Y.

All human X-linked genes known so far, except for the Xp/Yp pseudoautosomal genes, are conserved as a single linkage group on the murine X chromosome. We show that the interleukin-9 (IL-9) receptor gene (IL9R), which is located within the human Xq/Yq homology region, maps to the murine chromosome 11. The Xq/Yq pseudoautosomal region (Xq PAR) thus represents a second region on the human X chromosome which is not X linked in mice. Furthermore, we show that IL9R is absent on the Y of great apes. IL9R is thus exceptional among X/Y genes in that it is X linked in some mammals, but autosomal or pseudoautosomal in others. Genes located on the X and the Y generally escape X inactivation. An exception to this rule is SYBL1, a gene located in Xq PAR. SYBL1 is X inactivated and is inactive on the Y chromosome. In contrast, we show that IL9R expression does occur from the Y, the active and the inactive X chromosomes. This finding raises the question of how the transcriptional regulation of genes within Xq PAR occurs and how the X inactivation status of IL9R has evolved following the autosome to X and the X to X/Y translocation. The evolutionary analysis of the IL9R gene, which is located at 10 kb from the telomere, and its pseudogenes at several telomeres, also provides insight into the evolution of these loci and of subtelomeric regions in general.

Animals↗

Human BOULE gene rescues meiotic defects in infertile flies.

Defects in human germ cell development are common and yet little is known of genes required for germ cell development in men and women. The pathways that develop germ cells appear to be conserved broadly, at least in outline, in organisms as diverse as flies and humans beginning with allocation of cells to the germ cell lineage, migration of these cells to the fetal gonad, mitotic proliferation and meiosis of the germ cells, and maturation into sperm and eggs. In model organisms, a few thousand genes may be required for germ cell development including meiosis. To date, however, no genes that regulate critical steps of reproduction have been shown to be functionally conserved from flies to humans. This may be due in part to strong selective pressures that are thought to drive reproductive genes to high degrees of divergence. Here, we investigated the micro- and macro-evolution of the BOULE gene, a member of the human DAZ (deleted in azoospermia) gene family, within primates, within mammals and within metazoans. We report that sequence divergence of BOULE is unexpectedly low and that rapid evolution is not detectable. We extend the evolutionary analysis of BOULE to the level of phyla and show that a human BOULE transgene can advance meiosis in infertile boule mutant flies. This is the first demonstration that a human reproductive gene can rescue reproductive defects in a fly. These studies lend strong support to the idea that BOULE may encode a key conserved switch that regulates progression of germ cells through meiosis in men.

Animals↗

GB virus C/hepatitis G virus infection in southern China.

The prevalence of GB virus C/hepatitis G virus (GBV-C/HGV) infection among intravenous drug users (IVDUs), patients with liver diseases, and blood donors in Nanning, southern China was studied. GBV-C/HGV RNA was detected by reverse transcription polymerase chain reaction with primers derived from the 5'-untranslated region. GBV-C/HGV RNA was detected in 64 of 85 IVDUs, 20 of 80 persons with liver disease, and 1 of 50 blood donors. Among IVDUs, GBV-C/HGV infection was associated with antibodies to hepatis C virus (HCV) and with hepatitis B surface antigen (HBsAg). Eleven nucleotide sequences were determined and analyzed by molecular evolutionary analysis. In a phylogenetic tree, the isolates were grouped in three clusters with GBV-C and HGV grouped in two clusters. These data indicate that GBV-C/HGV infection is common in China among IVDUs but uncommon among persons with liver disease without HBsAg or anti-HCV and that there is a new group of GBV-C/HGV.

Blood Donors↗

Signal sequence conservation and mature peptide divergence within subgroups of the murine beta-defensin gene family.

Beta-defensins are two exon genes which encode broad spectrum antimicrobial cationic peptides. We have analyzed the largest murine cluster of these genes which localizes to chromosome 8. Using hidden Markov models, we identified six beta-defensin exon 2-like sequences and subsequently found full-length expressed transcripts for these novel genes. Expression was high in brain and reproductive tissues. Eleven beta-defensins could be grouped into two clear subgroups by virtue of their position and high signal sequence (exon 1 encoded) identity. In contrast, however, there was a very low level of sequence conservation in the exon 2 region encoding the mature antimicrobial peptide. Examination of the gene sequences of orthologs in other rodents also revealed an excess of nucleotide changes that altered amino acids in the mature peptide region. Evolutionary analysis revealed strong evidence that following gene duplication, exon 1 and surrounding noncoding DNA show little divergence within subgroups. The focus for rapid sequence divergence is localized in the DNA encoding the mature peptide and this is driven by accelerated positive selection. This mechanism of evolution is consistent with the role of this gene family as defense against bacterial pathogens and the sequence changes have implications for novel antibiotic design.

Animals↗

The mitochondrial genome of Phoronis architecta--comparisons demonstrate that phoronids are lophotrochozoan protostomes.

The proper reconstruction of the relationships among the animal phyla is central to interpreting patterns of animal evolution from the genomic level to the morphological level. This is true not only of the more speciose phyla but also of smaller groups. We report here the nearly complete DNA sequence of the mitochondrial genome of the phoronid Phoronis architecta, which has a gene arrangement remarkably similar to that of a protostome animal, the chiton Katharina tunicata. Evolutionary analysis of both gene arrangements and inferred amino acid sequences of these taxa, along with those of three brachiopods and other diverse animals, strongly supports the hypothesis that lophophorates are part of the large group that includes mollusks and annelids-i.e., the Lophotrochozoa-and solidly refutes the alternative of their being deuterostomes.

Animals↗

Nonneutral evolution of the transcribed pseudogene Makorin1-p1 in mice.

Pseudogenes are nonfunctional relics of formerly functional genes and are thought to evolve neutrally. In some pseudogenes, however, the molecular evolutionary patterns are atypical of neutrally evolving sequences, exhibiting sequence conservation, codon-usage bias, and other features associated with functional genes. Makorin1-p1 is a transcribed pseudogene first identified in the mouse Mus musculus. The transcript of Makorin1-p1 can regulate the stability of the transcript of its paralogous functional gene Makorin1. Specifically, the half-life of Makorin1 mRNA increases significantly in the presence of Makorin1-p1 transcript, and targeted deletion of Makorin1-p1 is lethal in mice. Here, we show that Makorin1-p1 originated after the separation of Mus and Rattus but before the divergence of M. musculus and M. pahari. The transcribed region of Makorin1-p1 exhibits rates of point and indel substitutions that are two to four times lower than those in the untranscribed region, suggesting that the transcribed region is under functional constraint and is not evolving neutrally. Although the transcript of Makorin1-p1 likely functions by its sequence similarity to Makorin1, we find no evidence of gene conversion between them, indicating that functional conservation alone is sufficient to maintain their coordinated evolution. A duplication-degeneration model is proposed to explain how Makorin1-p1 was co-opted into the regulatory system of Makorin1. There are over 10,000 pseudogenes in a typical mammalian genome, and it is plausible that many functional but untranslatable pseudogenes exist. Our results illustrate the potential of using evolutionary analysis to identify such pseudogenes from genome sequences.

Animals↗

Positive selection in the carbohydrate recognition domains of sea urchin sperm receptor for egg jelly (suREJ) proteins.

A wealth of evidence shows that protein-carbohydrate recognition mediates the steps of gamete interaction during fertilization. Carbohydrate-recognition domains (CRDs) comprise a large family of ancient protein modules of approximately 120 amino acids, having the same protein fold, that bind terminal sugar residues on glycoproteins and polysaccharides. Sea urchin sperm express three suREJ (sea urchin receptor for egg jelly) proteins on their plasma membranes. suREJ1 has two CRDs, whereas suREJ2 and suREJ3 both have one CRD. suREJ1 binds the fucose sulfate polymer (FSP) of egg jelly to induce the sperm acrosome reaction. The structure of FSP is species specific. Therefore, the suREJ1 CRDs could encode molecular recognition between sperm and egg underlying the species-specific induction of the acrosome reaction. The functions of suREJ2 and suREJ3 have not been explored, but suREJ3 is exclusively localized on the plasma membrane over the sperm acrosomal vesicle and is physically associated with sea urchin polycystin-2, a known cation channel. An evolutionary analysis of these four CRDs was performed for six sea urchin species. Phylogenetic analysis shows that these CRDs were already differentiated in the common ancestor of these six sea urchins. The CRD phylogeny agrees with previous work on these species based on one nuclear gene and several mitochondrial genes. Maximum likelihood shows that positive selection acts on these four CRDs. Threading the suREJ CRDs onto the prototypic CRD crystal structure shows that many of the sites under positive selection are on extended loops, which are involved in saccharide binding. This is the first demonstration of positive selection in CRDs and is another example of positive selection acting on the evolution of gamete-recognition proteins.

Animals↗

Adaptive evolution of fertilization proteins within a genus: variation in ZP2 and ZP3 in deer mice (Peromyscus).

Rapid evolution of reproductive proteins has been documented in a wide variety of taxa. In internally fertilized species, knowledge about the evolutionary dynamics of these proteins between closely related taxa is primarily limited to accessory gland proteins in the semen of Drosophila. Investigation of additional taxa and functional classes of proteins is necessary in order to determine if there is a general pattern of adaptive evolution of reproductive proteins between recently diverged species. We performed an evolutionary analysis of 2 egg coat proteins, ZP2 and ZP3, in 15 species of deer mice (genus Peromyscus). Both of these proteins are involved in egg-sperm binding, a critical step in maintaining species-specific fertilization. Here, we show that Zp2 and Zp3 gene trees are not consistent with trees based on nonreproductive genes, Mc1r and Lcat, where species formed monophyletic clades. In fact, for both of the reproductive genes, intraspecific amino acid variation was extensive and alleles were sometimes shared across species. We document positive selection acting on ZP2 and ZP3 and identify specific amino acid sites that are likely targets of selection using both maximum likelihood approaches and patterns of parallel amino acid change. In ZP3, positively selected sites are clustered in and around the region implicated in sperm binding in Mus, suggesting changes may impact egg-sperm binding and fertilization potential. Finally, we identify lineages with significantly elevated rates of amino acid substitution using a Bayesian mapping approach. These findings demonstrate that the pattern of adaptive reproductive protein evolution found at higher taxonomic levels can be documented between closely related mammalian species, where reproductive isolation has evolved recently.

Adaptation, Biological↗

Mitochondrial genomes and avian phylogeny: complex characters and resolvability without explosive radiations.

We improve the taxon sampling for avian phylogeny by analyzing 7 new mitochondrial genomes (a toucan, woodpecker, osprey, forest falcon, American kestrel, heron, and a pelican). This improves inference of the avian tree, and it supports 3 major conclusions. The first is that some birds (including a parrot, a toucan, and an osprey) exhibit a complete duplication of the control region (CR) meaning that there are at least 4 distinct gene orders within birds. However, it appears that there are regions of continued gene conversion between the duplicate CRs, resulting in duplications that can be stable for long evolutionary periods. Because of this stable duplicated state, gene order can eventually either revert to the original order or change to the new gene order. The existence of this stable duplicate state explains how an apparently unlikely event (finding the same novel gene order) can arise multiple times. Although rare genomic changes have theoretical advantages for tree reconstruction, they can be compromised if these apparently rare events have a stable intermediate state. Secondly, the toucan and woodpecker improve the resolution of the 6-way split within Neoaves that has been called an "explosive radiation." An explosive radiation implies that normal microevolutionary events are insufficient to explain the observed macroevolution. By showing the avian tree is, in principle, resolvable, we demonstrate that the radiation of birds is amenable to standard evolutionary analysis. Thirdly, and as expected from theory, additional taxa breaking up long branches stabilize the position of some problematic taxa (like the falcon). In addition, we report that within the birds of prey and allies, we did not find evidence pairing New World vultures with storks or accipitrids (hawks, eagles, and osprey) with Falconids.

Animals↗

Nucleotide sequences of two serine tRNAs with a GGA anticodon: the structure-function relationships in the serine family of E. coli tRNAs.

We have determined the nucleotide sequence of the major species of E. coli tRNASer and of a minor species having the same GGA anticodon. These two tRNAs should recognize the UCC and UCU codons, the most widely used codons for serine in the highly expressed genes of E. coli. The two sequences differ in only one position of the D-loop. Neither tRNA has a modified adenosine in the position 3'-adjacent to the anticodon. This can be rationalized on the basis of a structural constraint in the anticodon stem and may be related to optimization of the codon-anticodon interaction. Comparison of all E.coli serine tRNAs (and that encoded by bacteriophage T4) reveals characteristic (possibly functional) features. Evolutionary analysis suggests an eubacterial origin of the T4 tRNASer gene and the existence of a recent common ancestor for the tRNASerGGA and tRNASerGUC genes.

Anticodon↗

Evolution of the SNF2 family of proteins: subfamilies with distinct sequences and functions.

The SNF2 family of proteins includes representatives from a variety of species with roles in cellular processes such as transcriptional regulation (e.g. MOT1, SNF2 and BRM), maintenance of chromosome stability during mitosis (e.g. lodestar) and various aspects of processing of DNA damage, including nucleotide excision repair (e.g. RAD16 and ERCC6), recombinational pathways (e.g. RAD54) and post-replication daughter strand gap repair (e.g. RAD5). This family also includes many proteins with no known function. To better characterize this family of proteins we have used molecular phylogenetic techniques to infer evolutionary relationships among the family members. We have divided the SNF2 family into multiple subfamilies, each of which represents what we propose to be a functionally and evolutionarily distinct group. We have then used the subfamily structure to predict the functions of some of the uncharacterized proteins in the SNF2 family. We discuss possible implications of this evolutionary analysis on the general properties and evolution of the SNF2 family.

Amino Acid Sequence↗

The integrase family of tyrosine recombinases: evolution of a conserved active site domain.

The integrases are a diverse family of tyrosine recombinases which rearrange DNA duplexes by means of conservative site-specific recombination reactions. Members of this family, of which the well-studied lambda Int protein is the prototype, were previously found to share four strongly conserved residues, including an active site tyrosine directly involved in transesterification. However, few additional sequence similarities were found in the original group of 27 proteins. We have now identified a total of 81 members of the integrase family deposited in the databases. Alignment and comparisons of these sequences combined with an evolutionary analysis aided in identifying broader sequence similarities and clarifying the possible functions of these conserved residues. This analysis showed that members of the family aggregate into subfamilies which are consistent with their biological roles; these subfamilies have significant levels of sequence similarity beyond the four residues previously identified. It was also possible to map the location of conserved residues onto the available crystal structures; most of the conserved residues cluster in the predicted active site cleft. In addition, these results offer clues into an apparent discrepancy between the mechanisms of different subfamilies of integrases.

Amino Acid Sequence↗

Crystal structure of Bacillus subtilis TrmB, the tRNA (m7G46) methyltransferase.

The structure of Bacillus subtilis TrmB (BsTrmB), the tRNA (m7G46) methyltransferase, was determined at a resolution of 2.1 A. This is the first structure of a member of the TrmB family to be determined by X-ray crystallography. It reveals a unique variant of the Rossmann-fold methyltransferase (RFM) structure, with the N-terminal helix folded on the opposite site of the catalytic domain. The architecture of the active site and a computational docking model of BsTrmB in complex with the methyl group donor S-adenosyl-L-methionine and the tRNA substrate provide an explanation for results from mutagenesis studies of an orthologous enzyme from Escherichia coli (EcTrmB). However, unlike EcTrmB, BsTrmB is shown here to be dimeric both in the crystal and in solution. The dimer interface has a hydrophobic core and buries a potassium ion and five water molecules. The evolutionary analysis of the putative interface residues in the TrmB family suggests that homodimerization may be a specific feature of TrmBs from Bacilli, which may represent an early stage of evolution to an obligatory dimer.

Bacillus subtilis↗

Positive and negative selection in the DAZ gene family.

Because a microdeletion containing the DAZ gene is the most frequently observed deletion in infertile men, the DAZ gene was considered a strong candidate for the azoospermia factor. A recent evolutionary analysis, however, suggested that DAZ was free from functional constraints and consequently played little or no role in human spermatogenesis. The major evidence for this surprising conclusion is that the nonsynonymous substitution rate is similar to the synonymous rate and to the rate in introns. In this study, we reexamined the evolution of the DAZ gene family by using maximum-likelihood methods, which accommodate variable selective pressures among sites or among branches. The results suggest that DAZ is not free from functional constraints. Most amino acids in DAZ are under strong selective constraint, while a few sites are under diversifying selection with nonsynonymous/ synonymous rate ratios (d(N)/d(S)) well above 1. As a result, the average d(N)/d(S) ratio over sites is not a sensible measure of selective pressure on the protein. Lineage-specific analysis indicated that human members of this gene family were evolving by positive Darwinian selection, although the evidence was not strong.

Amino Acid Substitution↗