Biomedical subjects

Joseph P Bielawski

Publications and source records attributed to Joseph P Bielawski.

17 recordsLinked to original sources

PubMed2006

Large-scale analyses of synonymous substitution rates can be sensitive to assumptions about the process of mutation.

A popular approach to examine the roles of mutation and selection in the evolution of genomes has been to consider the relationship between codon bias and synonymous rates of molecular evolution. A significant relationship between these two quantities is taken to indicate the action of weak selection on substitutions among synonymous codons. The neutral theory predicts that the rate of evolution is inversely related to the level of functional constraint. Therefore, selection against the use of non-preferred codons among those coding for the same amino acid should result in lower rates of synonymous substitution as compared with sites not subject to such selection pressures. However, reliably measuring the extent of such a relationship is problematic, as estimates of synonymous rates are sensitive to our assumptions about the process of molecular evolution. Previous studies showed the importance of accounting for unequal codon frequencies, in particular when synonymous codon usage is highly biased. Yet, unequal codon frequencies can be modeled in different ways, making different assumptions about the mutation process. Here we conduct a simulation study to evaluate two different ways of modeling uneven codon frequencies and show that both model parameterizations can have a dramatic impact on rate estimates and affect biological conclusions about genome evolution. We reanalyze three large data sets to demonstrate the relevance of our results to empirical data analysis.

Amino Acid Substitution↗

PubMed2006

Evolutionary rate variation among vertebrate beta globin genes: implications for dating gene family duplication events.

A comprehensive dataset of 62 beta globin gene sequences from various vertebrates was compiled to test the molecular clock and to estimate dates of gene duplications. We found that evolution of the beta globin family of genes is not clock-like, a result that is at odds with the common use of this family as an example of a constant rate of evolution over time. Divergence dates were estimated either with or without assuming the molecular clock, and both analyses produced similar date estimates, which are also in general agreement with estimates reported previously. In addition we report date estimates for seven previously unexamined duplication events within the beta globin family. Despite multiple sources of rate variation, the average rate across the beta globin phylogeny yielded reasonable estimates of divergence dates in most cases. Exceptions were cases of gene conversion, where it appears to have led to underestimates of divergence dates. Our results suggest (i) the major duplications giving rise to the paralogous beta globin genes are associated with significant evolutionary rate variation among gene lineages; and (ii) genes arising from more recent gene duplications (e.g., tandem duplications within lineages) do not appear to differ greatly in rate. We believe this pattern reflects a complex interplay of evolutionary forces where natural selection for diversifying paralogous functions and lineage-specific effects contribute to rate variation on a long-term basis, while gene conversion tends to increase sequence similarity. Gene conversion effects appear to be stronger on recent gene duplicates, as their sequences are highly similar. Lastly, phylogenetic analyses do not support a previous report that avian globins are members of a relic lineage of omega globins.

Animals↗

PubMed2006

Adaptive evolution of multicolored fluorescent proteins in reef-building corals.

Here we investigate the evolutionary scenarios that led to the appearance of fluorescent color diversity in reef-building corals. We show that the mutations that have been responsible for the generation of new cyan and red phenotypes from the ancestral green were fixed with the help of positive natural selection. This fact strongly suggests that the color diversity is a product of adaptive evolution. An unexpected finding was a set of residues arranged as an intermolecular binding interface, which was also identified as a target of positive selection but is nevertheless not related to color diversification. We hypothesize that multicolored fluorescent proteins evolved as part of a mechanism regulating the relationships between the coral and its algal endosymbionts (zooxanthellae). We envision that the effect of the proteins' fluorescence on algal physiology may be achieved not only through photosynthesis modulation, but also through regulatory photosensors analogous to phytochromes and cryptochromes of higher plants. Such a regulation would require relatively subtle, but spectrally precise, modifications of the light field. Evolution of such a mechanism would explain both the adaptive diversification of colors and the coevolutionary chase at the putative algae-protein binding interface in coral fluorescent proteins.

Adaptation, Biological↗

PubMed2006

Proposed standard nomenclature for the alpha- and beta-globin gene families.

The globin family of genes and proteins has been a recurrent object of study for many decades. This interest has generated a vast amount of knowledge. However it has also created an inconsistent and confusing nomenclature, due to the lack of a systematic approach to naming genes and failure to reflect the phylogenetic relationships among genes of the gene family. To alleviate the problems with the existing system, here we propose a standardized nomenclature for the alpha and beta globin family of genes, based on a phylogenetic analysis of vertebrate alpha and beta globins, and following the Guidelines for Human Gene Nomenclature.

Animals↗

PubMed2006

Prevalence and evolution of core photosystem II genes in marine cyanobacterial viruses and their hosts.

Cyanophages (cyanobacterial viruses) are important agents of horizontal gene transfer among marine cyanobacteria, the numerically dominant photosynthetic organisms in the oceans. Some cyanophage genomes carry and express host-like photosynthesis genes, presumably to augment the host photosynthetic machinery during infection. To study the prevalence and evolutionary dynamics of this phenomenon, 33 cultured cyanophages of known family and host range and viral DNA from field samples were screened for the presence of two core photosystem reaction center genes, psbA and psbD. Combining this expanded dataset with published data for nine other cyanophages, we found that 88% of the phage genomes contain psbA, and 50% contain both psbA and psbD. The psbA gene was found in all myoviruses and Prochlorococcus podoviruses, but could not be amplified from Prochlorococcus siphoviruses or Synechococcus podoviruses. Nearly all of the phages that encoded both psbA and psbD had broad host ranges. We speculate that the presence or absence of psbA in a phage genome may be determined by the length of the latent period of infection. Whether it also carries psbD may reflect constraints on coupling of viral- and host-encoded PsbA-PsbD in the photosynthetic reaction center across divergent hosts. Phylogenetic clustering patterns of these genes from cultured phages suggest that whole genes have been transferred from host to phage in a discrete number of events over the course of evolution (four for psbA, and two for psbD), followed by horizontal and vertical transfer between cyanophages. Clustering patterns of psbA and psbD from Synechococcus cells were inconsistent with other molecular phylogenetic markers, suggesting genetic exchanges involving Synechococcus lineages. Signatures of intragenic recombination, detected within the cyanophage gene pool as well as between hosts and phages in both directions, support this hypothesis. The analysis of cyanophage psbA and psbD genes from field populations revealed significant sequence diversity, much of which is represented in our cultured isolates. Collectively, these findings show that photosynthesis genes are common in cyanophages and that significant genetic exchanges occur from host to phage, phage to host, and within the phage gene pool. This generates genetic diversity among the phage, which serves as a reservoir for their hosts, and in turn influences photosystem evolution.

Base Composition↗

PubMed2005

Evolutionary rates and expression level in Chlamydomonas.

In many biological systems, especially bacteria and unicellular eukaryotes, rates of synonymous and nonsynonymous nucleotide divergence are negatively correlated with the level of gene expression, a phenomenon that has been attributed to natural selection. Surprisingly, this relationship has not been examined in many important groups, including the unicellular model organism Chlamydomonas reinhardtii. Prior to this study, comparative data on protein-coding sequences from C. reinhardtii and its close noninterfertile relative C. incerta were very limited. We compiled and analyzed protein-coding sequences for 67 nuclear genes from these taxa; the sequences were mostly obtained from the C. reinhardtii EST database and our C. incerta EST data. Compositional and synonymous codon usage biases varied among genes within each species but were highly correlated between the orthologous genes of the two species. Relative rates of synonymous and nonsynonymous substitution across genes varied widely and showed a strong negative correlation with the level of gene expression estimated by the codon adaptation index. Our comparative analysis of substitution rates in introns of lowly and highly expressed genes suggests that natural selection has a larger contribution than mutation to the observed correlation between evolutionary rates and gene expression level in Chlamydomonas.

Animals↗

PubMed2005

Adaptive diversification of vomeronasal receptor 1 genes in rodents.

The vomeronasal receptor 1 (V1R) are believed to be pheromone receptors in rodents. Here we used computational methods to identify 95 and 62 new putative V1R genes from the draft rat and mouse genome sequence, respectively. The rat V1R repertoire consists of 11 subfamilies, 10 of which are shared with the mouse, while rat appears to lack the H and I subfamilies found in mouse and possesses one unique subfamily (M). The estimations of the relative divergence times suggest that many subfamilies originated after the split of rodents and primates. The analysis also reveals that these clusters underwent an expansion very close to the split of mouse and rat. In addition, maximum likelihood analysis showed that the nonsynonymous and synonymous rate ratio for most of these clusters was much higher than one, suggesting the role of positive selection in the diversification of these duplicated V1R genes. Because V1R are thought to mediate the process of signal transduction in response to pheromone detection, we speculate that the V1R genes have evolved under positive Darwinian selection to maintain the ability to discriminate between large and complex pheromonal mixtures.

Adaptation, Physiological↗

PubMed2005

Potential photosynthesis gene recombination between Prochlorococcus and Synechococcus via viral intermediates.

Genes (psbA and psbD) encoding for photosynthetically important proteins were recently found in a number of cultured cyanophage genomes. This phenomenon may be a beneficial trait to the viruses or their photosynthetic cyanobacterial hosts, or may represent an untapped pool of genes involved in the formation of the photosynthetic apparatus that are prone to lateral gene transfer. Here we show analyses of psbA genes from uncultured environmental viruses and prophage populations. We observe a statistically significant separation between viral genes and their potential Synechococcus hosts' genes, and statistical analyses under models of codon evolution indicate that the psbA genes of viruses are evolving under levels of purifying selection that are virtually indistinguishable from their hosts. Furthermore, our data also indicate the possible exchange and reshuffling of psbA genes between Synechococcus and Prochlorococcus via phage intermediates. Overall, these observations raise the possibility that marine viruses serve as a potential genetic pool in shaping the evolution of cyanobacterial photosynthesis.

Bacterial Proteins↗

PubMed2004

Darwinian adaptation of proteorhodopsin to different light intensities in the marine environment.

Proteorhodopsin, a retinal-binding protein, represents a potentially significant source of light-driven energy production in the world's oceans. The distribution of photochemically divergent proteorhodopsins is stratified according to depth. Here, we present evidence that such photochemical diversity was tuned by Darwinian selection. By using a Bayesian method, we identified sites targeted by Darwinian selection and mapped them to three-dimensional models of proteorhodopsins. We suggest that spectral fine-tuning results from the combined effect of amino acids that directly interact with retinal and those that influence the confirmation of the retinal-binding pocket.

Acclimatization↗

PubMed2004

A maximum likelihood method for detecting functional divergence at individual codon sites, with application to gene family evolution.

The tailoring of existing genetic systems to new uses is called genetic co-option. Mechanisms of genetic co-option have been difficult to study because of difficulties in identifying functionally important changes. One way to study genetic co-option in protein-coding genes is to identify those amino acid sites that have experienced changes in selective pressure following a genetic co-option event. In this paper we present a maximum likelihood method useful for measuring divergent selective pressures and identifying the amino acid sites affected by divergent selection. The method is based on a codon model of evolution and uses the nonsynonymous-to-synonymous rate ratio (omega) as a measure of selection on the protein, with omega = 1, < 1, and > 1 indicating neutral evolution, purifying selection, and positive selection, respectively. The model allows variation in omega among sites, with a fraction of sites evolving under divergent selective pressures. Divergent selection is indicated by different omega's between clades, such as between paralogous clades of a gene family. We applied the codon model to duplication followed by functional divergence of (i) the epsilon and gamma globin genes and (ii) the eosinophil cationic protein (ECP) and eosinophil-derived neurotoxin (EDN) genes. In both cases likelihood ratio tests suggested the presence of sites evolving under divergent selective pressures. Results of the epsilon and gamma globin analysis suggested that divergent selective pressures might be a consequence of a weakened relationship between fetal hemoglobin and 2,3-diphosphoglycerate. We suggest that empirical Bayesian identification of sites evolving under divergent selective pressures, combined with structural and functional information, can provide a valuable framework for identifying and studying mechanisms of genetic co-option. Limitations of the new method are discussed.

Amino Acids↗

PubMed2004

Gene conversion and functional divergence in the beta-globin gene family.

Different models of gene family evolution have been proposed to explain the mechanism whereby gene copies created by gene duplications are maintained and diverge in function. Ohta proposed a model which predicts a burst of nonsynonymous substitutions following gene duplication and the preservation of duplicates through positive selection. An alternative model, the duplication-degeneration-complementation (DDC) model, does not explicitly require the action of positive Darwinian selection for the maintenance of duplicated gene copies, although purifying selection is assumed to continue to act on both copies. A potential outcome of the DDC model is heterogeneity in purifying selection among the gene copies, due to partitioning of subfunctions which complement each other. By using the d(N)/ d(S) (omega) rate ratio to measure selection pressure, we can distinguish between these two very different evolutionary scenarios. In this study we investigated these scenarios in the beta-globin family of genes, a textbook example of evolution by gene duplication. We assembled a comprehensive dataset of 72 vertebrate beta-globin sequences. The estimated phylogeny suggested multiple gene duplication and gene conversion events. By using different programs to detect recombination, we confirmed several cases of gene conversion and detected two new cases. We tested evolutionary scenarios derived from Ohta's model and the DDC model by examining selective pressures along lineages in a phylogeny of beta-globin genes in eutherian mammals. We did not find significant evidence for an increase in the omega ratio following major duplication events in this family. However, one exception to this pattern was the duplication of gamma-globin in simian primates, after which a few sites were identified to be under positive selection. Overall, our results suggest that following gene duplications, paralogous copies of beta-globin genes evolved under a nonepisodic process of functional divergence.

Animals↗

PubMed2003

Widespread adaptive evolution in the human immunodeficiency virus type 1 genome.

We investigated variable selective pressures among amino acid sites in HIV-1 genes. Selective pressure at the amino acid level was measured by using the nonsynonymous/synonymous substitution rate ratio (omega = dN/dS). To identify amino acid sites under positive selection with omega > 1, we applied maximum likelihood models that allow variable omega ratios among sites to analyze genomic sequences of 26 HIV-1 lineages including subtypes A, B, and C. Likelihood ratio tests detected sites under positive selection in each of the major genes in the genome: env, gag, pol, vif, and vpr. Positive selection was also detected in nef, tat, and vpu, although those genes are very small. The majority of positive selection sites is located in gp160. Positive selection was not detected if omega was estimated as an average across all sites, indicating the lack of power of the averaging approach. Candidate positive selection sites were mapped onto the available protein tertiary structures and immunogenic epitopes. We measured the physiochemical properties of amino acids and found that those at positive selection sites were more diverse than those at variable sites. Furthermore, amino acid residues at exposed positive selection sites were more physiochemically diverse than at buried positive selection sites. Our results demonstrate genomewide diversifying selection acting on the HIV-1.

Amino Acid Sequence↗

PubMed2003

Novel Proteorhodopsin variants from the Mediterranean and Red Seas.

Proteorhodopsins, ubiquitous retinylidene photoactive proton pumps, were recently found in the widespread uncultured SAR86 bacterial group in oceanic surface waters. To survey proteorhodopsin diversity, new degenerate sets of proteorhodopsin primers were designed based on a genomic proteorhodopsin gene sequence originating from an Antarctic fosmid library. New proteorhodopsin variants were identified in Red Sea samples that were most similar to the original green-light absorbing proteorhodopsins found in Monterey Bay California. Unlike green-absorbing proteorhodopsins however, these new variants contained a glutamine residue at position 105, the same site recently shown to control spectral tuning in naturally occurring proteorhodopsins. Different proteorhodopsin variants were also found in the Mediterranean Sea. These proteorhodopsins formed new and distinctive proteorhodopsin groups. Phylogenetic analyses show that some of the new variants were very different from previously characterized proteorhodopsins, and formed the deepest branching groups found so far among marine proteorhodopsins. The existence of these varied proteorhodopsin sequences suggests that this class of proteins has undergone substantial evolution. These variants could represent functionally divergent paralogous genes, derived from the same or similar species, or orthologous proteorhodopsins that are distributed amongst divergent planktonic microbial taxa.

Amino Acid Sequence↗

PubMed2003

Maximum likelihood methods for detecting adaptive evolution after gene duplication.

The rapid accumulation of genomic sequences in public databases will finally allow large scale studies of gene family evolution, including evaluation of the role of positive Darwinian selection following a duplication event. This will be possible because recent statistical methods of comparing synonymous and nonsynonymous substitution rates permit reliable detection of positive selection at individual amino acid sites and along evolutionary lineages. Here, we summarize maximum-likelihood based methods, and present a framework for their application to analysis of gene families. Using these methods, we investigated the role of positive Darwinian selection in the ECP-EDN gene family of primates and the Troponin C gene family of vertebrates. We also comment on the limitations of these methods and discuss directions for further improvements.

Amino Acid Substitution↗

PubMed2002

Ancestral polymorphism and adaptive evolution in the trichothecene mycotoxin gene cluster of phytopathogenic Fusarium.

Filamentous fungi within the Fusarium graminearum species complex (Fg complex) are the primary etiological agents of Fusarium head blight (scab) of wheat and barley. Scab is an economically devastating plant disease that greatly limits grain yield and quality. In addition, scabby grain is often contaminated with trichothecene mycotoxins that act as virulence factors on some hosts, and pose a serious threat to animal health and food safety. Strain-specific differences in trichothecene metabolite profiles (chemotypes) are not well correlated with the Fg complex phylogeny based on genealogical concordance at six single-copy nuclear genes. To examine the basis for this discord between species and toxin evolution, a 19-kb region of the trichothecene gene cluster was sequenced in 39 strains chosen to represent the global genetic diversity of species in the Fg complex and four related species of Fusarium. Phylogenetic analyses demonstrated that polymorphism within these virulence-associated genes is transspecific and appears to have been maintained by balancing selection acting on chemotype differences that originated in the ancestor of this important group of plant pathogens. Chemotype-specific differences in selective constraint and evidence of adaptive evolution within trichothecene genes are also reported.

DNA, Fungal↗

PubMed2002

Mutation patterns of mitochondrial H- and L-strand DNA in closely related Cyprinid fishes.

Mitochondrial genome replication is asymmetric. Replication starts from the origin of heavy (H)-strand replication, displacing the parental H-strand as it proceeds along the molecule. The H-strand remains single stranded until light (L)-strand replication is initiated from a second origin of replication. It has been suggested that single-stranded H-strand DNA is more sensitive to mutational damage, giving rise to substitutional rate differences between the two strands and among genes in mammalian mitochondrial DNA. In this study, we analyzed sequences of the cytochrome b, ND4, ND4L, and COI genes of cyprinid fishes to investigate rates and patterns of nucleotide substitution in the mitochondrial genome. To test for strand-asymmetric mutation pressure, a likelihood-ratio test was developed and applied to the cyprinid sequences. Patterns of substitution and levels of strand-asymmetric mutation pressure were largely consistent with a mutation gradient between the H- and L-strand origins of replication. Significant strand bias was observed among rates of transitional substitution. However, biological interpretation of the direction and strength of strand asymmetry for specific classes of substitutions is problematic. The problem occurs because the rate of any single class of substitution inferred from one strand is actually a sum of rates on two strands. The validity of the likelihood-ratio test is not affected by this problem.

Animals↗

PubMed2002

Accuracy and power of bayes prediction of amino acid sites under positive selection.

Bayes prediction quantifies uncertainty by assigning posterior probabilities. It was used to identify amino acids in a protein under recurrent diversifying selection indicated by higher nonsynonymous (d(N)) than synonymous (d(S)) substitution rates or by omega = d(N)/d(S) > 1. Parameters were estimated by maximum likelihood under a codon substitution model that assumed several classes of sites with different omega ratios. The Bayes theorem was used to calculate the posterior probabilities of each site falling into these site classes. Here, we evaluate the performance of Bayes prediction of amino acids under positive selection by computer simulation. We measured the accuracy by the proportion of predicted sites that were truly under selection and the power by the proportion of true positively selected sites that were predicted by the method. The accuracy was slightly better for longer sequences, whereas the power was largely unaffected by the increase in sequence length. Both accuracy and power were higher for medium or highly diverged sequences than for similar sequences. We found that accuracy and power were unacceptably low when data contained only a few highly similar sequences. However, sampling a large number of lineages improved the performance substantially. Even for very similar sequences, accuracy and power can be high if over 100 taxa are used in the analysis. We make the following recommendations: (1) prediction of positive selection sites is not feasible for a few closely related sequences; (2) using a large number of lineages is the best way to improve the accuracy and power of the prediction; and (3) multiple models of heterogeneous selective pressures among sites should be applied in real data analysis.

Amino Acid Substitution↗