Neutralism and selectionism: the molecular clock

Identification, evolution, and expression of GDSL-type Esterase/Lipase (GELP) gene family in three cotton species: a bioinformatic analysis

BACKGROUND

GDSL esterase/lipases (GELPs) play important roles in plant growth, development, and response to biotic and abiotic stresses. Presently, an extensive and in-depth analysis of GELP family genes in cotton is still not clear enough, which greatly limits the further understanding of cotton GELP function and regulatory mechanism.

RESULTS

A total of 389 GELP family genes were identified in three cotton species of Gossypium hirsutum (193), G. arboreum (97), and G. raimondii (99). These GELPs could be classified into three groups and eight subgroups, with the GELPs in same group to have similar gene structures and conserved motifs. Evolutionary event analysis showed that the GELP family genes tend to be diversified at the spatial dimension and certain conservative at the time dimension, with a trend of potential continuous expansion in the future. The orthologous or paralogous GELPs among different genomes/subgenomes indicated the inheritance from genome-wide duplication during polyploidization, and the paralogous GELPs were derived from chromosomal segment duplication or tandem replication. GELP genes in the A/D subgenome underwent at least three large-scale replication events in the evolutionary process during the period of 0.6-3.2 MYA, with two large-scale evolutionary events between 0.6-1.8 MYA that were associated with tetraploidization, and the large-scale duplication between 2.6-9.1 MYA that occurred during diploidization. The cotton GELPs indicated diverse expression patterns in tissue development, ovule and fiber growth, and in response to biotic and abiotic stresses, combining the existing cis-elements in the promoter regions, suggesting the GELPs involvements of functions to be diversification and of the mechanisms to be a hormone-mediated manner.

CONCLUSIONS

Our results provide a systematic and comprehensive understanding the function and regulatory mechanism of cotton GELP family, and offer an effective reference for in-depth genetic improvement utilization of cotton GELPs.

Uncovering adaptive evolution in the human lineage

BACKGROUND

The recent increase in human polymorphism data, together with the availability of genome sequences from several primate species, provides an unprecedented opportunity to investigate how natural selection has shaped human evolution.

RESULTS

We compared human branch-specific substitutions with variation data in the current human population to measure the impact of adaptive evolution on human protein coding genes. The use of single nucleotide polymorphisms (SNPs) with high derived allele frequencies (DAFs) minimized the influence of segregating slightly deleterious mutations and improved the estimation of the number of adaptive sites. Using DAF ≥ 60% we showed that the proportion of adaptive substitutions is 0.2% in the complete gene set. However, the percentage rose to 40% when we focused on genes that are specifically accelerated in the human branch with respect to the chimpanzee branch, or on genes that show signatures of adaptive selection at the codon level by the maximum likelihood based branch-site test. In general, neural genes are enriched in positive selection signatures. Genes with multiple lines of evidence of positive selection include taxilin beta, which is involved in motor nerve regeneration and syntabulin, and is required for the formation of new presynaptic boutons.

CONCLUSIONS

We combined several methods to detect adaptive evolution in human coding sequences at a genome-wide level. The use of variation data, in addition to sequence divergence information, uncovered previously undetected positive selection signatures in neural genes.

Universal pacemaker of genome evolution in animals and fungi and variation of evolutionary rates in diverse organisms

Gene evolution is traditionally considered within the framework of the molecular clock (MC) model whereby each gene is characterized by an approximately constant rate of evolution. Recent comparative analysis of numerous phylogenies of prokaryotic genes has shown that a different model of evolution, denoted the Universal PaceMaker (UPM), which postulates conservation of relative, rather than absolute evolutionary rates, yields a better fit to the phylogenetic data. Here, we show that the UPM model is a better fit than the MC for genome wide sets of phylogenetic trees from six species of Drosophila and nine species of yeast, with extremely high statistical significance. Unlike the prokaryotic phylogenies that include distant organisms and multiple horizontal gene transfers, these are simple data sets that cover groups of closely related organisms and consist of gene trees with the same topology as the species tree. The results indicate that both lineage-specific and gene-specific rates are important in genome evolution but the lineage-specific contribution is greater. Similar to the MC, the gene evolution rates under the UPM are strongly overdispersed, approximately 2-fold compared with the expectation from sampling error alone. However, we show that neither Drosophila nor yeast genes form distinct clusters in the tree space. Thus, the gene-specific deviations from the UPM, although substantial, are uncorrelated and most likely depend on selective factors that are largely unique to individual genes. Thus, the UPM appears to be a key feature of genome evolution across the history of cellular life.

Universal pacemaker of genome evolution

A fundamental observation of comparative genomics is that the distribution of evolution rates across the complete sets of orthologous genes in pairs of related genomes remains virtually unchanged throughout the evolution of life, from bacteria to mammals. The most straightforward explanation for the conservation of this distribution appears to be that the relative evolution rates of all genes remain nearly constant, or in other words, that evolutionary rates of different genes are strongly correlated within each evolving genome. This correlation could be explained by a model that we denoted Universal PaceMaker (UPM) of genome evolution. The UPM model posits that the rate of evolution changes synchronously across genome-wide sets of genes in all evolving lineages. Alternatively, however, the correlation between the evolutionary rates of genes could be a simple consequence of molecular clock (MC). We sought to differentiate between the MC and UPM models by fitting thousands of phylogenetic trees for bacterial and archaeal genes to supertrees that reflect the dominant trend of vertical descent in the evolution of archaea and bacteria and that were constrained according to the two models. The goodness of fit for the UPM model was better than the fit for the MC model, with overwhelming statistical significance, although similarly to the MC, the UPM is strongly overdispersed. Thus, the results of this analysis reveal a universal, genome-wide pacemaker of evolution that could have been in operation throughout the history of life.

Molecular cloning, characterization of copper/zinc superoxide dismutase and expression analysis of stress-responsive genes from Eisenia fetida against dietary zinc oxide

The full length cDNA of copper/zinc superoxide dismutase (Cu/Zn-SOD) from Eisenia fetida (E. fetida) was cloned (GenBank accession no. JN579648). Sequence characterization revealed that the cDNA contained characteristic Cu/Zn-SOD family signatures ((45)GFHVHEFGDNT(55) and (138)GNAGGRLACGVI(149)), cysteines (Cys-58 and-146) predicted to form one disulphide bond, Cu-binding (His-47, -49, -64 and -120) and Zn-binding (His-64, -72, -81 and Asp-84). They were essential for the structure and function of Cu/Zn-SOD. Differential expression of stress-responsive genes like Cu/Zn-SOD, catalase (CAT), heat shock protein 70 (Hsp70) and metallothionein (MT) was applied as potential biomarkers to assess their efficacy for the ecotoxicological effects of dietary zinc oxide (ZnO) on E. fetida. The results showed that the expression of Cu/Zn-SOD and MT increased to reach the highest levels of 6.22 and 7.68 fold in a dose-dependent manner at day 10 respectively. The highest expression of 3.03 fold of CAT was registered at day 10. The transient expression of Hsp70 without consistent time- or/and dose-dependent was observed. It implied that the transcriptional patterns of Cu/Zn-SOD, CAT and MT could serve as early warning signals in ecotoxicological assessment of dietary ZnO on earthworms while the expression of Hsp70 was not well done, which is helpful to monitoring and regulation of ZnO in veterinary application.

Molecular cloning, expression, and immunolocalization of the NAD(+)-dependent glycerol 3-phosphate dehydrogenase (GPD) from Clonorchis sinensis

Glycerol 3-phosphate dehydrogenase (GPD) plays an important role in the energy metabolism and nutrition metabolism. In order to know about the biological functions of GPD of Clonorchis sinensis (C. sinensis), we identified a complete gene coding GPD from C. sinensis metacercaria cDNA library. This novel cDNA sequence contains 1,056 bp with a putative open reading frame of 351 amino acids and shares 74% identity with GPD from Schistosoma mansoni. Recombinant CsGPD was expressed and purified from Escherichia coli BL21 (DE3). Western blot analysis displayed that recombinant CsGPD can be recognized by anti-CsGPD serum and C. sinensis-infected serum. RT-PCR and immunolocalization analysis confirmed that GPD expressed both at the stage of adult worm and metacercaria of C. sinensis and immunolocated at the tegument of adult worm, tegument and tegumentary cells of metacercaria. Our current study has paved the way for the further researches about the biological functions involved in the growth of C. sinensis.

Sequence variability and sequence evolution: an explanation of molecular polymorphisms and why many molecular structures can be preserved although they are not predominant

The existence of many processes that regulate RNA expression poses a challenge to the idea that the cell is the culmination of a highly efficient interplay of individual proteins, each with specific, highly specialized functions. It will be demonstrated here the extent to which the cell may undergo evolutionary processes that also occur in the macrocosmos, specifically with reference to the rules of mutation and preservation. These molecular evolutionary processes could facilitate a better understanding of the development of molecular structures and the functioning of the cell and could give an explanation of the molecular polymorphisms and also explain why many molecular structures can be preserved although they are not predominant.

Molecular characterization and mRNA expression during metal exposure and thermal stress of copper/zinc- and manganese-superoxide dismutases in disk abalone, Haliotis discus discus

Complementary DNAs encoding copper/zinc superoxide dismutase (Cu/Zn-SOD; SOD1) and manganese superoxide dismutase (Mn-SOD; SOD2) were isolated from disk abalone, Haliotis discus discus. The open reading frame sequences of Cu/Zn- and Mn-SODs encoded 154 and 226 amino acids, respectively. Multiple sequence alignments using the deduced amino acid sequences revealed that both abalone SODs showed considerable sequence similarities with their orthologues from diverse aerobic organisms, in which the amino acid residues forming metal ligands were highly conserved. All phylogenetic trees for both SOD genes inferred from maximum likelihood and Bayesian inference analyses presented the monophyletic status of Teleostei and Aves/Tetrapoda clades, and recovered relatively close genetic affiliation of H. discus discus with some molluscan species. Expression of both SODs at mRNA levels were highly modulated in various tissues (gill, muscle and hepatopancreas from juveniles, and haemocytes from adults) by experimental exposures to heavy metals (copper, zinc and cadmium) and also by thermal treatments (elevation of temperature). The mRNA levels of both SODs were increased in general during the metal or thermal treatments; however, the transcriptional responses of SOD genes were quite variable depending upon isoforms and tissues based on semi-quantitative and/or real-time RT-PCR assays.

Plant pathogen forensics: capabilities, needs, and recommendations

A biological attack on U.S. crops, rangelands, or forests could reduce yield and quality, erode consumer confidence, affect economic health and the environment, and possibly impact human nutrition and international relations. Preparedness for a crop bioterror event requires a strong national security plan that includes steps for microbial forensics and criminal attribution. However, U.S. crop producers, consultants, and agricultural scientists have traditionally focused primarily on strategies for prevention and management of diseases introduced naturally or unintentionally rather than on responding appropriately to an intentional pathogen introduction. We assess currently available information, technologies, and resources that were developed originally to ensure plant health but also could be utilized for postintroduction plant pathogen forensics. Recommendations for prioritization of efforts and resource expenditures needed to enhance our plant pathogen forensics capabilities are presented.

Molecular clock-like evolution of human immunodeficiency virus type 1

The molecular clock hypothesis states that the rate of nucleotide substitution per generation is constant across lineages. If generation times were equal across lineages, samples obtained at the same calendar time would have experienced the same number of generations since their common ancestor. However, if sequences are not derived from contemporaneous samples, differences in the number of generations may be misinterpreted as variation in substitution rates and hence may lead to false rejection of the molecular clock hypothesis. A recent study has called into doubt the validity of clock-like evolution for HIV-1, using molecular sequences derived from noncontemporaneous samples. However, after separating their within-individual data according to sampling time, we found that what appeared to be nonclock-like behavior could be attributed, in most cases, to noncontemporaneous sampling, with contributions also likely to derive from recombination. Natural selection alone did not appear to obscure the clock-like evolution of HIV-1.

Ancestral gene fusion in cellobiose dehydrogenases reflects a specific evolution of GMC oxidoreductases in fungi

Cellobiose dehydrogenases (CDHs) are extracellular hemoflavoenzymes that are thought to be involved in the degradation of two of the most abundant biopolymers in the biosphere, cellulose and lignin. To date, these enzymes, consisting of a cytochrome domain and a flavin domain, have been detected and sequenced exclusively in the kingdom of fungi. Independent phylogenetic analyses of two distinct domains of CDH genes reveal that they evolved in parallel as fused genes. Whereas the cytochrome domains are unique sequence motifs, the flavin domains clearly belong to the glucose-methanol-choline (GMC) oxidoreductase family--an evolution line of widespread flavoproteins extending from the Archae to higher eukaryotes. The most probable unrooted phylogenetic tree obtained from our analysis of 52 selected GMC members reveals five principal evolutionary branches: cellobiose dehydrogenase, cholesterol oxidase (COX), hydroxynitrile lyase, alcohol oxidase (AOX)/glucose oxidase (GOX)/choline dehydrogenase, and a branch of dehydrogenases with various specificities containing also an Archaeon open reading frame (ORF). Cellobiose dehydrogenases cluster with cholesterol oxidases and the clade of various specificities, whereas hydroxynitrile lyases are closely related to glucose oxidases, alcohol oxidases, and choline dehydrogenases. The results indicate that the evolutionary line from a primordial GMC flavoprotein to extant cellobiose dehydrogenases was augmented after an early acquisition of the cytochrome domain to form two distinct branches for basidiomycetes and ascomycetes. One ascomycetous evolutionary line of CDHs has acquired a carbohydrate-binding module (CBM) of type 1, the sequence of which is similar to that of corresponding domains in several glycosidases. This is the first attempt towards a comprehensive phylogenetic analysis of cellobiose dehydrogenases.

Molecular clock and gene function

Molecular phylogenies based on the molecular clock require the comparison of orthologous genes. Orthologous and paralogous genes usually have very different evolutionary fates. In general, orthologs keep the same functions in species, whereas, particularly over a long time span, paralogs diverge functionally and may become pseudogenes or get lost. In eukaryotic genomes, because of the degree of redundancy of genetic information, homologous genes are grouped in gene families, the evolution of which may differ greatly between the various organisms. This implies that each gene in a species does not always have an ortholog in another species and thus, due to multiple duplication events following a speciation, many orthologous clades of paralogs are generated. We are often dealing with a one-to-many or many-to-many relationship between genes. In this paper, we analyze the evolution of two gene families, the p53 gene family and the porin gene family. The evolution of the p53 family shows a one-to-many gene relationship going from invertebrates to vertebrates. In invertebrates only a single gene has been found, while in vertebrates three members of the family, namely p53, p63, and p73, are present. The evolution of porin (VDAC) genes (VDAC1, VDAC2, and VDAC3) is an example of a many-to-many gene relationship going from yeast to mammals. However, the porin gene redundancy found in invertebrates and possibly in some fishes may indicate a tendency to duplicate the genetic material, rather than a real need for function innovation.

The apparent clock-like evolution of Escherichia coli in glucose-limited chemostats is reproducible at large but not at small population sizes and can be explained with Monod kinetics

To follow and model evolution of a microbial population in the chemostat, parameters are needed that give an indication of the absolute extent of evolution at a high resolution of time. In this study the evolution of the maximum specific growth rate ( micro (max)) and the residual glucose concentration was followed for populations of Escherichia coli K-12 under glucose-limited conditions at dilution rates of 0.1 x h(-1), 0.3 x h(-1) and 0.53 x h(-1) during 500-700 h in continuous culture. Whereas micro (max) improved only during the initial 150 h, the residual glucose concentration decreased constantly during 500 h of cultivation and therefore served as a convenient parameter to monitor the evolution of a population at a high time resolution with respect to its affinity for the growth-limiting substrate. The evolution of residual glucose concentrations was reproducible in independent chemostats with a population size of 10(11) cells, whereas no reproducibility was found in chemostats containing 10(7) cells. A model based on Monod kinetics assuming successive take-overs of mutants with improved kinetic parameters (primarily K(s)) was able to simulate the experimentally observed evolution of residual glucose concentrations. Similar values for the increase in glucose affinity of mutant phenotypes (K(s(mutant)) approximately equal 0.6 x K(s(parent))) and similar mutation rates per cell per generation leading to these mutant phenotypes (1-5 x1 0(-7)) were estimated in silico for all dilution rates. The model predicts a maximum rate of evolution at a dilution rate slightly below micro (max)/2. With increasing and decreasing dilution rates the evolution slows down, which also explains why in special cases a selection-driven evolution can exhibit apparent clock-like behaviour. The glucose affinity for WT cells was dependent on the dilution rate with highest values at dilution rates around micro (max)/2. Below 0.3 x h(-1) poorer affinity was mainly due to the effects of rpoS.

Structure, molecular evolution, and gene expression of primate superoxide dismutases

Mn- and Cu,Zn-superoxide dismutase (SOD) cDNAs of eight primate species, Pan troglodytes, Pongo pygmaeus, Hylobates lar, Macaca fuscata, Macaca fascicularis, Macaca mulatta, Cebus apella, and Callithrix jacchus, were cloned. The whole protein-coding sequences were covered, comparing 198 and 153 (or 154) amino acids, for Mn- and Cu,Zn-SODs, respectively. Residues forming metal ligands were completely conserved in the two primate SODs and nucleotide/amino acid substitutions were more frequent in Cu,Zn-SODs than in Mn-SODs. Molecular evolutionary analyses showed Mn-SOD to have evolved at a constant rate and its phylogenetic tree well reflected primate phylogeny. Cu,Zn-SOD was shown to have evolved differently between primate lineages. The significant high ratio of a non-synonymous/synonymous rate was found in the lineage leading to great apes and humans, showing that this lineage underwent positive Darwinian selection. Southern hybridization suggested that the genes for primate Mn- and Cu,Zn-SOD exist as single copies. Northern analysis in various Japanese monkey tissues showed Mn- and Cu,Zn-SOD expression to be high in the liver, kidneys, and adrenal glands.

Mutator dynamics in fluctuating environments

Populations with high mutation rates (mutator clones) are being found in increasing numbers of species, and a clear link is being established between the presence of mutator clones and drug resistance. Mutator clones exist despite the fact that in a constant environment most mutations are deleterious, with the spontaneous mutation rate generally held at a low value. This implies that mutator clones have an important role in the adaptation of organisms to changing environments. Our study examines how mutator dynamics vary according to the frequency of environmental fluctuations. Although recent studies have considered a single environmental switch, here we investigate mutator dynamics in a regularly varying environment, seeking to mimic conditions present, for example, under certain drug or pesticide regimes. Our model provides four significant new insights. First, the results demonstrate that mutators are most prevalent under intermediate rates of environmental change. When the environment oscillates more rapidly, mutators are unable to provide sufficient adaptability to keep pace with the frequent changes in selection pressure and, instead, a population of 'environmental generalists' dominates. Second, our findings reveal that mutator dynamics may be complex, exhibiting limit cycles and chaos. Third, we demonstrate that when each beneficial mutation provides a greater gain in fitness, mutators achieve higher densities in more rapidly fluctuating environments. Fourth, we find that mutators of intermediate strength reach higher densities than very weak or strong mutators.

Erratic overdispersion of three molecular clocks: GPDH, SOD, and XDH

The neutrality theory predicts that the rate of neutral molecular evolution is constant over time, and thus that there is a molecular clock for timing evolutionary events. It has been observed that the variance of the rate of evolution is generally larger than expected according to the neutrality theory, which has raised the question of how reliable the molecular clock is or, indeed, whether there is a molecular clock at all. We have carried out an extensive investigation of three proteins, glycerol-3-phosphate dehydrogenase (GPDH), superoxide dismutase (SOD), and xanthine dehydrogenase (XDH). We have observed that (i) the three proteins evolve erratically through time and across lineages and (ii) the erratic patterns of acceleration and deceleration differ from locus to locus, so that one locus may evolve faster in one than another lineage, whereas the opposite may be the case for another locus. The observations are inconsistent with the predictions made by various subsidiary hypotheses proposed to account for the overdispersion of the molecular clock.